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I. INTRODUCTION 


A. GENERAL BACKGROUND 

The classical electronic countermeasure (ECM) techniques against incoming 
active missiles attack is to employ noise and deceptive jammers to capture the missile 
sensors range gate and angular gate and Steer it toward some false target generated by 
passive ECM devices. Such an ECM technique works well when the distance between 
missile and jammer is large and the jamming signal to missile return signal ratio is high. 
On the other hand, modern missiles are designed to counter classical ECM techniques, 
such as the use of a monopulse tracking technique which makes angular gate stealing 
very difficult, and the home-on-jam technique to home onto the jammer when the 
jamming signal is very high. As the missile closes in on the jammer, the jamming 
signal to missile return signal ratio decreases and the jamming effectiveness drops. By 
the time the distance between the missile and the jammer is reduced to less than 
_ 10 Km, there is even less chance that the jammer will achieve its mission. Hence, there 
is a need to design another form of ECM technique that is effective at the short range. 

One solution to the short range ECM problem is to use a high power narrow 
beam device, concentrating very high energy towards the missile at ranges up to 
10 Km. As most sensors are designed to have high sensitivity for target acquisition at 
far distances, they usually do not have very wide dynamic range for handling very high 
input power density. If the power density produced by the ECM device at 10 Km 
exceeds the breakdown linut of the missile sensor, the sensor may be damaged or 


malfunction, causing the missile to lose tracking ability. 


B. SYSTEM SPECIFICATIONS 

The Flyswatter Array is a high power, vehicle mounted, ECM system capable of 
delivering up to 5 GW of output peak power within a very narrow pencil beam. This 
creates a very high power density of the order of 10-20 W/cm? at a distance of 10 Km. 
This amount of energy when radiated onto a seeker head of an incoming missile will 
damage or degrade the highly sensitive sensor of the missile, thereby causing it to lose 


track of its target. The specifications of such a system are listed below : 


Frequency 10 GHz 
Total output peak power 5 GW 
Maximum aperture E-field 10.6 KV/cm 


Reflector size Slime 27 kOe mm 


Array size Cink 2.0m 

Array mounting plate size 3m x 4.5m 

Array type endfeed broadside array 
Element orientation length in parallel with the 


long side of the plate 


Full gain elevation coverage 0° to 45° 
Azimuth coverage 180° 

Pulse width 100ns 

Pulse repetition rate 10pps 
Required power density at 10 Km 10 - 20 W/cm? 
Overall structure (Scone Iola) 


The above specifications were derived from a simplified model consisting of uniformly 
distributed elements. More detailed description of such a simplified model will be 
described in Chapter IV. 


C. PROJECT PURPOSE 

Given that we have such a high power device, one can ask “what is the radiation 
intensity in the vicinity of the system? Will it cause a radiation hazard to the operators 
of the vehicle?” The purpose of this project is to simulate the ECM system using a 
more realistic TE,, mode slot array rather than uniformly distributed elements. The 
far field patterns are recorded to ensure that the power density requirement of 
10-20 W/cm? at 10 Km is met. Then the near field radiation patterns are plotted to 


determine the radiation hazard regions, if any. 
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Figure 1.1 Flyswatter Array Moun 


Il. RADIATION HAZARD EVALUATION PROCEDURE 


A. PROBLEM AREA IDENTIFICATION 
The project can be divided into two main areas. The first one is to determine the 
array configuration to meet system requirements and the second is to determine anv 
radiation hazard zones. 
1. Array Configuration Determination 
The system specifications listed in Chapter I, section B do not define the array 
configuration. The only restrictions are on the array type, element orientation and 


overall array dimension. Other parameters that must be chosen are as follows: 


e the dimension of element, 
e the spacing between elements, 
e the tvpes of field tapering and 


e the number of elements. 


By varying all these parameters, we hope to design a TE,, mode slot array 
that will meet the system requirements; namely, deliver up to 5 GW of output power to 
achieve 10-20 W/cm? of power density at 10 Km distance without exceeding the 
maximum aperture E-field breakdown limit of 10.63 KV/cm. At the same time, the 
sidelobes and backlobes must be Kept small so that the near field radiation level outside 
of the main beam is low enough not to cause any hazard to the operators of the 
system. 

2. Radiation Hazard Regions Determination 

After the array configuration has been defined and the system specifications 
are met, the near field patterns of the ECM system are then computed at various 
reflector tilting angles and array azimuthal pointing directions. The radiation hazard 


regions are then recorded. 


B. EVALUATION PROCEDURE 
The evaluation procedure used to accomplish the task is as follows : 
Step 1 : Compute the far field patterns of a simple plate mounted array with fewer 
elements and no reflector plate. The purpose is to determine the element 


spacing and element sizes so as to reduce unwanted grating lobes. 
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Step 2: Compute the far field and near field patterns of the array with the reflector 
plate, with different tapering types. This will determine the type of tapering 
that will give high gain, narrow main beam beamwidth and low sidelobe 
and backlobe levels. 

Step 3: Compute the far field patterns of the flyswatter array with full input power 
of 5 GW. By adjusting the number of elements and the type of tapering, 
one can deternune if the required power density at 10 Km distance is 
achieved. At the same time the maximum aperture E-field must also be 
checked to ensure that it is within the limit specified. 

Step 4: Compute the far field patterns of the array in the presence of the modeled 
vehicle and determine if the requirements are still met at various elevation 
and azimuthal pointing angles. Vary the number of elements and tapering 
type if necessary. 

Step 5: Compute the near field patterns of the array in the presence of the modeled 
vehicle at various elevation and azimuthal pointing angles and determine 


the radiation hazard zones. 


C. RADIATION EXPOSURE LIMITS 

The radiation exposure limits of radio frequency electromagnetic spectrum for 
humans used in this report is based on the recommendations made by the National 
Council on Radiation Protection and Measurements (NCRP) and American National 
Standards Institute (ANSI). The NCRP report [Ref. 1] suggested two different limits 
on the far field power density, averaged over 6 min. and 30 min. respectively as shown 
in Fig. 2.l(a). The average-over-6-min. curve is meant for occupational exposure 
criteria which is similar to the limits given in ANSI C95.1-1982 protection guide 
[Ref. 2]. The average-over-30-min. curve is one-fifth that of the occupational criterion 
and is the criterion for the general population. These limits apply to both continuous 
wave and pulsed signals. No specific peak power limitation was given. 

The ANSI C95.1-1982 protection guide suggested that the radiation intensity be 
nieasured at a distance of 5 cm or more from any object in the field and averaged over 
a 6 minute duration as shown in Fig. 2.1(b). Hence, it can be applied to both near 
field and far field radiation conditions. At 10 GHz, the power density limit 1s 
5 mW/cm7”, the same as that given in NRCP report. Therefore, we will use this as the 


occupational exposure limit. However, the Flyswatter Array has a duty cycle of 10°, 
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hence the linut on peak power density based on an average of 5 mW/cm* is equivalent 
ae 4 
to 5 KW,cmy or 77 dB W, m*. 
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" ——-— Average for 30 min 
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Frequency, MHz 


Equivalent . 
Frequency (/) ; (Magnetic 
aie pace (Electric field)? field)? 
MHz mW /cm? V?/m? A?/m? 
0.3-3 100 4x 10° 2.9 
3.0-30 900/f? 4 x 10° (900/f ?) 0.025(900/f ?) 
30-300 1.0 4x 10° 0.025 
300-1500 {/300 4 x 10° (f/300) 0.025 (//300) 
1500-100,000 5.0 2 10" O:125 


*From ANSI (1982). 
> measured 5 cm or greater from any object in the field and averaged for any 0.1 h (6 
min). 


(a) NCRP Proposed Limits. 
‘(Electric Field)?/1200x or 12z7(magnetic field)?, whichever is greater. 


(b) ANSI C95.1-1982 Protection Guides. 


Figure 2.1 Radiation Exposure Linuts. 
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Il. SOFTWARE SUPPORT 


The simulation of the Flyswatter Array and vehicle requires a software code that 
is capable of modeling the whole system structure and computing the near field and far 
field radiation patterns at X-band. The Basic Scattering Code of the Numerical 
Electromagnetic Code (NEC-BSC) series is chosen for its capabilities in computing 
radiation patterns of complex structures using the Geometrical Theory of Diffraction 
(GTD) method. This code provides the required calculation and analysis power and 
outputs the data in tabular form. Additional routines were then written to support the 
code. A brief description of each of these codes is included in this chapter. The names 
enclosed in the parenthesis are the file names of the codes. The flow charts and listings 
of the four supporting software codes are attached. in Appendix A. A detailed 
description of the Basic Scattering Code and its applications can be found in the 


references [Refs. 3.4]. 


A. BASIC SCATTERING CODE (GTD2P51) 


} 


1. Capabilities. 
e Frequency Range: 
* UHF and above. 
e Modeling: 
* Infinite Ground Plane : Perfect conductor or dielectric filled plate. 
* Plates : Perfect conductor plates or thin dielectric slabs. 
* Cylinders : Elliptical shape with slanted end caps. 
* Antenna : 
- Dzipole or slot. 
- Single element or array. 
- Uniform, piecewise sinusoid or TE), distribution. 
- Electric or magnetic excitation. 
a.) Receiver: 
- Dipole or slot. 
- Single element or array. 
- Uniform, piecewise sinusoid or TE), distribution. 


- Electric or magnetic excitation. 
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e Computation : 

* Far field : E-field in either @ or @ direction or versus frequency. 

* Near field : E-field, H-field and Power density computations with one of 
the following options : 
- Spherical coordinates. 
- Rectangular coordinates. 
- Frequency sweeping. 

* Antenna Coupling : Magnitude, Phase and dB values of reaction between 
the antennas with one of the following options : 
- Mutual impedance between the two antennas. 
- Ratio of the power-out over power-in using a modified Frii’s 

transmission formula for conjugate matched loads. 


- Maximum gains between the antennas using the Linville method. 


2. Limitations. 

e Plates should have edges at least a wavelength long. 

e Thin dielectric slabs should be at least a wavelength away from any source. 

: Incident fields should not strike the slab too close to grazing. 

e Each antenna element should be at least a wavelength from all edges. 

e The code delivers up to a second order answer, ie., it can include fields for ray 
paths that interact with at most two scattering centers. Hence, the plots are 
accurate for the first 30 dB only. 

e Double diffracted fields are not included. 

In many cases, the one wavelength limits may be reduced to a quarter wavelength for 


engineering purposes. 


3. Modifications. 
The GTD2P51 code has been modified to suit our requirements for this 
project. The modifications are carried out largely in the main program and subroutine 
OUTPUT which are listed below: } 


e Added a CT card to read in the title of input data set. This title will be printed 
out together with the data set and passed over to the plotting routine to be 


included as the heading of the plots. 
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e Modify the subroutine OUTPUT to print out the title, abbreviated codes and 
output data at specific formats for output unit No. 8. The abbreviated code 
consists of 5 characters : 

« Character 1: F, N for Far field or Near field. 

“ Character 2: F, 3S for Fixed or Scanned trequenc - 

* Character 3: R, S for Rectangular or Spherical coordinates. 
* Character 4: E, H, P for E-field, H-field or Power density. 

Character 5: 1, 2, 3 defines the variable parameters; 

- X,Y or Z for rectangular coordinates, 

- R, 0 or @ for spherical coordinates. 

e Maximum dimensions for plates, cylinders and sources have been increased as 
follows : 


* Maximum number of plates = 50. 


* Maximum number of cvlinders = 100. 


* Maximum number of sources = 10100. 


B. INTERFACING CODE (GIDRWSEL) 

The output from the Basic Scattering Code consists of a string of abbreviated 
control codes followed by several columns of data. The interfacing code (GTD Read, 
Write, SELect) reads in these abbreviated codes and converts them into useful 
information. It then displays this information for the user to select the type of plots 
that he wants. These selections include, 


° The tvpe of field, ie. E-field, H-field or power density in either X, Y, Z or R, 
9, @ directions or the total field. 
’ The type of plots, ie. magnitude plot or phase plot, linear or dB scale for Y 


axis values, linear or log plot. 


After the selection, the GFDRWSEL program will then generate another set of 
control codes to be used by the plotting program. These codes and the labels for each 
plot are placed in a control file using output unit No. 3. The program also chooses the 
correct data sets and places them in the output data file using output unit No. 2. 
These two files will be used by the plotting program (GTDPLOT). 


C. OUTPUT DATA PLOTTING CODE (GTDPLOT) 

This program is modified from an existing program called PLOTNF (Near Field 
Plot). Additional features were added to allow for all options in the interfacing code 
(GIDRWSEL). It reads in the control codes from the control file and obtains the 
data from output data file and plots them using DISSPLA subroutines. 


D. GEOMETRIES PLOTTING CODE (GIDGEOMP) 

This program reads in the same set of input data that is to be fed into the NEC- 
BSC code. It then selects only the geometrical data and plots out its 3-D geometries 
using DISSPLA subroutines. It is used as an aid to check the geometrical data of the 
input data set. 

The GIDGEOMP code selects the geometrical data and plots the following 


items: 


e Infinite ground plane, 
e Plates, 
e Cylinders and 


® Sources. 
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It also reads in and executes the following commands that are used in GTD2P51 : 


at Print title. 

e NC Remove all cylinders. 

e NG _ Remove ground plane. | 
e NP Remove all plates. 

e NS Remove all sources. 

e NX _ Start next set of data. 


e RT _ Rotate and translate coordinates. 


The user may choose to view the object from any viewing angle in order to give a 
better perspective of the structure. This code, however, does not have the capability to 


blank out hidden lines. 


E. SA CARD GENERATING CODE (GTDSA) 

The code SA stands for Source-Array. This is an interactive program where the 
user is asked to enter a series of data regarding the geometries of the array, the type of 
element used, the tapering type and the total input power. It then computes and 
generates a data file in the format specified by the NEC-BSC code for the SA card. 
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This program is capable of generating a slot antenna array with a rectangular 
grid arranged on the X,, Z, plane, as used in NEC-BSC code. All elements of the array 
are assumed to have the same dimensions and orientation. It can simulate either a 
TE,g mode or Uniform distribution for each element. There are 5 different tapering 
types to choose from, for each of the directions (along X, or Z, axis). Namely, 
Uniform , Cosine , Cosine-square , Triangular and Binomial . One can also choose to 
add a pedestal height to the tapering. Note that even without any pedestal (ie. 
pedestal height equal zero), the outer most elements will still have some excitation 
because the algorithm always includes one more element at each side of the array, 
having the same spacing as the rest, to be the one with Zero excitation. 

The distribution used in each element is assumed to be perfect so that the 


relationship between the E-field and the output power for TE,, mode excitation is 


abE* 
dn 





A 


where E is the maximum E-field across the width b of the element, a is the length of 
the element and n is the intrinsic wave impedence. For an array of N elements with 


total output power P the above equation becomes 


total’ 


_ ab 5 
Protal = 4n DE, ’ 


with E, being the maximum E-field of the ae 


element. Knowing the output power, 
element dimensions and the tapering tvpes, the excitation for each element E? can 
be computed. 

For a Uniform distribution, the relationship between the input power and the 


excitations is given below : 


ab a 
= 2 
Protal = " y Ey 


Note that given the same input power, element dimensions, number of elements and 
tapering types, the element excitation of each element for uniform distribution 1s only 
one half of that of the TE) mode. Hence, a Uniform distribution will have much 


lower aperture E-field than TE,, mode. 


tJ 
to 


NALYSIS AND RADIATION HAZARD 
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IV. ARRAY CONEIGURAIIG! 


The evaluation of radiation hazard is carried out according to the 5 steps spelled 


out in Chapter [I under Evaluation Procedure. 


A. ELEMENT SPACING AND SIZE VARIATION 

In order not to complicate the evaluation, we assume that the elements are 
arranged on a rectangular grid and all elements are of the same rectangular size and 
Orientation as shown in Fig. 4.1. The axes X, and Z, orientate the array. (The array 


has no fixed relationship with the global co-ordinates X, Y and Z until specified.) 


E field 


element spacing 
tn Zt afrection 





Figure 4.1 Plate mounted rectangular grid array. 
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There are three main factors that impose restrictions on the element spacing and 
sizes. Namely, supporting the TE,, mode only, suppressing grating lobes and sidelobes 
and reducing maximum aperture E-field. 

1. TE10 mode Propagation 

In order to support the TE,, mode, the waveguide elements must have a 
length of at least half a wavelength (A/2). Similarly, for the TE), mode to propagate, 
the width of the element must be at least half a wavelength, too. Hence, we have to 
restrict the width to less than half a wavelength. As the length increases, and exceeds 
one wavelength, the TE,, mode will begin to propagate. Therefore, in order to support 


only the TE,, mode, the element size must be restricted to the following values : 


1.0X > element length (a) 2 0.54 
0.54 > element width (b) 


2. Grating Lobes and Sidelobes Suppression 
Grating lobes will appear whenever the spacing between elements is greater 
than or equal to one wavelength. This can be seen in Fig. Bl of Appendix B. These 
three plots show the far field patterns of a simple plate mounted array with fixed 
spacing of 1.0A, but different element lengths of 0.54, 0.754 and ‘1.04. There is no 


reflector plate used and the other parameters that are common to all three plots are as 


follows : 
Array mounting plate size 4.5m x 3m 
Number of elements in X, direction 10 
Number of elements in Z, direction 10 
Element width 0.25% 
Amplitude Tapering in X, direction Uniform 
Amplitude Tapering in Z, direction Uniform 
Output power 100 W 
Frequency 10 GHz 
Distance 10 Km 


The array is placed on the X-Y plane with the X, axis parallel to the Y-axis. 
These figures show that grating lobes always occur at @ = 90° and 270° 
irrespective of the element length. The main difference between the three plots is the 


gain at 8 = 0°. It is highest for element length = 1.0A and reduces as the length 
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decreases. Hence, it may be desired to have the length of elements as large as possible 
to ensure higher gain, but not too large to create grating lobes. 

Fig. B2 shows the effects when the spacing between elements in the Z, 
direction is reduced. Fig. B2(a) and Fig. B2(b) are for element spacing of 0.754 while 
Fig. B2(c) shows the radiation pattern plot when the spacing is further reduced to 0.5). 
Note that at 0.54 spacing, the antenna gain pattern at 8 = 90° and 270° is totally 
cancelled as one would expect. It is very important to ensure that the radiation level 
near the region where 8 = 90° and 270° be kept to its minimum because these 
directions could be in line with the operator. 

Similar deductions can be drawn for the width of the element and the spacing 
in X, direction. Fig. B3 shows 3 plots of fixed element length and spacing in both the 
X, and Z, directions, but with different element widths of 0.5, 0.3754 and 0.2524. The 
wave in the element is still propagating in the TE,, mode and hence the E-field 
strength is uniform in the X, direction. This results in broader beamwidth and even at 
a spacing of 0.75A in the X, direction, the grating lobes appear to be quite significant. 
As the spacing in the X, direction is reduced to 0.5, as shown in Fig. B4, the grating 
lobes disappear. Further reduction in the spacing causes the radiation level at 8 = 90° 
to increase as shown in Fig. B5. Hence, the spacing in the X, direction should be fixed 
at 0.5A. 

3. Maximum Aperture E-field Radiation 
The equation below shows the relationship between the input power and the 


maximum E-field of the array, assuming a perfect TE,, mode of propagation. 


abE? 
maxX 
P otal =| ————___ 5) Eee (n) 
4Zy n=1 


From this equation, it is obvious that by increasing the element size a and b, the 


maximum E-field Ey will be reduced since E (n) is only a function of N and the 


ratio 
tapering type. Hence, by making the elements as large as possible, without affecting 
the mode of propagation and the grating lobes, the maximum aperture E-field mav be 
reduced. 
4. Effects of varying the element spacing and size 
From the above investigation, the effects due to the variation of the element 


spacing and size are summarized below. 


° Spacing between elements in the Z, direction should be less than one 
wavelength in order to reduce grating lobes. It may be desired to operate 
with spacing equal to 0.54 to reduce the radiation at regions near 
@ = 90° and 270°. 

° Spacing between elements in X, direction must be 0.5 or less to reduce 
grating lobes. At 0.54 spacing, the radiation level at 8 2 90° is very low 
and is preferred. 

° Element size affects the gain of the main lobe as well as the maximum power 
flux at array output for a fixed number of elements. The larger the element 


size, the higher the gain and the lower the maximum power flux. 


From the above observations, one realizes that the spacing between elements 
in both Z, and X, directions are restricted to 1.0A and 0.5A respectively in order to 
support only the TE,, mode of propagation and at the same time suppress the grating 
lobes. Hence, by operating with minimum gaps between elements in both directions so 
that the element size is as large as possible, one can achieve high array gain at low 
aperture E-field strength. Alternately, if radiation hazard is a more critical issue, then 
one would prefer to operate with element size of 0.54 x 0.254 and with element spacing 
of 0.5% in both directions to ensure minimum direct radiation for 9 between 90° and 
270°. 

Before we conclude this section, let us consider the case where the array area 
is fixed and we want to fill the area with as many elements as possible so that the peak 
excitation is reduced. Fig. B6 to B9 show that a fixed array area of 7.54 x 3.75A is to 
be filled by elements of length of either 0.75 or 0.54 and width of either 0.3754 or 
0.254. The plots in Fig. B6 and B7 are for Uniform tapering and they show identical 
radiation pattern in both Eg and EO cuts. Even the maximum excitations are all equal 
to 0.2442 x 10+ V/m with 100W of input power. The plots in Fig. B8 and B9 are for 
Cosine tapering and again the radiation patterns are almost identical except for some 
variation in the sidelobe level. However, the maximum excitation varies with the plots, 
with 0.4349 x 107 V/m for 0.75 x 0.375A, 0.4462 x 10* V/m for both 0.75\ x 0.252 
and 0.54 x 0.3754 and 0.4578 x 10* V/m for 0.54 x 0.254. Hence, the peak excitation 
for Cosine tapering increases as the element size is reduced, while that of the Uniform 


tapering remains unchanged. 
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B. TAPERING TYPES VARIATION 


Five commonly used tapering types are selected for comparison. These are: 


e Uniform tapering 

e Cosine tapering 

¢ Cosine-Square tapering 
e Triangular tapering and 


e Binomial tapering 


Fig. B10 to B14 show the far field and near field patterns of the array with 


different taperings for each figure. The tapering in X, and Z, directions for each plot is 


the same to simplify the evaluation. The array is placed on the X-Y plane with the X 


aot 


axis parallel to the Y-axis and the reflector plate is tilted at 45° to give maximum 


radiation pointing towards 8 = 90°. The other parameters chosen are as follows: - 


Reflector plate (single plate) 

Array Mounting plate 

Number of elements in X, direction 
Number of elements in Z, direction 
Element length 

Element width 

Element spacing in X, direction 
Element spacing in Z, direction 
Output power 


Frequency 


Fig 4.2 shows the geometry of the array. 


1. Far Field Pattern Comparison 


5.1m x 2.7m 


4.5m x 3m 


For far field radiation patterns, we are interested in the gain and beamwidth. 


Hence, the plots are taken with +45° coverage about the boresight for cuts in both 


elevation and azimuth. The origin of the coordinate system is located at the center of 


the reflector plate and the plots are computed at a distance of 10 Km. Table 1 shows 


the comparison of gain and beamwidth. 


From Table 1, we have the following relationship : 


Gain: Uniform > Cosine > Triangular > Cosine-Square > Binomial 


Beamwidth : Binomial < Cosine-Square < Triangular < Cosine < Uniform 


2/ 





Figure 4.2 Simplified Array Geometry. 


Note that the reduction in gain and the increase in beamwidth for Binomual tapering as 
compared to the others are quite significant. Hence the Binomial tapering may not be 
suitable. 

2. Near Field Pattern Comparison 

For near field comparisons, we are interested in radiation levels in all 
directions. Hence, 0° to 360° coverage in both vertical and horizontal cuts are 
obtained. These patterns are computed at 5m distance from the center of the reflector 
plate. 

The plots show that with Uniform tapering, the sidelobe levels are very high 
and they occupy quite a wide area, hence, they are likely to cause a radiation hazard 
in a Wider region. With Cosine tapering, the mainlobe beamwidth 1s slightly broadened 
but the sidelobe level is reduced by about 10 dB. The mainlobe beamwidth is even 


wider for Cosine-Square and Triangular tapering. However, the sidelobe level is also 
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TABLE 1 
COMPARISON OF GAIN AND BEAMWIDTH 


Taping . . ne 
Uniform Cosine Cosine Triangular} Binomial 
Parameter 


Max. E-field 
“27.15 -28.29 -29.90 -28.95 -32.47 
in dB V/m 
BW in 
Se ea) cut 10.25 1250 15,50 1351.0 0 21.00 
BW in 
horizontal 10.00 127250 15.00 sa 5 0 20.00 
cut (deg) 





proportionally reduced. As for Binomial tapering, the mainlobe beamwidth is the 
broadest and theoretically, it should offer the lowest sidelobe level among the 5 
tapering types. But due to the broad beamwidth and the limited reflector plate size, 
there is spillover at 8 near 40°, making the sidelobe level appear to be much higher 
than for Cosine-Square and Triangular tapering. However, as the number of elements 
increases, the beamwidth will be reduced and Binomial tapering should offer the lowest 
sidelobe level. 

Notice that in the vertical cut near field plots, there is a considerable amount 
of radiation at @ near to 180° and 300° for Uniform, Cosine and Triangular tapering. 
These additional sidelobes also appear in the far field pattern but are not shown in our 
plots due to the limited angular coverage. The first lobe near 180° is part of the 
radiation pattern continuing from the main lobe. The sudden drop in radiation level 
prior to this lobe is due to the presence of the array mounting plate. By simply 
increasing the plate size, this lobe can be eliminated. The second lobe near 9 = 300° is 
the result of one of the sidelobes generated by the array itself. The presence of the 
reflector causes distortion on this sidelobe. For Cosine-Square and Binomial tapering, 
the sidelobe levels of the array are so low that no significant sidelobes appear at the 


near field patterns beyond 8 = 180°. 
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3. Preferred Tapering Type 
The far field pattern comparison shows that Uniform tapering gives the 
highest gain and narrowest mainlobe beamwidth. On the other hand, the near field 
pattern comparison suggests that Cosine-Square or Binomial tapering offer the lowest 
sidelobe level. A compromise between gain and sidelobe levels would make Cosine- 


Square tapering the preferred choice. 


C. OVERALL ARRAY CONFIGURATION ANALYSIS 

In this section, we will attempt to determine the number of elements in both the 
X, and is directions, element sizes, as well as the tapering type in order to meet the far 
field requirement of 10-20 W/cm? of power density at 10 Km, subject to the following 


constraints as spelled out in the system specifications. 


e Maximum total output power of 5GW. 
e Maximum aperture E-field of 10.6 KV/cm. 


e Maximum array area of 2.5m x 2.0m. 


1. Array Model Simplification 

The array area is given as 2.5m x 2.0m. In order not to exceed the maximum 
aperture E-field constraint, one would expect to have the whole array area filled with 
array elements. Since we are operating at 10 GHz, and the element spacing is 
restricted to less than one wavelength in Z, direction and half a wavelength in X, 
direction, the minimum number of elements needed will be N = 84 x 133 = 11,172 as 
shown in Fig. 4.3(a). With such a large number of elements, it 1s impossible to execute 
the GID2P5 program. Hence, a simplified array model is needed. 

The GTD2P5 code can model the TE,, mode or uniformly distributed 
elements regardless of the element length and width. If we assume that tapering along 
the X, direction is Uniform, and there is no gap between elements in X, direction, then 
we may use a simplified array that has the length of each element unchanged but the 
width extended to include all elements in DG direction. That is, each element’s width 1s 
equal to 2m to cover the whole width of the array area as shown in Fig. 4.3(b). Such a 
simplified array model will produce exactly the same solution as one using the actual 
model. 

Alternatively, if we wish to vary the tapering type along 4, direction, we may 


replace the elements along the Z, direction with a single element of equivalent length 
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(a) Actual Array 


NZ=84 
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(c} Simplified Array Model 2 


Figure 4.3 Array Model Simplification. 


on 


while the width remains unchanged. This configuration is shown in Fig. 4.3(c). 
Because the GID2P5 code can only simulate either Uniform or Cosine (TE,,) 
distribution within the element, we are restricted to using either Uniform or Cosine 
tapering along Z, direction. Together with the use of either a uniform distribution or a 


TE, mode excitation for each element, we have the following four configurations. 


(a) Uniformly tapered TE,, linear array simulated by a single uniformly 
distributed element. 

(b) Cosine tapered TE,, linear array simulated by a single TE,, element. 

(c) Cosine tapered uniformly distributed linear array simulated by a single TE), 
element. 

(d) Uniformly tapered, uniformly distributed linear array simulated by a single 


uniformly distributed element. 


The aperture E-field distribution of these configurations are shown in Fig. 4.4. 

Except for the last configuration, the simplified model will not yield an exact 
solution. Fortunately, the near and far field radiation patterns for the simplified 
models are very similar to those of the actual arravs. The relationship of the gain and 
maximum aperture E-field between the actual and simplified models are derived in 
Appendix C. These derivations allow us to accurately predict the performance of the 
actual array from the data obtained from the simplified model. 

To illustrate the similarity between the actual array and the simplified model, 
5 diagrams consisting of far field and near field pattern plots are shown in Fig. BIS to 
B19. These plots are measured at 10 Km away from the center of the reflector. The 
distance of 10 Km ensures that the field generated by any element within the array area 
of 2.5m x 2.0m will arrive almost in phase, so that direct addition of the individual field 
patterns is possible. The input power is 25W, the mounting plate is 4.5m x 3m and the 
reflector plate is titled at 45°. The elements are arranged in a straight line with a 
constant width of 0.254 along Z, axis. The other parameters of each plot are listed 


below: 


e Model 1 : Actual array for configuration (a) (Fig. B15). 
= Tapering tvpe = Uniform 
= Aperture field distribution = TE,, 
i No. of elements = 20 


= Element size = 0.54 x 0.254 


Actual Array Simplified Model 


Eo Ea 
E, 
HA 
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Ze Zt 
(a) Configuration (a). 
ES Ea 
E. 
Ea 
NZ 
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(b) Configuration (b). 
Eo Fa a 
Anes 
Es 
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(c) Configuration (c). 
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E,. E, 
24 PZ + 


(d)} Configuration (d). 


Figure 4.4 Aperture E Field (E,) for Both the Actual Array 
and the Simplified Model. 


Element spacing = 0.54 
¢ Model 2: Simplified model for configuration (a) (Fig. B16). 
S Aperture field distribution = Uniform 
* No. of elements = 1 
c Element size = 10.04 x 0.254 
e Model 3: Actual array for configuration (b) (Fig. B17). 
7 Tapering tvpe = Cosine 
* Aperture field distribution = TE,, 
t No. of elements = 20 
‘ Element size = 0.5 x 0.252 
, Element spacing = 0.5i 
¢ Model 4: Simplified model for configuration (b) & (c) (Fig. B18). 
¥ Aperture field distribution = TE), 
2 No. of elements = 1 
s Element size = 10.54 x 0.254 
¢ Model 5: Actual array for configuration (c) (Fig. B19). 
cs Tapering type = Cosine 
¥ Aperture field distribution = Uniform 
: No. of elements = 20 
r, Element size = 0.54 x 0.254 


: Element spacing = 0.54 


Fig. BIS and B16 allow the comparison of radiation patterns between the 
actual array and the simplified model of configuration (a). Fig. B17, B18 and B19 
show the radiation patterns of the configurations (b) and (c). Except for the absolute 
magnitude, both the near and far field radiation patterns of the simplified models 
appear to be very similar to those of the actual array. Table 2 lists the maximum 
aperture E-field, boresight far field strength and boresight near field strength of each 
model. It also shows computed ratio of the values of the actual array to those of the 
simplified model. The computed ratios of the maximum aperture E-field and the 
boresight far field strength turn out to be very close to the derived value. Even the 
ratios for the near field field strength also appear to be closely related to the derived 
values for far field. Hence, we can use the data obtained from the simplified model to 


predict the actual array performance accurately. 


TABLE 2 


MAXIMUM APERTURE E-FIELD AND BORESIGHT FIELD STRENGTH COMPARISON 


Max. Oe Boresight Far Boresight Near 
E Fiel Field Field Field Field 
— Strength 


| Model. 1 4095.0 | 4095.0 V/m | 0.0192 0.0192 V/m 23.30 V/m 


2047.0 V/m | 0.0181 V/n 18.30 V/n 


Model 1 
Ratio=—————— 2.0005 eer! 5 ees 3 
Model 2 


Derived Result 
faoG Con- 2.0000 72/32 a 
, figuration (a) 


Model 3 5635.0 V/m 0.0176 V/m 21.40 V/m 


Model 4 3996.0 V/m 0.0196 V/m 23.80 V/m 


Model 3 
Ratio=-—————_- 1.4102 0.8980 0.899 
Model 4 
poneved Result 
for Con- 1.4142 0.8986 
Bieration (b) 





Model 5 2818.0 | 2818.0 V/m 0.0139 | 0.0139 V/m 16.80 | 16.80 V/m 


Model 5 
Ratio=——————_— 0)’; 7051 0.7092 0.706 
Model 4 
Derive result 
Buon Con- 0.7071 0.7058 
figuration (c) 





2. Uniform Field Realization 

The system specifications are derived based on the simplified model of the 
array. It assumes that the field distribution within each element is uniform (not TE, ,) 
and there is no gap between adjacent elements. The area is then occupied bv 96 
elements each of length 0.86874 and width 2m with the tapering along the Z, axis being 
Triangular (linear). This simplified model generates 10.47 KV/cm of maximum 
aperture E-field. The far field patterns generated by this array are attached in Fig. B20, 
and the peak E-field of this array is about 78.81 dB V/m or equivalent to 20.2 W/cm" 
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of power density. The near field patterns are attached in Fig. B21. Fig. B21l(a) shows 
the vertical cut radiation pattern measured at 3.5m. It shows very high radiation level 
in both forward lobes and backlobes. Fortunately, the azimuthal beamwidth is very 
narrow as shown in Fig. B21{b) and the hazardous region is confined to within + 10° 
of the boresight, and is much less for the backlobe since the peaks are much lower. The 
array will be mounted on a vehicle 2m high, hence, it is not possible for the mainlobe 
or the backlobe to point directly at any operator. Fig. B21(c) and B21(d) shows the 
horizontal radiation plots of the array at two sidelobe peaks near the antenna 
mounting plate where the operators could be located. At the front (@ = 90°) and back 
(@ = 270°) of the array, the radiation is still above the hazardous level and hence 


operators should be kept away from the front and back of the array. 


3. TE10 mode Realization 

Although a uniform distribution is easier to simulate and it also simplifies the 
computations, it is not realizable in a slot antenna. Instead, a perfect TE,,) mode is 
chosen to make the simulation more realistic, bearing in mind that the peak E-field in 
each element of the TE,,) mode excitation is much higher than that generated by the 
uniform distribution. 

In order to minimize the aperture field strength, we will fill the 2.5m x 2.0m 
array area with elements close to each other. Using the knowledge we gained in 
Section A of this chapter, we select the element size of 0.8687 x 0.5A so as to reduce 
the aperture field strength while at the same time keeping the sidelobes as low as 
possible. Once the element size is chosen, the number of elements is then fixed at 
96x 133. With this configuration, we can then use the two suggested simplified model 
methods to determine the type of tapering that will meet the requirements. 

After several trials with TE, mode excitations, we find out that by using 
5 GW of input power, the far field criterion can be met by most tapering tvpes. But at 
the same time, the aperture field strength constraint is exceeded. Hence, the input 
power must be lowered. It turns out that when the input power is lowered, only 3 
tapering types will satisfy both the far field criterion and near field constraint. These 
are Uniform - Uniform tapering, Cosine - Uniform tapering and Uniform - Cosine 
tapering. (With the tapering type along Z, direction positioned before the hyphen and 
the tapering type along X, direction placed after it.) 
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Table 3 lists the maximum aperture field strength and the far field boresight 
field strength computed at 10 Km for the 3 tapering types mentioned above. The far 
field field strength is also converted to power density for ease of comparison. Two sets 
of data are obtained for each tapering type. One with 5 GW of input power and the 
other with 1.88 GW of input power to meet the requirements. As can be seen from the 
table, Uniform - Uniform tapering offers the highest far field power density and lowest 
maximum aperture field strength. The Cosine - Uniform and Uniform - Cosine 


taperings provide approximately the same performance. 








TABLE 3 


COMPARISON OF MAXIMUM APERTURE FIELD STRENGTH AND 
FAR FIELD POWER DENSITY FOR DIFFERENT TAPERING TYPE 

























. Max. Aperture Boer te la etree ne ce 
Tapering Field Strength | Boresight Boresight — 
‘ Types (V/m) Field Strength | Power Density , 
at 10 Km (Vm) at 10 Km (W/cm*) 








Fig. B22, B23 and B24 show the near field radiation patterns of the array 
using Uniform - Uniform tapering, Cosine - Uniform tapering and Uniform - Cosine 
tapering respectively. These plots are measured at 3.5m away from the center of the 
reflector plate with 1.88 GW of input power. Note that the NEC-BSC code suggested 
the output is only accurate for the first 30 dB of the pattern, due to the absence of the 
higher order diffraction/reflection terms. For these 3 plots, however, we are only using 
2 flat plates, the antenna mounting plate and the reflector plate, with the reflector plate 


slanting at an angle of 45° from the antenna mounting plate. Furthermore, most of 


Sil 


the energy is concentrated towards the center of the plates and then reflected 
horizontally out of the plates, causing very little interaction within the structure. 
Therefore, we can safely say that the accuracy range limit of 30 dB does not apply to 
these plots. Unfortunately, it is not within the scope of this project to construct an 
actual model for direct measurement and comparison. We will bear in mind the 
accuracy problem and place additional safety margins on the radiation hazard level. 

For Uniform - Uniform tapering, the radiation level is very high at both the 
forward sidelobes and backlobes as can be seen in Fig. B22(a). There may be a 
radiation hazard at 8 around 140° and 230°, assuming that the array is placed on the 
vehicle 2m above ground. The radiation level at the sides of the array (9 = 0° & 180°) 
is also quite high as shown in Fig. B22(b), but it drops down to below 70 dB W/m? for 
@ between 140° and 230°. Fig. B22(c) shows the horizontal radiation pattern at 
§@ = 90°. It shows that the main lobe beamwidth is very narrow and is less than 1° at 
the -10 dB point. Similarly, for the backlobe at @ = 270°, its beamwidth is as narrow 
as the mainlobe. Hence, the radiation hazard regions at 8 around 140° and 230° can 
be expected to be very narrow. For Cosine - Uniform tapering, the problem remains at 
the front and backlobes as in the case of Uniform - Uniform tapering. Its radiation 
level at the sides of the array is very low due to the Cosine tapering along Z, direction. 
For Uniform - Cosine tapering, because Cosine tapering is used along X, direction, the 
vertical cut radiation pattern at @ = 90° shows very low sidelobe and backlobes. 
Hence, there is no possibility of causing a radiation hazard outside of the mainlobe. 
The radiation level along the sides of the array is quite high due to Uniform tapering 
along the Z direction, but it decreases faster than that of the Uniform - Uniform 
tapering and is below 60 dB W/m? for © between 120° and 240°, therefore not causing 
any radiation hazard. 

From the above comparisons, we see that with Uniform tapering along the 4, 
direction, there is very likely to be a radiation hazard due to high sidelobe and 
backlobe levels. Uniform tapering along the Z, direction also causes high radiation 
lobes at the sides of the array. But because the array area is 0.5m (or 16.7 wavelength) 
longer along the Z, direction than along the X, direction, its sidelobe level decreases 
faster and hence is less likely to cause a radiation hazard. Therefore Uniform - Cosine 
tapering type is preferred. 

Fig. B25 and B26 show the far field radiation patterns of the array using 


Uniform - Cosine tapering with the reflector plate tilted at 45° and 22.5° respectively. 
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Although the patterns are not quite as one would expect, the maximum E-field of both 
plots are equal to 0.724E4 V/m. Hence, the tilting of the reflector plate does not affect 


the gain of the array and the far field criterion 1s still met. 


D. VEHICLE MODELING AND FAR FIELD REQUIREMENT 

Fig. 1.1 shows the drawing of an armored vehicle carrying the Flyswatter Array. 
We need to model this vehicle to determine the effects of the presence of the vehicle on 
the radiation pattern of the array and then determine the radiation hazard region cause 
by the entire system. Fig. 4.5 and 4.6 show the computer model input to the NEC- 
BSC code, drawn by GITDGEOMP program. Note that the wheels of the vehicle are 
not shown on the model because they do not affect the radiation pattern significantly. 
The supporting poles and pivots for the array reflector plate are formed by cylinders 
while the rest of the structure is formed by plates of different shapes. Fig. 4.5 is not 
the true plot from GIDGEOMP. The hidden structures have been removed manually 
for easy comparison between the computer generated model and the actual vehicle as 
shown in Fig. 1.1. Fig. 4.6 on the other hand, shows the actual GTDGEOMP plot 
with the array rotated to face the front of the vehicle. The darkened region is the array 
area filled by small elements. ; 

‘The full listing of the input data for Fig. 4.6 with far field pattern plot commands 
is attached in Appendix D, and the far field radiation plots are shown in Fig. B27. 
Fig. B28 shows the far field radiation plots of the same system but with the reflector 
plate tilted at 22.5°. Comparing them with Fig. B25 and B26, suggests that the 
mainlobe and the near in sidelobes have not been affected by the presence of the 
vehicular structure, and the gain also remain unchanged. By rotating the array at 
various azimuth angles and changing the tilt angle of the reflector plate within the full 
gain range of 22.5° to 45°, the gain of the model remain unaltered. Hence, the model 


has met all the requirements specified. 


E. NEAR FIELD RADIATION HAZARD 

Before we evaluate the near field patterns, we need to consider again the accuracy 
of these plots. Unlike previous discussions, we now have many plates and cylinders 
being placed around the array. The vehicular structure plates are placed mostly below 
the array mounting plate. With most of the energy radiating in a narrow beam out of 
the structure, we can be quite sure that there is not much reflected and diffracted field 


from these plates. The cylinders on the other hand are placed above the antenna 


a 





Figure 4.5 Computer Generated System Model with Hidden Structures Removed. 


mounting plate and are subjected to direct radiation. From the previous plot 
Fig. B24(b), we see that there is quite a considerable amount of energy radiating along 
the sides of the array. These radiated energy will be dispersed by the cylinders. Hence, 
we can expect more higher order diffraction and reflection terms to appear and the 
accuracy of the plot is reduced. Using the NEC-BSC code guide line of 30 dB accuracy 
range, We will need to apply more safety margin to the radiation hazard evaluation. 

Fig. B29 and B30 show the near field radiation patterns of the full model with its 
array reflector plate tilted at 45° and 22.5° respectively. The array radiates towards the 
front of the vehicle. From these plots, it is obvious that most of the energy 1s 
concentrated near the main beam. The only region that the radiation is beyond 
hazardous level is within +5° in the elevation angle and +20° in the azimuth angle 
near boresight. With the center of the array reflector erected at about 5m above the 


ground, the safety margin for operator standing at ground level is more than 27 GB. 
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Figure 4.6 Computer Generated System Model. 


Fig. B31 and B32 show the near field radiation patterns of the full model with the array 
rotated to point at @ = 45° and 90° respectively. The reflector plate is tilted at 45° for 
, these plots. Again, there is no sign of radiation hazard near the ground level. Hence, 
for a Uniform - Cosine tapered array, there is a safety margin of 30 dB below the 


allowable radiation hazard level. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


From the preceding investigation, we see that although Uniform tapering 
provides high gain and narrow mainbeam beamwidth, it generates very high sidelobe 
and backlobe levels and is likely to cause a radiation hazard. Cosine tapering on the 
other hand, reduces the sidelobe and backlobe levels by a considerable amount at the 
expense of lower gain and broader mainbeam beamwidth. In order to meet the far 
field requirement of producing 10-20 W/cm? of power density at 10 Km distance and at 
the same time be subjected to the constraint of 10.6 KV/cm of maximum aperture E- 
field strength, we are then restricted to only 3 tapering types, namely, Uniform - 
Uniform, Cosine - Uniform and Uniform - Cosine. After comparing the near field 
radiation patterns, we see that Uniform - Uniform tapering produces the highest 
sidelobe and backlobe levels and is very likely to cause a radiation hazard to the 
Operators working nearby. Cosine - Uniform tapering reduces the radiation level 
towards the sides of the array. But the radiation near the front and back of the array 
is still quite high and 1s also likely to cause a radiation hazard. A Uniform - Cosine 
tapered array operating at 1.88 GW of input power, is capable of producing 
13.91 W/cm? of far field power density at 10 Km range, with a maximum aperture E- 
field strength of 10.6 KV/cm. The near field patterns show that the radiation level near 
ground level is very much lower than with the other two tapering types. With the 
vehicular structure included, we see that there is a safety margin of more than 27 dB. 
below the recommended radiation hazard limit near ground level. Although the NEC- 
BSC code suggested that the accuracy range for the radiation pattern is only up to 
30 dB from the peak value, the contributions from the higher order diffracted and 
reflected terms cannot bring the near field level up by 27 dB. Hence, we can safely say 
that there is no radiation hazard to the operator operating nearby when the Uniform - 
Cosine tapering type is used. 

From this evaluation, we see that the NEC-BSC code is very useful for 
electromagnetic radiation analysis. However, its capability in dealing with complex 
structures is very restricted by its limitation in handling only up to second order 
reflected and diffracted terms excluding doubly diffracted field. It will be able to yield 


more accurate results if higher order terms can be included. The supporting software 
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programs, due to the limited time frame, are written mainly for this project. Further 
improvement will make them fully compatible with the NEC-BSC code. These 


modifications are listed below : 


8 Modify NEC-BSC code to output phase and antenna coupling values. 

9 Modify GIDRWSEL and GIDPLOT codes to handle phase and antenna 
coupling output. 

° Improve the cylinder plotting algorithm of the GTDGEOMP code. 

* Improve the GFDGEOMP code to remove hidden objects. 

° Include a phase variation subroutine in the GIDSA code. 
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APPENDIX A 
PROGRANI FLOW CHARTS AND LISTINGS 


INTERFACING CODE (GIDRWSEL) 
a. Flow Chart 





selects the 
type of olot 
RAMA 


yes CTC 
1YPARMeS ? 
(no aor e 

Alot 


NO 


'§Q. yf 
mes contro 
cCeta fille for 

GTDPLOT 


Ss tA 66 pt 

woltee output 

Gate fille for 
GIOPLOT 
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b. Program Listing 


Cc KAKKAKKKAKKAKAKKRAK KARA 
¢ * GTDRWSEL * 
C * BY KOH WEE JIN * 
s * ON MAY 26 1986 * 
C * AT NPS * 
C KAKAKKAKRAARKEAK RAK RK AK 
C+++ 
Git This Boonen reads in the abbreviated codes and data created 
oe ay GTD2P 
Sate tt then s@lects the plots needed, convert the abbreviated codes 
Ctt++ into useful information and put them into a Control file. 
Crt++ It also selects the appropriate data set and put them into a 
Sait. Data file _ 
aa These two file will be used by GTDPLOT for plotting 
+++ 
REAL DATAIN(1000,5) 
INTEGER IPLOT(4) 
CHARACTER*1 COND(4) 
CHARACTER*3 CONTRL 
CHARACTER*6 YHEAD(4) 
10 CONTINUE 
BOes0 I=1,5 
erat 
C+++ Initialize data set 
C+++ 
DO 20 J=1,1000 
DATAIN(J,1I)=0. 
20 CONTINUE 
30 CONTINUE 
DO 40 I=1,4 
IPLOT(I)=0 
40 CONTINUE 
C++ 
C+++ Read in control and comment data and 
oat Determine if job end is encountered 
+++ 
CALL READC (COND,IXPARM,NPTS) 
CONTRL = COND(1)//COND(2)//COND(3) 
. PecCONIRL .EO. “END') GOTO 1000 
+++ 
eet Convert abbreviated code to appropriate variables 
+++ 
7 CALL CONVRT (COND, IXPARM,ICOORD, YHEAD) 
+++ 
Cft+ Read in plotting data 
ent 
CALL READD (DATAIN,IXPARM,NPTS) 
50 CONTINUE 
ttt 
ew Select the type of plot needed 
+++ 
CALL CSTART eae ches orotate. 
CALE OBLECT (YHEAD, EYPARM, IPLOT, COND 
Pie EveARMasnO- 5) GOTO 10 
Gr++ 
eet Write control data to control file 
+++ 
CALL WRITEC (NPTS,COND) 
ert+ 
Bort Write output data to gtdpdata file 
+++ 


GOTO 
1000 WRITE 
STOP 


CALL WRITED (DATAIN,COND,NPTS, IYPARM) 
(37 '°(a3)') CONTRL 


END 
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BE 
10 CONTINUE 


KRREKKRKKRKRKRKRRRRKRKRRKRRKRREKKRER 


SUBROUTINE CONVRT * 


KREKKKKRKRKRKRKERKRKRKRKKKRKRRRRKKKK 


This subroutine converts the abbreviated codes 
to appropriate variables 


SUBROUTINE CONVRT (COND, IXPARM, ICOORD, YHEAD) 
CHARACTER*1 COND(4) 

CHARACTER*3 YDIM, XAXIS, YAXIS 

CHARACTER*5 COORD (3 2)/'X','¥','Z' SR! ' THETA! 'PHI!/ 
CHARACTER*6 XQ, ¥0 Q, XPARM, YPARM, YHEAD(4 

COMMON /DPARMS/XO,YQ,XPARM, YPARM, YDIM, XAXIS , YAXIS 

IF (COND(3) .EQ. 'R') ICOORD = 1 

IF (COND 3 .EO. 'S') ICOORD 
IF COND {2 2), .EO. 'S') THEN 


XPARM = 'FREQ! 
IXPARM = 4 


END 
IF (Con (2) EQ. ‘F!) THEN ; 
XPARM = 'DIS 
IF (COND (3) 750. 'S' .AND. IXPARM .EQ. 2) XPARM = 'ANGLE! 
IF (COND(3) .EO. 'S' .AND. IXPARM .EO. 3) XPARM = 'ANGLE' 
XO = COORD(IXPARM, ICOORD) 
ENDIF 
DO 10 I = 1,4 
F (1 EQ. . YHEAD(T) = coun 4) BERR K! 
N YHEAD(I) = COND(4)//COORD(I, ICOORD) 





RETURN 
END 


KRRKRKKKKKKRKRKRKRKKRKRKKRKEKKEKKER 


SUBROUTINE CSTART * 


KRKEKKKKKKKRKRKRKRKRKRKRKRRRRKRRKEKK 


This subroutine displays the comment statements 
as a start for the selections that follows 


SUBROUTINE CSTART (COND , ICOORD , IPLOT,, YHEAD) 

INTEGER IPLOT(4) 

CHARACTER*1 COND(4) 

CHARACTER*5 COORD(3,2)/'S' ,'Y' 2 ee 
CHARACTER*6 YHEAD(4 

CHARACTER*60 LABEL1,LABEL2 ,LABEL3,LABEL4, CNTRL 

COMMON /START/FSTART, Raoe TYS, PZS 

COMMON /LABELS/LABEL1, LABELZ, ,LABEL3, LABEL4 , CNTRL 

CALL FRTCMS('CLRSCRN ') 


PRINT = 

PRINT *, 'LABEL1 = ',LABEL1 

PRINT * 

ee eon EQ. a PRINT *, 'FAR FIELD PATTERN PLOT' 

Fe ae it) EO. NOSE HID * 'NEAR FIELD PATTERN PLOT' 

PRUNE 2 \ jae, Thera CONDITIONS OF PLOTS au 

PRINT *, '! ER EQ. = ', FSTART, ',COORD(1, TCOORD), prea 

PRINT oe COORD (2. ICOORD) ,' = "]iya- “COORD (3, ICOORD), 

= ',P2S 

PRINT * 

DO 10 I = 


1,4 
IF (IPLOT(I) .EQ. 1) PRINT *, YHRAD(2) 3" IS BEING FROGEsser: 


10 CONTINUE 
PRINT * 


RETURN 
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- KRKEAEKKRRKKRKRKRRKRRKKRKKK 
c * SUBROUTINE READC * 
C KRRERKKRKKRRRREKRRRKEKKRERE 
C++ 


C+++ This subroutine reads in the comments and control cards 


ot 
SUBROUTINE READC (COND,IXPARM,NPTS) 
CHARACTER*1 COND(4) 
CHARACTER*2 CT 
CHARACTER*60 LABEL1,LABEL2,LABEL3,LABEL4,CNTRL 
COMMON /START/FSTARI ,RXS,TYS,PZS 
COMMON /LABELS/LABEL1 , LABEL2 , LABEL3 , LABEL4 ,CNTRL 
READ (1,100,END=1000) CT,LABEL1 
100 FORMAT (1X,A2,1X,A60) 
ieee «NE. 'CT') LABELI='S' 
READ (1,110) COND,IXPARM,NPTS,FSTART,RXS,TYS,PZS 
110 eee el ziaie as), 2b2-2) 


1000 CONTINUE ° 
e+ 
at Indicate the end of the input data set 
sae 
COND(1) = 'E! 
COND(2) = 'N'! 
COND(3) = 'D' 
20 CONTINUE 
RETURN 
END 
c 
Cc KAKKKKARKKKKKRKKRKRKKRKRKKK 
G x SUBROUTINE READD x 
C KAEKKKKKKKKKKKKRKKKKRKKKKKK 
C+t++ 


C+++ This subroutine reads in the data 


SUBROUTINE READD (DATAIN,IXPARM,NPTS) 
REAL X(4) ,DATAIN(1000,5) 
ICHANG = Q 
Be lO I = 1 NPIS 
READ (1,100) X(4),X(1),X(2),X(3), (DATAIN(I,J) ,J=2,5) 
100 MORMAT (lx, FP? .3,3F9.3,4E10.3 
IF ee sem 1) Xi2)o= 360.8 =5x(2) 
IF (IXPARM .EQ. 2 .AND. X(2) .EQ. 180.) ICHANG = 1 
DATAIN (I,1) = X(IXPARM) 
10 CONTINUE 
URN 


Ret 
END 
s2uaeanaeweae eee eee eee eee eeaeee eee eee eeeeeeeee eee eeeeee eee eee eet ee ee eee e == 
Cc KERAAKAKKKKKKKKAKKKKKKKKKK 
e * SUBROUTINE SELECT * 
Cc KAARKAKKKKKAKKKKAKKKKKKAKKE 
Cr+ 


C+++ This subroutine displays the various options 
C+++ and reads in the choices 


SUBROUTINE SELECT (YHEAD,IYPARM, IPLOT, COND) 

INTEGER IPLOT(4) 

CHARACTER*1 COND(4 

CHARACTER*3 YDIM,XAXIS, YAXIS 

CHARACTER*6 X (10, SPARM , YPARM, YHEAD (4) 

CHARACTER*8 YUNITS(2,3)/'V/M ' 'DBV/M =! , 'A/M 
'DBA/M  ','W/SOM_ ','DBW/SQ M'/ 

COMMON /DPARMS/XQ,YQ,XPARM, YPARM, YDIM, XAXIS , YAXIS 

10 CONTINUE 
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Ctt+ 


C+++ Select the type of plots 


Ct++ 


TS 


PRINT *, 'SELECT ONE OF THE FOLLOWINGS : | 


PRINT. ~,. de Plot ee vHEAD GE 
PRINT *, '2) PLOT ', YHEAD(2 
PRINT *, '3) PLO vhbeaos 
PRINT *, '4) PLOT ', YHEAD(4 
PRINT *, '5) NOP 


RI "ENTER THE NUMBER OF THE PLOT.’ 
READ (5,'(1I1)') IYPARM 
LF aaa -~LT. 1 .OR. LYPARM 2GT. 5) GOunemre 
COND(1) .EQ. 'F' .AND. IYPARM .EQ. 1) THEN 
CALL FRTCMS('CLRSCRN ') 
PRINT *, 'THERE IS NO !',YHEAD(1),'! FOR FAR FIELD PLOT' 
PRINT *, 'RE-ENTER PLOT SELECTION' 


ENDIF 

IF (IYPARM .EQ. 5) GOTO 1000 

Ory BEAD eae 

IF (IPLOT(IYPARM) .EQ. 1) THEN 2 
CALL FRTCMS('CLRSCRN ') 
0 


LS: b= 1G 
IF (IPLOT(I) .EQ. 1) PRINT *, YHEAD(I),' HAS BEEN PROCESSED' 
CONTINUE : 

PRINT *, 'RE-ENTER PLOT SELECTION! 

PRINT * 

GOTO 10 


END IF 
IPLOT(IYPARM) = 1 


20 CONTINUE 


C+++ 
C+++ 
CH++ 
C+++ 


C+++ 
C+++ 


Select magnitude or phase plot . 
(Presently, this code only caters for magnitude plot) 


CALL FRTCMS('CLRSCRN ') 

PRINT * 

PRINT * 

PRINT *, 'SELECT ONE OF THE FOLLOWINGS : ' 


PRINT *, oy MAGNITUDE PLOT' 

PRINT *, '2) PHASE PLOT! 

PRINT *, 'ENTER THE NUMBER OF THIS PLOT.' 
ISELEC = 0 

READ (5,'(1I1 LSELEC 


mt) 
IF (ISELEC .LT. 1 JORa ISELEG) Gtnez) pcowomee 
LF ( ISELEC EO. 5 YPARM = ‘MAG! 

ER (LSELECG Wr EOee 2) eCeann 'PHASE' 


Select the Y-axis! unit. 


IF (COND(4) .EQ. 'E') IYUNIT = 1 

IF (COND(4) .EO. 'H') IYUNIT = 2 

IF (COND(4) .EO. 'P') IYUNIT = 3 
30 CONTINUE 

CALL FRTCMS('CLRSCRN ') 

PRINT * 

PRINT * 


C++ 


PRINT *, ‘SELECT THE UNITS TO BE USEDE 
PRINT *, 2 ee enon 

PRINT *, '2) ' ,YUNDTTSt2nDyUNtT 

PRINT *, 'ENTER THE NUMBER OF YOUR CHOICE.' 
PoELEC = 0 

READ (5, (01)*) 1Seube 

IF (ISELEC .LT. 1 .OR. ISELEC .GReeZ economy 
TES Cis leew bo. 3 YDIM = 'LIN' 

IF (ISELEC..EO. 2) e7p Mig DE. 


40 CONTINUE 
- 
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C+++ Select linear or log scale plot 


C+++ 
CALL FRTCMS('CLRSCRN ') 
PRINT * 
PRINT * 
Print oso SELECT THES@yreno OF PLOT:: ' 
PRINT *, '1) X : LINEAR, Y : LINEAR' 
PRINT *, '2) X ;: LINEAR, Y : LOG 
PRINT *, '3) X : LOG, Y : LINEAR' 
PRINT *, '4) X : LOG, Y : LOG 
PRINT *, 'ENTER THE NUMBER OF YOUR CHOICE.' 
ISELEC = 0 
READ (5,'(I1 


J SP rORLEE 

iM tSELEG .LT. 1 .OR. ISELEC .GT. 4) GOTO 40 
XAXIS = 'LIN' 

YAXIS = ‘'LIN' 

Pie lsnLee .EQO. 2 .OR. ISELEC .EQ. i YARIS 
TeeonLee .EO. 3 .OR. ISELEC .EO. 4) YAXIS 


1000 CONTINUE 
RETURN 


ET 
END | 

C 

c RAKKKKKAKKKAKRKKKKKRAKRAKRKRRK 

¢ x SUBROUTINE. WRITEC  * 

Cc KRARAEKKRKRKKKKKRKRKRRKKKKKKRARKK 


C+tt+ This subroutine prepares and writes the various labels 
C+++ and fields needed to create the control file 


SUBROUTINE WRITEC (NPTS,COND) 
CHARACTER*1 ND(4 
CHARACTER*3 XAXIS,YAXIS,YDIM 
CHARACTER*4 UNITSS 
CHARACTER*5 NPTS 

CHARACTER*6 RXCS,TYCS,PZCS,XQ,YQ,XPARM, YPARM, FIELDS 
CHARACTER*7 S 

CHARACTER*9 RXS,TYS,PZS 

CHARACTER*60 LABEL1,LABEL2,LABEL3, LABEL4, CNTRL 
COMMON /START/FSTART,RXS,TYS,PZS | 

COMMON /DPARMS/XQ,YQO,XPARM,YPARM, YDIM, XAXIS,YAXIS 
COMMON /LABELS/LABEL1,LABEL2,LABEL3, LABEL4, CNTRL 


OQ 
O 


rrj 
— 


C+++ 
C+++ Converts numerical variables to characters 


Seatt 
OPEN (4,STATUS='NEW' ) 
WRITE (4,100) NPTS,FSTART,RXS,TYS,PZS 
meOehORMAT {1X,14,1X,F7.3,3F9.3) 
CLOSE (4) 


OPEN (4,STATUS='OLD' ) 
Renee 4,110) NPIS,FSS,RXS,TYS,PZS 
110 FORMAT(1X,A5,A7,3A9) 
CLOSE (4,STATUS='DELETE' ) 
eet 
eet Sreate the control statement 
+44 
CNIRL = 'NEW '//NPTS//X0//' '//XPARM//' '//¥Q//' '//Y¥PARM//' '// 
: AYDIM//' '//XAXIS//' YARIS 
+++ 
eae Frepare the labels 
ere+ aver displays the title of the plot - 
eet Label 2 to 4 display the initial conditions 
at 
ENGtoS =" MM! 
RXCS = '(X)= ! 
ues = 'CY)= | 
PZCS = '(Z)= ! 
ie GGlibas ))-bO. 1R*) GOTO 10 
UNITSS = ' DEG' 
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RKC$ = '()R = ! 
TYcs = 38 all 
PZCS = !'? = aul 
10 CONTINUE 
LE ee EQ. i FIELDS = ' (F)AR' 
IF (COND(1) .EQ. 'N') FIELDS = ' (N)EAR' 
LABEL2=FIELDS//' FIELD (GTD) PLOT WITH Serene Sie seu FOLLOWSS ' 
LABEL3=' FRE L//ESS/ 7) GE 2 ne S | 7 ee 
nee LABEL4=' '//TYC$//TY$//UNITSS//' WB ACS/ IPD /Honttss//' S! 
+++ 
C+++ Writes the control statement and Labels ontosthe controle. 
Cia 
WRITE (3,'(A60)') CNTRL 
WRITE (3,'(A60)') LABEL1 
WRITE (3,'(A60)') LABEL2 
WRITE (3,'(A60)') LABEL3 
WRITE (3,'(A60)') LABEL4 
RETURN 
END 
C 
C KAAKKARKAKKARKARAAKAKKAKA 
c * SUBROUTINE WRITED * 
C KRAKKKRKKKKARKARAKKRKRKAKR 
Cars : 


C+++ This subroutine writes the appropriate data to GIDPDATA file 


SUBROUTINE WRITED (DATATN COND ,NPTS , LYPARM) 
REAL DATAIN(1000,5) 

CHARACTER*1 COND(4) 

CHARACTER*3 YDIM, XAXIS,YAXIS 

CHARACTER*6 KQ YO /¥Q, XPARM, YPARM. 

COMMON  /D /DPARMS/XO, 9, XPARM , YPARM, YDIM, XAXIS, YAXIS 


ic TYPARM +1 
DO 10 I = 1,NPTS 
(YD 














LE IM (COS Je tHen 
IF (DATAIN(I,IY) .LE. 1.E-20) DATAIN(I,IY)=1.E-20 
IF (COND(4) .EQ. 'E') DATA=20.*ALOG1O(DATAIN( L, seg 
IF (COND(4) .EQ. 'H') DATA=20.*ALOG1O(DATAIN(I,IY 

a IF (COND(4) .EQ. 'P') DATA=10.*ALOG1O(DATAIN(I,IY 
DATA=DATAIN(I,IY) 

ENDIF 


WRITE (2, 00 be 1), DATA, DUMMY 
100 FORMAT (5X,3 3) 
10 CONTINUE 
RETURN 
END 


30 


OUTPUT DATA PLOTTING CODE (GTDPLOT) 
a. Flow Chart 


Mafn Program 





Flow Chart For GTOPLOT (cont!nue) 
Subroutine PARMS 


start 


reecs {n the 
contro] 

{nformation 

& ecnos !t 






CTLe°ENO’, eats 
END=. TRUE. 


return 





extrects @xe3 

levels, ecales, 
dimensions & 

types of olot 
usfndq XTRACT 


reece In 
taoe ls & 
ecros tres 


return 


Ca 
tJ 


AAARAANAAN 


aAAAN 
++++ 
++++ 
t++++ 


ett 
C+++ 


C+++ 
C+++ 
C+t++ 


C+++ 
ert + 
er+t 


at 
cra + 
C+++ 


1010 
C+++. 


C+++ 
C+++ 


b. Program Listing 


KRAAKKRKKRKKKRKKRKKKRKKKRKKKKKKRKKEK 


* GTDPLOT * 
* MODIFIED FROM NFPLOT * 
* BY KOH WEE JIN * 
x ON MAY 29 1986 * 
* * 
* * 


AT NPS 
KAKKKKKKKKKRKARKARKRRKRKERKKKK 


Bo rogram reads in the Control and Data files created by 
GTDRWSEL. 
Tt then select the SPE LCR tnoeS labels and scales and plot types 
and plot them using Disspla subroutines 

LOGICAL*4 STRTED 

LOGICAL*4 END /.FALSE./ 

COMMON /DISPLA/ STRTED 
Turn off the start indicator 

STRTED = .FALSE. 
Initialize Disspla 

CALL DISSET 
Plot until ‘END! 


CONTINUE 


control statement is detected 


Read control statements and labels to be plotted. 


CALL PARMS(END) 
IF (END) GO TO 1020 


Plot the graph 
CALL PLOT 
Indicate that we have completed at least one plot. 


STRTED = .TRUE. 
GO TO 1010 


End of plot. Otherwise, 
routine. 


CONTINUE 
CALL DISEND 
SLOP 


program abort through an error handling 


KRAKKKKKKKKKKRRKRKKKKKKKKE 


* SUBROUTINE ALGPLT * 
KAKAKAKKAAKKKAAKKKAKKRKKK 


as er outine replaces the Disspla subroutine ALGPLT. 
Bae oUt des a rounded axis origin and log cycle length in 
ine es per cycle. 


Modification is provaded to eliminate truncated axes syipeomatic 
Mielodicithmie plots. Only full log cycles are returne 


SUBROUTINE ALGPLT(AMIN, AMAX, AXLEN, AO, ACY) 
REAL*4 AMIN 


5 


Cate 
C+++ 
C+++ 
C+++ 


1020 


REAL*4 AMAX 
REAL*4 AXLEN 
REAL*4 AO 
REAL*4 ACY 
REAL*8 DMIN 
REAL*8 DMAX 
REAL*8 DXLEN 
REAL*8 DAO 
REAL*8 OMAGO 
REAL*8 OMAG1 
INTEGER*4 MAGDIF 
Do it all in double precision 
DMIN = AMIN 
DMAXK = AMAX 


DXLEN = AXLEN 
Get the order of magnitude of the minimum value 
OMAGO = DLOG10(DMIN) 


If the minimum 1s an even power of 10, then use it as the lower 
limit. Otherwise, go to next lower order of magnitude. 


IF (IDINT(OMAGO) .NE. OMAGO) OMAGO = IDINT(OMAGO - 1.) 
Do the same for the maximum. 
OMAG1 = DLOG10(DMAX) 


If the minimum is an even power of 10, then use it as the upper 
limit. Otherwise, go to next higher order of magnitude. 


IF (IDINT(OMAG1) .NE. .OMAG1) OMAG1 = IDINT(OMAG1 + 1.) 
Compute the origin 


CONTINUE 
DAO = 10.D0**IDINT(OMAGO) 
IF (DAO .GT. DMIN) THEN 
-OMAGO = IDINT(OMAGO) - 1. 
GO TO 1010 
ENDIF 


Compute the difference in order of magnitude and the number of 
user units (usually in inches) per cycle of log. 


CONTINUE 

MAGDIF = (OMAG1 - OMAGO) - 1. 

IF ((DAO 10.DO**MAGDIF) .LT. DMAX) THEN 
OMAG1 = OMAG1 + 1. 
GO TO 1020 

ENDIF 

AO = DAO 

ACY = DXLEN / MAGDIF 

RETURN 

END 


RAEKAARKAKKARARAARARARARRARARARAARKRKRE 


* SUBROUTINE BLSRCH * 
KAAKKKKAKARKAAKAAKKAAKKK 


This subroutine scans for the first non-blank character and 
returns with the pointer (I) set to point to the character. 


SUBROUTINE Se ENa| cae I) 
CHARACTER*1 LINE(80) 
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E CHARACTER*1 BLANK /' '/ 
+++ 
C+++ If we have gone beyond the end of the line, complain and abort. 
a 
IX =I 
IF (I .GE. 81) GO TO 8010 
Crt 
C+++ Scan for a non-blank 
er 
DO 1010 I = IX, 80 
IF (LINE(I) .NE. BLANK) RETURN 
1010 CONTINUE 
oot 
Gia Error message. 
or 
8010 CONTINUE 
CALL ERRMSG(IX, 'Unexpected end of line detected.$') 
CALL DISEND 
STOP 400 
END 
C eeeeeeeeeew eee eeeee eee eee ee ee eee ee ee eee eee wee eee eee eee ee eee eee ee ee eee ee = 
5 RAKKAKKAARKKAKAKKEKKAKAKK 
C * SUBROUTINE DISEND * 
Cc KAKKAKKARKKKKKKKKKRKRKKKK 
at 
ear End of Disspla 
ot 
: SUBROUTINE DISEND 
LOGICAL*4 STRTED 
COMMON /DISPLA/ STRTED 
ek tne CALL DONEPL — 
IF (.NOT. STRTED) WRITE(6,9010) 
RETURN 
9010 gob (ONO PLOTS GENERATED... PROGRAM ABORT. ') 
Cw aw wo ww wee ewww www ewww eee eee ce ee ee we ww eew ewe n ween wneaaasonoae 
: RARKAAKAKKAKAKKAKRKAKAAKK ° 
C * SUBROUTINE DISSET * 
Cc KAKKARKKKAAKAKKKRKKKAKK 
Cassa 
oa This subroutine establishes the parameters for Disspla. 
+++ 
SUBROUTINE DISSET 
CALL COMPRS 
See CALL SETCPR(O, 1, 0, 0) 
CALL PAGE(8.5, 11.0) 
art 
C+++ Lettering is FUTURA with upper case and lower case. 
ore, Upper case letters are delimited by ‘(' escape character. 
+++ 
CALL BLOWUP (.5) 
CALL FUTURA 
CALL BASALF('L/CSTANDARD' ) 
CALL MIXALF('STANDARD' ) 
CALL MX3ALF('INSTRUCTION', '|') 
CALL MX4ALF('L/CGREEK', '?') 
CALL NOBRDR 
C+++ 
et Integer (or rounded) axes with Y axis labelling at 0 degrees. 
+++ 
CALL INTAXS 
CALL YAXANG(O. ) 
RETURN 
END 
ell | 
Cc RAKKKAKAKKKKKRAKAKRAKKKRK 


> 


aul * SUBROUTINE ERRMSG * 
C KKAKKK 


KKK KAKKKKRKRRRRKA 
Cid 
Sak This subroutine writes an error message to the logical unit 6. 
+++ 
SUBROUTINE BRRUSEAN, MSG) 
CHARACTER*1 DOLLAR /'$'/ 
CHARACTER*1 MSG i 
CHARACTER*1 BAR /'_'/ 
CHARACTER*1 BLANK /' 
Ct+++ 


coe Indicate the error by writing the bar below it. 
+44 
WRITE(6,9010) (BLANK, I = 2, N), BAR 
C+t++ 
C+++ Find the end of the message. 
Crt 
DO 1010 [=1 


ee 
IF (MSG(T) -BO, sDOMEAR NGO TO 1020 
1010 CONTINUE 


I = 81 
crt 
C+++ Write the message. 
C+++ 


1020 CONTINUE , 
WRITE(6,9010) (MSG(J), J = 1, I) 
RETURN 

9010 FORMAT (1X,80A1) 


END 

G s 

C KRAKKAKKRAAKRKAR RAR KARR IK 

e * SUBROUTINE EXPONE * 

Cc RREKRRRKRRRKRKKAKRKRKAAKKKERE 

Gra 

oe This subroutine concatenates the exponent (if required). 

det 
SUBROUTINE EXPONE(INVEC,* EXPON) 
CHARACTER*1 INVEC(64) 
INTEGER*4 EXPON 
INTEGER*4 , FILL /ZOOOO0O00FO/ 
CHARACTER*1 HEXCHR(4) 
INTEGER*4 DUMMY 
SET CoE (HEXCHR(1) , DUMMY) 

(EXPON .EQ. 0) RETURN 
C+t++ 


C+++ Find the end of the input vector and generate the exponent. 


Gates 
DO ee 1 = 
F CINVEC(L) .EQ. '$') GO TO 1020 
1010 aetna 
CALL east 'Unable to concatenate. Vector full.$') 
1020 “ENVEC(),= | ; 


tee 2 
CALL TIMS10(INVEC(L), I) 
IXP = EXPON 
TR CIxe ure ft THEN 
IxP = a IXP) 
INV -! 
i °F t) 
ENDIF 
Oe = Oe 


Pe 
IF (IXX .NE. 0) THEN 
DUMMY = IXX + FILL 
INVEC(I) = HEXCHR(4) 
IXP = IXP - (IXK * 10) 
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ENDIF 
DUMMY = IXP + FILL 
pUeK 2) = HEXCHR(4) 





IL 
INVEC(I) = ‘8! 
RETURN 
END 
c 
Cc KAKKAKKKKKKAKAKKARKAKKKA RAK 
C * SUBROUTINE GETTOK * 
Cc KAKKKKKKKKEKKKKKKKKKKKA 
ert+ 
aut This subroutine extracts an 8 character token. 
+++ 
SUBROUTINE Se eee I, TOKEN) 
CHARACTER*1 LINE(80) 
CHARACTER*8 TOKEN 
CHARACTER*8 TOK 
CHARACTER*1 oe 
EQUIVALENCE (CHAR(1), TOK) 
CHARACTER*8 BLANKS /' My 
CHARACTER*1 BLANK /' '/ 
CHARACTER*1 COMMA /','/ 
C+++ 
Ct+t++ Clear the token. 
C+++ 
TOK = BLANKS 
Gr++ 
C+++ Find the first non-blank character. 
C+++ 
CALL BLSRCH(LINE, I) 
att 
Suet Extract the token (up to 8 characters). 
C+t++ | 
N= 1 
Ix =I 
DO 1010 I = IX, 80 
IF acta SEO ea GO TO 1030 
IF (LINE Eee COMMA) GO TO 1020 
IF (N .EQ. 9) GO To 8010 
Crr+ 
et Character Ok. So put it into the token buffer. 
$+< 
See, & LINE(1) 
1010 CONTINUE 
I = 81 
RETURN 
ert+ 
eat Exit if comma found after bumping the pointer. 
bo 
1020 CONTINUE 
iat + 1 
1030 CONTINUE 
TOKEN = TOK 
RETURN 
ert++ 
C+++ Error, token length is > 8. 
eat 


8010 CONTINUE 
CALL ERRMSG(I, 
‘Command word or token length exceeds 8 characters.$') 
CALL DISEND 
STOP 500 
END 


C 
Cc KEKEKKKKKKKKKKKKKKKAKRKE 


oF 


{I | 
c 


C+++ 
C+++ 
C+++ 


C+t++ 
C+++ 
C+++ 


C+++ 
C+++ 
C+++ 


A2LNANA 


* SUBROUTINE LINLIN * 
KRAKKKKKKKKKKKKKKKKKKKK 


This subroutine sets up DISSPLA for a Linear-Linear axis plot. 
SUBROUTINE LINLIN 


REAL*4 AMIN 
REAL*4 XMAX 
REAL*4 YMIN 
REAL*4 YMAX 


COMMON /RANGE/ XMIN, XMAX, YMIN, YMAX 
CALL GRAF(XMIN, 'SCALE', XMAX, YMIN, 
RETURN 

END 


'SCALE', YMAX) 


KERERKRKKKRKKARKRRERRRKKRKKK 


* SUBROUTINE LINLOG * 
KAAKKAKKKAKKAKKAKRKAKKA 


This subroutine sets up DISSPLA for a Linear-Log axis plot. 
SUBROUTINE LINLOG 
REAL*4 XMIN 


REAL*4 AMAX 
REAL*4 YMIN 
REAL*4 YMAX 


COMMON /RANGE/ XMIN, XMAX, YMIN, YMAX 
Get the proper scaling for both axes. 


CALL UCoL Cae XMAX, 6.0, XORIG, XSTEP, DUMMY) 
CALL ALGPLT(YMIN, YMAX, 6.0, YORIG, YCYCLE) 


Set up DISSPLA for a Log scale on the Y axis. 
CALL YLOG(XORIG, XSTEP, YORIG, YCYCLE) 


RETURN 
END 


REKKKRKKKRRKAKKRRRKKERKRK 


* SUBROUTINE LGLG * 
KEKKKAKKKAKKAKKAKKKERK 


This subroutine sets up DISSPLA for a Log-Log axis plot. 
SUBROUTINE LGLG 


REAL*4 AMIN 
REAL*4 XMAX 
REAL*4 YMIN 
REAL*4 YMAX 


COMMON /RANGE/ XMIN, XMAX, YMIN, YMAX 
Get the proper scaling for both axes. 


CALL ALTE aN, XMAX, 6.0, XORIG, ea 
CALL ALGPLT(zMmN, YMAX, 6.0, YORZG © YeverE 


Set up DISSPLA for a Log-Log scale plot. 
CALL LOGLOG(XORIG, XCYCLE, YORIG, YCYCLE) 


RETURN 
END 


RAAKRRRKKRRKARKARRRRKRKRKRKKE 


* SUBROUTINE LOGLIN * 
KRAKRKKKRAKKAKRAKKRAKARR 
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C+++ 
C+++ 
C+++ 


C+t++ 
C+++ 
C+t+ 


C+++ 
C+++ 
C+++ 


er. 
C+++ 
oats 


C+++ 
C+++ 
C+++ 


C+++ 
Crrt 
C+++ 


1010 


C+++ 
C+t+ 
C+t+ 


8010 


This subroutine sets up DISSPLA for a Log-Linear axis plot. 
SUBROUTINE LOGLIN 


REAL*4 AMIN 
REAL*4 XMAX 
REAL*4 Se 
REAL*4 


COMMON /RANGE/~ RMIN, XMAX, YMIN, YMAX 
Get the proper scaling for both axes. 


CALL Se a aE XMAX, 6.0, XORIG, XCYCLE) 
CALL AXSPLT(YMIN, YMAX, 6.0, YORIG, YSTEP, DUMMY) 


Set up DISSPLA for a Log scale on the X axis. 
CALL XLOG(XORIG, XCYCLE, YORIG, YSTEP) 
RETURN 


KRAKAKKKAKKRRKKKRKKKARRKRKAE 


* SUBROUTINE NUMBER * 
KRAKKKAKKKAKKAKKAKRAKRKE 


This subroutine extracts an integer from the data stream. 


SUBROUTINE Se I, N) 
CHARACTER*1 LINE a 
CHARACTER*1 PREG 
CHARACTER*1 ZERO /'0'/ 
CHARACTER*1 NINE on) 
CHARACTER*1 BLANK /! '/ 
INTEGER*4 ZIP 
CHARACTER*1 ocean /3*Z00,'0'/ 
neces (NONE(1), ZIP) 
INTEGER*4 NUM 
CHARACTER*1 BYTES (4) /4*Z00/ 
EQUIVALENCE (NUM, BYTES) 
Extract the token. 
CALL GETTOK(LINE, I, TOKEN) 
Get the number. 
DO 1010 J = 1, 8 
IF ens) .EQ. oo RETURN 
(TOKEN(J) .LT. ZERO) .OR. (TOKEN(J) .GT. NINE)) GO TO 8010 
Extract the digit and get it into the number. 
eek a ee 
+ (NUM - ZIP) 
a 
RETURN 


EMror. Not a digit. 


CONTINUE 

CALL ERRMSG(I, 
CALL DISEND 
STOP 


'Not a digit.$') 


KRAAKKKRKKKKKKKKKKKKKKKE 


oP 


alll * SUBROUTINE PARMS * 

Cc KAKKKKKAKKKKKKKAKKKKAKKA 

Crt 

C+++ This subroutine reads control parameters and data 

C+++ The control statement syntax lis: 

C++4 

C+++ ine lL) == 

Grit 

CHt+ <STATUS><NPTS><X Q><X PARM><Y¥ Q><¥ PARM><Y UNIT><X AXIS><Y¥ AXIS> 
et+ 


C+t++ where. <STATUS> is "OLDUNEWio) OF sunious 
C+++ <NPTS> is the number of points to be plotted (omit if 
C+++ 'OLD' of 'END' is selected), 
Crt <X Q> Loa oe a) | Ra ‘THETA! , Ht’ or "PREG 
C+++ <X PARMS rs Dist ea MAG, 'ANGLE' or 'FREQ! 
ea+¢+ <¥ Q> Sig.en Ey. 'EPHI', "Eo, EloeweA' ia 
ea + 'EPEAK', 'HR', 'HX', 'HPHI', 'HY', 'HTHETA", 
er++ 'HZ', 'HPEAK', 'PR', 'PX', 'PPHI', 'PY'!, 
C+++ 'PTHETA', "P42! or ‘ERR 
Cae <Y PARM> is 'DIST!, 'MAG', or 'PHASE' (<X PARM> and 
C+++ <¥ PARM> are mutually exclusive), 
C+++ <Y¥ UNII> 1s "LIN or YDB! 
C+++ <X AXIS> is 'LIN' or 'LOG', and 
Girt: <Y Seiiseis “LIN! or ‘hoc’. 
C+++ 
Ct+++ line 2 to 5 are <LABEL>'s 
Catt 
C+t++ where <LABEL> is a label up to 60 characters (letters to 
C+t++ be capitalized must be enclosed in parentheses) 
oT terminated by the '$' character. 
s+ 

SUBROUTINE Paes 

REAL*4 ATA(1000,3) 


INTEGER*4 NPTS 

INTEGER*4 Pl 

INTEGER*4 P2 

CHARACTER*1 LABEL1(80 

CHARACTER*1 LABEL2(80 

CHARACTER*1 LABEL3(80 

CHARACTER*1 LABEL4(80 

INTEGER*4 XTYPE 

INTEGER*4 YTYPE 

INTEGER*4 XQUAL 

INTEGER*4 YOUAL 

INTEGER*4 Y 

COMMON /DPARMS/ DATA, NPTS, Pl, P2, LABEL1, LABEL2, LABEL3, 
LABEL4, XTYPE, YTYPE, XQUAL, YQUAL, YU 

LOGICAL*4 


END 
CHARACTER*1 CTL(80) 
CHARACTER*8 SELECT(5) is te cee 


PHASE! , 
eae 


— XQUALS (7) je REQ"/ 
WANE 

'R! 
| 


CHARACTER*8 YUNIT(2) /\LIN’, 
CHARACTER*8 YQUALS(21) 1 EB agi 'ETHETA', 'EY', 'EPHI', 'EZ', 

'HR', 'HX!', 'HTHETA!, 'HY', 'HPHI', 'HZ', 
'PR! ASTER! 'PTHETA', 'PY!', 'PPHI', 'PZ!, 
'PPEAK! / 


60 


C+++ 
C+++ 
oo 


C+++ 
C+++ 
Crtt+ 


C+t++ 
C+++ 
C+++ 
C+++ 
oa 
ot 
C+++ 


C+++ 
C+++ 
C+++ 


1010 


C+++ 
C+++ 
C+++ 


C+++ 
C+t+ 
C+t+ 


1020 


C+++ 
C+++ 
ont 


1030 


C+++ 
C+++ 
C+++ 


C+++ 
C+t+ 
C+++ 


C+++ 
C+++ 
C+++ 


CHARACTER*8 COMAND(3) /'NEW', 
OLD’, 


'EN 

CHARACTER*8 TYPES(2) eet 
LOGICAL*4 DATAIN /.FALSE./ 

Read the control information and echo it. 


READ(5,9010, eee) Cre 
WRITE (6, 9020) C 


Parse the control information 
I=1 


Data source: 


OLD -- use data already read in. . 
NEW -- read a new block of <NPTS> points 
END -- end of job 


eX TRAGIKCIL, 1, 


CAL COMAND , 3, JUMP) 
Gouto (1010, UMP 


1030, Fi0Z0)r 
'NEW! paaected: 
CONTINUE 
INUMER = I 


CALL NUMBER(CTL, I, NPTS) 

only 1000 or less points allowed. 
IF (NPTS .GT. 1000) GO TO 8020 
CALL READIN 

GO TO 1030 

'END' selected. 

CONTINUE 
END =. 
RETURN 


TRUE. 


Process the X and Y PARM specifications 


CONTINUE 
GAM ATRACT(CIL, I, x UALS, 7, XQUAL) 
GhemeXTRACTV GEL, I, SELECT, 5! PI) 
Gham XTRACT( ETL, I, HQUALS, 2 21, a 
CALL XTRACT(CTL, I, SELECT 2 
PRINT *, eSuaES YOURE} kguALs{xouaL 
PRINT *, 'YOQUALS(YOUAL)=' , YOUALS(YOUAL 

ny “Sa 

'YSELECT(P2)=! 


PRINT x! 


a SEL Bee 
PRONT *, 


“SELECT P2 
Select the Y axis unit. 


CALL XTRACT(CTL,I,YUNIT,2,YU) 
PRINT *, 'YUNIT(YU)=', YUNIT(YU) 


Select the axis type 
CALL oaefen Ga i) TYPES ac eat yre 
Chil TRAGCL (elm, L, TYPES, Z, YIYPE 
PRINT *, Bet es Sie ce carpe, 
PRINT * | Pie bow CL PE j—a hye ks (YTYPE 


Read in the labers 
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C+++ 

C+++ 

C+++ 
8010 


8020 


READ(5,9010,END= B03 0) LABEL1 
WRITE (6, 9020) LABEL1 
READ(5,9010,END=8030) LABEL2 
WRITE (6 9020) LABEL2 
READ(5,9010,END=8030) LABEL3 
WRITE (6 9020) LABEL3 

READ ( ,9010, END=8030) LABEL4 
WRITE (6, 9020) LABEL4 

RETURN 


Error messages 


CONTINUE 

CALL ERRMSG(1, 
CALL DISEND 
SLOP 100 
CONTINUE 

CALL ERRMSG(INUMER, 
CALL DISEND 
SLOnaaoe 
CONTINUE 

CALL ABEOF('PARMS ', 5) 
See e ee 
FORMAT (1X, 80A1) 

END 


'No END statement found.$' ) 


‘Number of points exceeds 1000.$') 


KAAKKAKKKKKAKAKKKKKKKKAARKK 


*x SUBROUTINE PLOT * 
KKAKKKKAAKAKKAKARKAKKK 


This subroutine plots the data. 
SUBROUTINE PLOT 


REAL*4 DATA(1000,3) 
INTEGER*4 NFTs 
INTEGER*4 Pl 
INTEGER*4 EZ 
CHARACTER*1 LABEL1(80 
CHARACTER*1 LABEL2(80 
CHARACTER*1 LABEL3(80 
CHARACTER*1 LABEL4(80 
INTEGER*4 ALLE E 
INTEGER*4 Leer 
INTEGER*4 XQUAL 
INTEGER*4 YOQUAL 
INTEGER*4 YU 
COMMON /DPARMS/ DATA, NPTS, Pl, P2, LABEL1, LABEL2, LABEL3, 
LA BEL4, XTYPE, ‘YTYPE, XQUAL, YQUAL, YU 
INTEGER*4 EXPONX 
INTEGER*4 EXPONY 
REAL*4 AMIN 
REAL*4 AMAX 
REAL*4 Sl 
REAL*4 


YMAX 
COMMON /RANGE/ XMIN, XMAX, YMIN, YMAX 
CHARACTER*64 NAMVEC 


CHARACTER*12 NAME(5) /'(D 


CHARACTER*16 XUNITS(7) / 


ISTANCE |, 
AGNITUDE ', 
HASE ', 
NGLE ! 
REQUENCY 


= 


- -mpry 
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VUUU 


, EXPONX, XMIN, XMAX) 


NPTS 
63 


) 


| / /YUNITS (YQUAL, YU) 


| / /RUNITS (XQUAL) 
'NAMVEC=! , NAMVEC 


NAMVEC 
Lei, 
, 100) 


a 
a 


| NAMVEC 


NAMVEC, EXPONX 


LABEL4, 


C+++ Label the axes 


C+++ 


LABEL1, 
LABEL2, 
LABEL3, 
DATA(1 

(P2)// 


NAMVEC 





| 


NAY 


'NAMVEC='! 
NAME 


NAME(P1)//' 


a 


CALL HEADIN 
CALL HEADIN 
CALL HEADIN 
CALL HEADIN 
CALL SCLDAT 
CALL EXPONE 
ERONT *, 
CALL ANAME 


PRING * 
NAMVEC 
PRINT *, 


C+++ Write the headers 
NAMVEC 


eas + 
So 
C+++ 


CALL ERPONE (NAMVEC. ,1P2), NPTS, EXPONY, YMIN, YMAX) 
CALL EXPONE(NAMVEC EXPONY ) 


PRINT *, 'YNAME= »NAMVEC 
CALL YNAME (NAMVEC, 100) 
é PRINT *, 'XTYPE=',XTYPE, ' YTYPE=" ,YTYPE 
tre 
Ct+++ Dispatch to the correct routine. 
C+++ 
IF XTYPE .EQ. 1) . AND. (7a vEre cee CALL LINLIN 
rs XTYPE .EQ. 1) . AND. (oye onc CALL LINLOG 
IF MIYPE .EQ. 2) -AND. (7iveE eo. CALL LOGLIN 
IF XTYPE «EQ. 2) SANDR (vit Gee ec CALL LGLG 
C+++ 


C+++ Draw the curve 


C+t+ 
CALL CURVE(DATA(1,IP1), DATA(1,IP2), NPTS, 0) 
CALL RSCDAT(DATA(1, IPL), NPTS, EXPONX, XMIN, RMAX} 
CALL RSCDAT(DATA(1,IP2), NPTS, EXPONY, YMIN, YMAX 
CALL ENDPL(0) 
RETURN 
END 
G 
C KAEKKKKKKKKKKKKKEKKKREKK 
E * SUBROUTINE READIN * 
C KAKKKEKKKRKKERKKKKKEKKKKKKK 
C+++ 
C+t++ This subroutine reads in the data. 
C+t+ 
SUBROUTINE READIN 
REAL*4 DATA(1000,3) 
INTEGER*4 NPTS 
INTEGER*4 Pl 
INTEGER*4 P2 
CHARACTER*1 LABEL1(80 
CHARACTER*1 LABEL2(80 
CHARACTER*1 LABEL3(80 
CHARACTER*1 LABEL4(80 
INTEGER*4 XTYPE 
INTEGER*4 YTYPE 
INTEGER*4 RQUAL 
INTEGER*4 
COMMON /DPARMS DATA. NPTS, Pl, P2, LABEL1, LABEL2, LABEL3, 
- | LABEL4, XTYPE, ‘YTYPE , XQUAL, YQUAL 
+++ 
C+t+t+ Read the data 
C+++ 
DO) LOLS i =a 
READ(1,9010, END= 73010) (DATA(I,J), =e) Se) 
LORG CONTINUE 
RETURN 
C+++ 
C+t+t+ Abends 
C+++ 


8010 CALL ABEOF('READIN', 1) 
9010 FORMAT(5X,3F10.3) 


END 
C 
Cc KRAAKKKAKKKKAKKKAKKKKAKKKA 
C * SUBROUTINE RSCDAT * 
C KREAKAKKKEKKKKKRKKKKKKKKKK 


C+++ This subroutine re-scales the data 


SUBROUTINE SS OH arene NPTS, EXPON, MIN, MAX) 
REAL*4 DATA(NPTS 

INTEGER*4 NELS 

INTEGER*4 EXPON 
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C+++ 
C+++ Generate scale factor and compute the data 


C+++ 
SF = 10**EXPON 
Domo t= 1, NPT 
DATA(I) = ‘DATA(T) * SF 
1010 CONTINUE 
MIN = MIN * SF 
MAX = MAX * SF 


RETURN 

END 
C 
C KAKKKAAAKAAAKKKAKKKKKRR 
C * SUBROUTINE SCLDAT * 
C KKKKKAKKRARKKRAAARKAKKKK 


C+++ This subroutine scales the data (if necessary) 


SUBROUTINE Sat er NPTS, EXPON, MIN, MAX) 
REAL*4 DATA(NPTS} 
INTEGER*4 NPTS r 
INTEGER*4 EXPON 
REAL*4 MIN 
REAL*4 MAX 
C+++ 
Boe Get order of magnitude of minimum and maximum 
$e 
- IF (MIN .NE. 0) T 
IOMMIN = 6610 (ABS (MIN) ) 
E 
IOMMIN = 0 
ENDIF 
IF (MAX .NE. 0) THEN 
IOMMAX = LOG10 (ABS (MAX) ) 
EDSE 
IOMMAX = 0 
ENDIF 
C+++ 


C+++ If the order of magnitude is within 10**(+/-3), then just set 
C+++ exponent to 0 and don't scale 


C+++ 
IF Se GE. i ~AND. aa ae a ~AND. 
IOMMAX .GE. -3) .AND. (IOMMAX .LE. 3))) THEN 
O — 


ELSE 
C+++ 
C+++ Compute a suitable scaling factor 


C+++ 
ee Boson) Si 3) THEN 
(TOMM MME) 3) * 3 
SF is = tol wk IOMMAX 
EXPON = IOMMAX 
ENDIF 


IF (ABS ( LOMMIN) rok. 2). T 
IOMMIN = (LOM MIN / 3) Ne 3 
SF = 10.**IOMMIN 
EXPON = IOMMIN 
ENDIF 


DO 1010 I = NPTS 
DATA (I) = DATA(I) / SF 
1010 CONTINUE 

MIN = MIN / SF 

MAX = MAX / SF 
ENDIF 
RETURN 
END 


c 











Cit 
Obese 
Ctr+ 


Crt 
Carat 
ort 


KRRAERKRRKRKRKKKKRKRRKRKRKRKKEKR 


* SUBROUTINE TIMS10 * 
KAAAKRKAKKRARKKRRAKKRRKK 


This subroutine puts the '*10**' string in 


SUBROUTINE TIMS10(INVEC, I) 
CHARACTER*11 INVEC 
INTEGER*4 I 

INVEC = 'X 10]EH.75)' 

ie = tel 

RETURN 

END 


RERRRRRKRKRRARRRRRRRKRKRRKRK 


* SUBROUTINE XTRACT * 
KAAKKAKKAAKKAKKRAKKRAKKK 


This subroutine scans for a token, checks it against the token 
table (word) and returns the index of the token in the table 


SUBROUTINE XTRACT(LINE, I, WORD, NWORDS, INDEX) 
CHARACTER*1 LINE(80 

CHARACTER*8 WORD(NWORDS ) 

CHARACTER*8 TOKEN 


Get the token 


DO 10 J=1,NWORDS 
PRINT *, WORD(J) 


10 ae 


Cana 
Cath 
Ca++ 
C+++ 


ERR = I 
CALL GETTOK(LINE, I, TOKEN) 


Scan the table for the token. If 1t Is fet Gere, write aneenccus 
message and terminate. 


DO 1010 INDEX = 1, NWORDS 
IF (WORD( INDEX) .EQ. TOKEN) RETURN 


1010 CONTINUE 


Crt. 
Cre 
ert 


Error message 


CALL ERRMSG(IERR, 
CALL DISEND 

SEOP 300 

END 


'Unidentifiable command or token.$') 
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3. GEOMETRIES PLOTTING CODE (GTDGEOMP) 
a. Flow Chart 


Mafn Program 


MCX 3; max. No. of cylinders 
NPX ; @ax. No. of pletes 


NSAX : max. No. of source éerrays 





NSX : max. No. Of sources 


eets Ww Slocx 
Cats 





reecs In 


e'fae requenty fn 


/ 


no 
yes L@=, TRUE. 
C2=°"@* 2 
é) 
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Fiow Chart For GTDGEOMP (continue) 
Mafn Program 





plete dets 











reecs ate 
Osts & Ost. 
o locatic 





forms exes 
trenmsformation 
matrix 










Flow Chart For GTDGEOMP (continue) 
Mafn Program 
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Flow Chart For GTDGEOMP (continue) 
Subroutine PLOT 





the cy] {noer 
4,° id a 
a BH 
Check Ciots qond 
UE PLU’ 
= plots pietes 


| Sus. ENOL O) | 
7 — 


~) 
> 


C+++ 
C+++ 
C+++ 
att 
C+++ 
C+++ 
Ct++ 
eat 
oot 
C+++ 
C+++ 
et 
C+++ 
C+++ 
C+++ 
at 
ort 
C+++ 
C+++ 
C+++ 
et 
oa 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
ert 


C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
ert 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+t+ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 


b. Program Listing 


KKEKKKKKKKEKKEKRERERKEKRERKRKERERKRKKA 
GTDGEOMP 


WRITTEN BY KOH WEE JIN 
Chips SEP 1986 


AT NPS 
KAAKAAKAKKKKARKKAKKAKKKKKKKKKKKKK 


THIS PROGRAM READS IN THE DATA FROM GTD 
IT THEN PLOTS THE FOLLOWING GEOMETRIES 
CYLINDER(S 
INFINITE GROUND PLANE 
PLATE(S) 
SOURCE(S) 
SOURCE ARRAY(S) 
THE PLOTS INCLUDES THE EFFECTS DUE TO 
AXES ROTATION 
ORIGIN TRANSLATION 
SCALING 


A A Ot 
HO OO OO 


MAXIMUM 
MAXIMUM 
MAXIMUM 
MAXIMUM 
MAXIMUM 
MAXIMUM 
MAXIMUM NUMBER 


PARAMETER 
PARAMETER 
PARAMETER 
PARAMETER 
PARAMETER 
PARAMETER (NSAX=5) 

PARAMETER (NESAX=500) 


DECLARE VARIABLES 


PLATE 
CORNER LOGALTIONS = XPEC,YPC,2PC 
Nee ORSELATES = IPG 
NO= OF CORNERS PER PLATE : MEP 


CYLINDER 
CUNBER 3s) sCCYL, YCCyYL, ZCCYL 
ENIDRGAP CENTERS, = 2CN,2CP 
END CAP LOCI 
END CAP CUT ANGLES : THIN, THTP 
AXES ROTATION MATRIX : RCYL 
AXES LENGTH FROM CENTER : AC,BC 
NOewOreerLINDER. : ICG 
SINGLE SOURCE 
CORNERVPLOACTIONS = XSC,7YSC,2ZSC 
NO. OF SOURCES : ISG 
SOURCE ARRAY 
CORNER LOACTIONS : XSAC,YSAC,ZSAC 
NO] Ob BoOURGES. : 254 
NO OF ELEMENTS PER SOURCE ARRAY 
TOTAL NO. OF ELEMENTS : MSAA&T 
INFINITE GROUND PLANE 
WITH GROUND PLANE : LGP = .TRUE. 
WITHOUT GROUND PLANE : LGP = .FALSE. 
& MAX. KX, Y & Z VALUES COMPARISON 
CLE ENPER =: -GOMUPCG 
PEATE | :, COMPPG 
SOURCE = COMPSG 
SOURCE ARRAY : 


NUMBER PLATES 

EDGES PER PLATES 

CYLINDERS 

POINTS FOR CYLINDER END CAPS 
SINGLE SOURCES 

SOURCE ARRAYS 


ELEMENTS PER SOURCE ARRAY 
RESO 

NEX=12 

NCX=100) 


NCECP=201) 
NSX=1000) 


MIN. 


COMPSA 
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: MSAX 


INPUT CARDS 


NP& 
NEX 


NGECF 
NSX 
NSAX 
NESAX 


eerur rec r1 Le yoer Le, wcrLN, YCYLN, Z2CYLN 


Crt OBSERVATION POINT IN GLOBAL PLANE 


C+++ ANGLE FROM X-AXIS : PHI 

Ctt+ ANGLE FROM Z-AXIS : THETA 

C+++ RADIUS FROM ORIGIN : RADIUS 

C+++ GLOBAL AXES ROTATION MATRIX : ROT 

C+++ GLOBAL AXES ORIGIN TRANSLATION : XCR,YCR,ZCR 
Cnet FREQUENCY IN GHZ : FRE 

C+++ SCALING FACTOR : SCALE 

C+++ INTEGER DEFINED BY UN CARD : IUN 

crt+ INTEGER DEFINED BY US CARD : IUS 

ae 


REAL COMPPG(3,2),COMPCG(3,2),COMPSG(3,2) ,COMPSA(3,2) 
REAL ROT eo 

REAL XPC(NPX,NEX) , YPC(NPX,NEX) ,ZPC(NPX,NEX) 

REAL REE Nea eae aan Ne SECO SCT Ne aieaee. 
REAL XCYLN(NCX ,NCECP) , YCYLN(NCX,NCECP) , ZCYLN(NCK ,NCECP 
REAL TE CaRTT NG Tot eee caren ea 
REAL XCCAPN(NCX) , YCCAPN(NCX) , ZCCAPN(NCX 
REAL XSC(NSX,4) , YSC(NSXK,4) ,ZSC(NSX,4 
REAL XSAC(NSAX ,NESAX,4) , YSAC(NSAX,NESAX,4) ,ZSAC(NSAX ,NESAX, 4) 
INTEGER MEP(NPX) ,MSAX(NSAX) 

CHARACTER*60 TITLE,LABEL 

CHARACTER*2 C2 

LOGICAL LGP,LGEOM 


C+++ 
on SET UP COMMON BLOCKS 
+++ 
COMMON /CGDATA/ ICG,COMPCG,XCCAPP ,YCCAPP,ZCCAPP, 
* XCCAPN , YCCAPN , ZCCAPN 
COMMON /OBDATA/ XOBS,YOBS,ZOBS,ROBS 
COMMON /OSDATA/ POSMAX,TOSMAX, PHIOS , THETOS 
COMMON /PGDATA/ IPG,MEP,COMPPG 
COMMON /PIDATA/ PI,RTD 
COMMON /RTDATA/ XCR,YCR,ZCR,ROT 
COMMON /SADATA/ ISA,MSAX,COMPSA,MSAXT 
COMMON /SCALER/ FREQ,IUN,IUS,SCALEF 
COMMON /SGDATA/ ISG,COMPSG 
COMMON /VUDATA/ PHI,THETA, RADIUS 
TEILE='5S! 
LGP=.FALSE. 
LGEOM=.FALSE. 
IST=0 
C+++ 
Ctt++ INITIALIZE DISsees 
C+++ 
CALL DISSIN 
C+++ 
Ce READ IN KEY LETTERS AND SELECT THE BRANCHES 
tt 
10 CONTINUE 
READ (1,'(A2,1X,A60)') C2,LABEL 
IF (€Z .EO. “Gh 3GolerZ0 


IF (C2 .EO. 'CG') GOTO 30 
IF (C2 .EO. 'EN') GOTO 40 
IF (C2 .EO. 'FR') GOTO 50 
IF (C2 .EO. 'GP') GOTO 60 
IF (C2 .£0, NG conor 70 
IF (C2 .EO. 'NG') GOTO 80 
IF (C2 .EO. 'NP') GOTO 90 
IF (C2 .EO. 'NS') GOTO 100 
IF (C2 .EO. "NK scone 110 
IF (C2 .EO. 'PG') GOTO 120 
IF (C2 .EO. "RE economies 
IF (C2 .EO. SA!) ecomoui4e 
IF (C2 .EOnSG ir Gomemi5u 
IF (C2 .EO. 'UF') GOTO 160 
IF (C2 .EO. 'UN') GOTO 170 
IF (C2 .EO. 'US') GOTO 180 
IF (C2 .EO. 'XQ') GOTO 190 


a2 


GOTO 10 


C+t+t+ 
C+++ CT : COMMENT TITLE CARD 
C+++ 
20 CONTINUE 
WRITE 3+ 
WRITE (2,%*) 'CT! 
TITLE=LABEL 
WRITE (2,%*) ! Pith ee=- LLLLE 
GOTO 10 
C++t+ 
C+++ CG : CYLINDER GEOMETRY CARD 
C+tt 
30 CONTINUE 
WRITE Bra 
WRITE (2,%*) 'CG' 
LGEOM=.TRUE. 
ICG=ICG+t1l 
WRITE (2,%) ! IcG =!',ICG 


IF (ICG .GT. NCX) GOTO 35 
CALL See een iage ut cer nN»YCYLN , ZCYLN } 
GOTO 1 
35 PRINT *, 'NO. OF CYLINDERS EXCEEDS NCK'! 
CALL DONEPL 
SLOP 


C+++ 
C+++ EN : END CARD 
C+++ 
40 CONTINUE 
WRITE (203 
WRITE (2,*) 'EN'! 
CALL DONEPL 
STOP 
C+++ 
Crom PR =: FREQUENCY CARD 
C+++ 
50 CONTINUE 
WRITE (2.5 
WRITE (2,*) ‘FR! ° 
LGEOM=.TRUE. 
READ (1,%*) FREQ 
WRITE (2,*) ! FREQ =',FREQ 
GOTO 10 
C+++ 
C+++ GP : GROUND PLANE CARD 
C+++ 
60 CONTINUE 
WRITE ra 
WRITE (2,*) 'GP! 
LGEOM=.TRUE. 
LGP=.TRUE. 
WRITE (2,*) ! EGE =" snGe 
READ (1,*) LSLAB 
IF (LSLAB Be 0) GOTO 10 
READ (1,*) DUMMY, DUMMY, DUMMY, DUMMY 
GOTO 10 
C+++ 
C+++ NC : NO CYLINDER CARD 
C+++ 
70 CONTINUE 
WRITE con 
WRITE (2,*) 'NC! 
LGEOM=.TRUE. 
ICG=0 
Gabe NOevL( XCYLP, YCYLP, ZCYLP, XCYLN, YCYLN,ZCYLN) 
GOTO 10 
C+++ 
C+++ NG : NO GROUND PLANE CARD 
C+t+ 
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80 CONTINUE 
WRITE Bo 
WRITE (2,*) 'NG! 
LGEOM=.TRUE. 
LGP=.FALSE. 

GOTO 10 


C+++ NP : NO PLATE CARD 


90 CONTINUE 
WRITE or 
WRITE (27) Ne 
LGEOM=.TRUE. 
IPG=0 
CALL Nee TE DESERT SNE 
CALL RESET2(NPX,MEP,1,DUMMY,1,DUMMY) 
GOTO 10 

C+++ 


C+++ Ns): NO SCURCE CARD 
C+++ 
100 CONTINUE 

WRITE ek 
WREIE@(Z,7) “Ns” 
. LGEOM=.TRUE. 
ISA=0 
ISG=0 
CALL RESED 1 
CALL RESET1 
CAEL RESET2Z 
GOTO 10 

Cr+ 


Staak NX : NEXT SET OF NEW DATA CARD 
b+ 
110 CONTINUE 
WRITE eo 
WRITE (2,*) 'NX! 
LGEOM=.FALSE. 
CALL NXT(XPC,YPC,ZPC , XCYLP, YCYLP ,ZCYLP, XCYEN  YGYLN Zevee 
*XSC,YSC,ZSC, XSAC, YSAC, ZSAC, LGP, TITLE) 
GOTO 10 
C+++ 


C+++ PG :,. PLATE CARD 
C+++ 
120 CONTINUE 
WRITE oe 
WRITE (2,*) 'PG! 
LGEOM=.TRUE. 
IPG=IPG+1 
WRITE (2 © ag! IPG =',IPG 
IF (IPG .GT. NPX) GOTO 125 
CALL PLATEG(XPC,YPC,ZPC) 
GOTO 10 
125 PRINT *, 'NO. OF PLATES EXCEEDS NPX'! 
CALL DONEPL 


1 ,NSX,4,XSC,YSC,ZSC, COMPSG) 
NSAX ,NESAX,4,XSAC ,YSAC,ZSAC,COMPSA) 
NSAX,MSAX,1,DUMMY,1,DUMMY) 





SLOP 
C+++ 
tm RT : ROTATE TRANSLATE AXES CARD 
+++ 
130 CONTINUE 
WRITE Sr 
WRITE (2,*%) 'RT' 
LGEOM=.TRUE. 
See XCR, VER AGCR 
READ(1,*) THZP,PHZP,THXP, PHXP 
CALL ROTATE(THZP,PHZP,THXP,PHXP,ROT) 
GOTO 10 
C+++ 
C+++ SA : SOURCE ARRAY CARD 
C+++ 
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140 CONTINUE 
W 


RITE 13 A 
Whine (2,7) oA! 
LGEOM=.TRUE. 

ISA=ISAt1 
WRITE (2,%*) | oA =',1SA 
IF (ISA .GT. NSAX) GOTO 145 
IST=ISA+ISG+MSAAT 

IF (IST .GT. NSX) GOTO 147 
Sire gene ere eae) 


OTO 
145 CONTINUE 


PRINT *, 'THE TOTAL NO. OF SOURCE ARRAYS = ',ISA 
PRINT x 'NO. OF SOURCE ARRAYS EXCEEDS NSAX'! 
CALL DONEPL 


SLOP 
147 CONTINUE 


C+++ 
C+++ 
C+++ 


PRINT *, 'THE TOTAL NO. OF SOURCES = ',IST 
PRINT *, 'NO. OF SOURCES EXCEEDS NSX' 

CALL DONEPL 

SLOP 


SG : SINGLE SOURCE CARD 


150 CONTINUE 


WRITE $3 ie 

BRE (2,7) 'SG' 
LGEOM=.TRUE. 

ISG=ISG+tl 

WRITE (2,%*) ! ISG =',1ISG 
IST=ISA+tISG+MSAAT 

ieers st §6.GT. NSX) "GOTO 155 
CALL SORCEG(XSC,YSC,ZSC) 
GOTO 10 


155 CONTINUE 2 


C+++ 


PRINT *, 'THE TOTAL NO. OF SOURCES = ',IST 
PRINT *, 'NO. OF SOURCES EXCEEDS NSX'! 

CALL DONEPL 

SOP 


C+++ UF ;: GEOMETRY SCALE FACTOR CARD 


er++ 


160 CONTINUE 


C+t++ 


WRITE tS x 

Weaoe (2,%) ‘UF! 

LGEOM=.TRUE. 

READ > SCAL Ee 

WRIT E (2 *) SCALEF =' ,SCALEF 
GOTO 10 


C+++ UN : PLATES AND CYCLINDERS GEOMETRY SCALING CARD 


C+++ 


170 CONTINUE 


et 


WRITE (2, R 

WRITE (2,%) ‘UN! 

LGEOM=. TRUE . 

READ(1,*) I 

WRITE (2,%) Hh IUN =', IUN 
GOTO 10 


C+t++ US : SOURCES GEOMETRY SCALING CARD 


Cir? 


180 CONTINUE 


WRITE os 
WRITE HuS * 


READ(1,*) IUS 


WRITE (2,%) ! IUS =',IUS 
Gano 10 


is 


Cate 
Ctt+ “AOR S exe CuPrOnN CARD 


C+++ 
190 CONTINUE 
WRITE 30 
WRITE (2,%* A 
IF (.NOT. LGEOM) GOTO 10 
CALL PLOT(XPC,YPC,ZPC, XCYLP , YCYLP , ZCYLP  XGYEN eye vclle2ey ee 
AXSC,YSC,Z5SC,XSAC ,YSAC,ZSAC,LGP, TITLE) 
LGEOM=.FALSE. 
GOTO 10 
END 
C se © eee es 2 eee eee e2e2e2eeaeene2e2aeee2eee2ee2eaeaenaeeneae2enaeaeaeaanenenaeae2eneneneaane2e2eaaeaeaenaeaaaaeaae & 
Gc see eet e2 ene eet een een e2ee2eenet eae een Z2e2eeneee2e2eeee2ae2nteaeaee2eeaeaee2eeeee ®e2eee 22 © © 2 @ @ @ @ @ 
CkkX 
cxxx SUBROUTINE BLANKC 
Crk 
C+++ 


C+++ THIS SUBROUTINE DETERMINES THE PORTION OF ELETPSE OF THE 
C+++ CYLINDER TO BE BLANKED OFF 


SUBROUTINE BLANKC(ICYL, XCYL,YCYL,ZCYL, ICBDY,PNS,LBLNK1,LBLNK2) 

PARAMETER ieute 

PARAMETER (NCECP=201 ) 

REAL ee GORE oa 

REAL XCYL(NCX,NCECP) , YCYL(NCX,NCECP) , ZCYL(NCX ,NCECP) 

REAL X(NCECP) ,Y(NCECP) ,Z(NCECP) 

REAL A CADN INNER UE CHO NCH Coc oneee 

REAL XCCAPN(NCX) , YCCAPN (NCX) ,ZCCAPN (NCX 

CHARACTER PNS 

LOGICAL LBLNK1 , LBLNK2 

COMMON /CGDATA/ ICG,COMPCG,XCCAPP,YCCAPP,ZCCAPP, 
* XCCAPN , YCCAPN ;ZCCAPN 

COMMON /OBDATA/ XOBS,YOBS,ZOBS,ROBS 

LBLNK1=.TRUE. 

LBLNK2=.TRUE. 

WRITE 3+ 

WRITE (2,%) 'SUBROUTINE BLANKC'! 
IF ene FEOG an WRITE 13 POSITIVE END CAP' 
IF (PNS .EO. 'N') WRITE (2,*) ' NEGATIVE END CAP! 
DO 10 I=1,NCECP 
AX(I)=XCYL(ICYL,I 
Y( T)=YCyYLt revi i 
20) =Z2CYLCreyYL,.1 

10 CONTINUE 

C+++ 
C+++ TRANSFER THE LOCUS OF THE END CAP TO OBSERVATION CO-ORDINATE 


C+++ 
CALL ROTRAN(NCECP,X,Y,Z,XOBS, YOBS,ZOBS ,ROBS) 
C+++ 


C+++ TRANSFER THE CENTER POINTS OF THE TWO END CAPS TO OBSERVATION 
C+++ CO-ORDINATE 


C+++ 

XCP=XCCAPP (ICYL 

YCP=YCCAPP( ICYL 

ZCP=ZCCAPP(ICYL 

XCN=XCCAPN( ICYL 

YCN=YCCAPN( ICYL 

ZCN=ZCCAPN( ICYL 

CALL ROTA a SReH ON SON CO CtEcE ae ee 
CALL ROTRAN(1,XCN ,YCN ,ZCN,XOBS, YOBS, ZOBS, ROBS 
+t+ 


C+++ CONSIDER ONLY THE X-Y PLANE 

C+++ DETERMINE THE INTERCEPTION POINTS BETWEEN THE ELLIPSE AND 
C+++ THE CENTER LINE FORMED BY JOINING THE CENTERS OF THE TWO 
C+++ END CAPS 


C+++ 
NHALF=(NCECP-1)/2 
ISLOPE=1 
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IF (ABS(XCP-XCN) .LT. 0.00001) GOTO 15 
SLOPE=(YCP-YCN) / (XCP-XCN) 
DO 20 I=ICBDY,ICBDY+NHALF 
7 
IF (I1 .GE. NCECP) I1=I1-NCECP+1 

= +1 
Water occa eta yanee) 
YI2=YCP+SLOPE*(X(I2)-XCP 
IF (Y a .GE. YI1 .AND. ae mis Ee es GOTO 30 
tie vem ee khee ile nlp. Y(I2) .GE. YI2) GOTO 30 


20 CONTINUE 
GOTO 80 
15 CONTINUE 


IF 
25 CONTINUE 
GOTO 80 


ISLOPE=0 

Fi ee I=ICBDY , ICBDY+NHALF 

IF (I1 .GE. NCECP) I1=I1-NCECP+1 

I2=I1itl 

IF oe .GE. XCP .AND. ae LE. ee GOTO 30 
Mme Lhoener Gay. Avic) «Gk. KEP) GOTO 30 


30 CONTINUE 


oat + 
Gert 
C+++ 
ot 


IBLNK1=I11 
IBLNK2=I1+NHALF 
IF (IBLNK2 .GT. NCECP) IBLNK2=IBLNK2-NCECP 


PetonLNE THE DISTANCES FROM THE INTERCEPTION POINTS TO THE CENTER 
FOUND OF THE OPPOSITE END CAP 


XC=XCP 
YC=YCP 

Ze=7CP 

IF (PNS .EQ. 'N') GOTO 40 
XC=XCN 

YC=YCN 

ZC=ZCN 


40 CONTINUE 


ert 
Cat 
Coat 
Cert 


D1=SORT ( (X(IBLNK1)-XC)*(X(IBLNK1)-XC)+ 
*(Y(IBLNK1)-YC)*(¥(IBLNK1)-YC 

D2=SORT ( (X(IBLNK2) -XC)*(X(IBLNK2)-XC)+ 
* (Y(IBLNK2)-YC)*(¥(IBLNK2)-YC) ) 


COMPARE THE TWO DISTANCES AND INCULDING IN THE DEPTH (Z CO-ORD) 
TO DETERMINE IF THE PORTION IS BEING BLOCKED OR NOT 


Ite (DL ce. D2) GOTO 50 
PeCrsLOre «EO. ©) GOTO 55 
TheCDe.GT. DZ) GOTO 60 

enn cet ns zee eEnIGe 
mee a ee ) .GT. ABS(ZCI)) LBLNK1=.FALSE. 


60 CONTINUE 


BMT ION Sone oo BRE 
ee eens ) .GT. ABS(ZCI)) LBLNK2=.FALSE. 


55 CONTINUE 


Ee. GT. DZ).GOTO 70 
See ee eee tees 
ie eee EEN ) .GT. ABS(ZCI)) LBLNK1=.FALSE. 


70 CONTINUE 


ee esa ae acco eee 
IF (ABS(Z(IBLNK2)) .GT. ABS(ZCI)) LBLNK2=.FALSE. 


50 CONTINUE 


ae Ce.” ) ‘SUBROUTINE BLANKC END' 


RETU 
80 PRINT *, 'THE INTERCEPTION POINT WAS NOT FOUND '! 


PRINT *, 'PROBABLY THE ALGORITHM IS NOT GOOD ENOUGH' 


ue 


PRINT *, ‘TRY OTHER METHOD! 
enae DONEPL 





Te 
END 
ce a G0 am a GO =O GD GD GD G&D & 20 ED OP OD ED ad 2D 2D OD GD GD OD GD OD GD GD GD GD GE GD GD GD GD MD GD GD aD a ad OD OO GD GD GD GD ab GD aD OD GD GD GD GD GD GD GD GD GD GD GD GD ae ae GD aS aD aD me > 
Cx 
c*x*xX BLOCK DATA 
CARR 
BLOCK DATA 
C+++ 
C+++ SET UP PARAMETERS 
C+++ 
PARAMETER (NPX=50) 
PARAMETER NCX=100) 
PARAMETER NSX=1000) 
PARAMETER (NSAX=5) 
C+++ 
C+++ DECLARE VARIABLES AND ARRAYS 
C+++ 
REAL COMPPG(3,2) ,COMPCG(3,2) ,COMPSG(3,2) , COMPSA(3,2) 
REAL ROT(3,3),ROBS(3,3) 
REAL XCCAPP(NCX RC EOWIERY cecatniet, 
REAL XCCAPN(NCX) , YCCAPN(NCX ) ,ZCCAPN(NCX 
P INTEGER MEP(NPX) ,MSAX(NSAX) 
+++ 
ai SET UP COMMON BLOCKS 
+++ 
COMMON /CGDATA/ ICG ,COMPCG ,XCCAPP , YCCAPP,ZCCAPP, 
* “ACCAPN , YCCAPN , ZCCAPN 
COMMON /OBDATA/ XOBS,YOBS,ZOBS,ROBS 
COMMON /OSDATA/ POSMAX,TOSMAX, PHIOS , THETOS 
COMMON /PGDATA/ IPG,MEP,COMPPG 
COMMON /PIDATA/ PLR 
COMMON /RTDATA/ XCR,YCR,ZCR,ROT 
COMMON /SADATA/ ISA,MSAX, COMPSA ,MSAXT 
COMMON /SCALER/ FREO; LUN, LUS ,SCALEE: 
COMMON /SGDATA/ ISG, COMPSG 
- COMMON /VUDATA/ PHI,THETA, RADIUS 
+++ 
C+++ INITIALIZE DATA SETS 
C+t++ 
DATA COMPCG /3*10.0E10,3%*-10.0E10/ 
DATA COMPPG /3*10.0E10,3*-10.0E10/ 
DATA COMPSA /3%*10.0E10,3%*-10.0E10/ 
DATA COMPSG /3%*10.0E10,3%*-10.0E10/ 
DATA XCR,YCR,ZCR /3*0.0/ 
DATA ROT /1.0,3%0.0,1.07320.0,31007 = 
DATA FREQ, IUN,IUS,SCALEF /0.2997925 ,1 Ge i. Oy 
DATA PHI,THETA,RADIUS /30.0,60.0,1.0/ 
DATA ICG,IPG,ISA,ISG,MSAXT /5%*0/ 
DATA PI,RTD /3.14159,57.29578/ 
DATA POSMAX , TOSMAX/4.0,9.0/ 
END 
CxARK 
C*x** SUBROUTINE COMP 
Cx 
C+++ 
C+++ THIS SUBROUTINE COMPARES THE OVERALL MIN. & MAX. X, Y & Z VALUES 
ae IT ALSO COMPUTES THE PLOTTING SIZE NEEDED BY DISSPLA 
+++ 


SUBROUTINE COMP(C,C1,C2,C3,C4,L) 

REAL a 
I 
I 
I 


; 'SUBROUTINE COMP! 

iE 

LE) 2h. Clee CUT). ) SO 

: (GiGi, I C(I,1)=C2 a 
I, 


<Glaeeoun, L CU lj) = 
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Melee Gianc4(1,1)) C(I,1)=C4(I,1) 
10 CONTINUE 


DO 20 I=1,3 
EMG hea el(i.2,).) C(T.2)=C1(I.,.2 
ie CCG eT. COUmeae C(T, 2)=Cc2¢T , 2 
[eet tomes (t. 2)) C(L.2)=C3(1,2 
ihe Genes meee, C4(1,2)) C(1,2)=C4(1,2 
L(t)=¢(1,2)-C(T, 1) 
fr (L(1) le 0:01 )0 THEN 
CENTER=(C(I,2)+C(I,1))/2.0 
@ £2 )=CENTERTO. 005 
C(I, 1)=CENTER-0.005 
L(1)=0.01 
ENDIF 


20 CONTINUE 
Were ee") ‘SUBROUTINE COMP END' 


RET 
END 
os  .  —_ ee, 
Cx** SUBROUTINE COMPAR 
Cx 
C+++ 
cttt THIS SUBROUTINE COMPARES THE MIN. & MAX. X, Y & Z VALUES 
++ 
SUBROUTINE COMPAR(N Cae) 
REAL c(3, 2),X(N),Y(N),Z(N) 
DO 10 I=1,N 
Tac 4 .GT. X(I)) C(1,1)=X(I 
IF (C(1,2) .LT. X(I)) C(1,2)=X(T 
inet 2,1) .GT. ¥irt)) €(2,1)=V(I 
Meee( 2,2) LT. Y(1)) C(2,2)=v(I 
aici 3.1) .GT, Z(1T)} c(3,1)=Z(1 
IF (c(3,2) .LT. Z(1)) C(3,23=Z(1 
10 CONTINUE 
RETURN 
END 
gg | | << 9 
C*** SUBROUTINE CYLBDY 
Cx 
C+++ 


fae tfl5 SUBROUTINE COMPUTES THE POSITION OF THE POINT ICBDY 
C+++ NEEDED TO DRAW THE TWO LINES JOINING THE ELLIPSES 


SUBROUTINE SEDI CYL, MGYLP, YG@vEr, gevLP, ICBDY) 
PARAMETER (NCX=100) 
PARAMETER (NCECP= ECP) 
REAL XCYLP(NCX,NCECP) , YCYLP(NCX,NCECP) , ZCYLP (NCX,NCECP) 
REAL ROBS(3,3),COM 
REAL RCCAPP (NCR) ECCAEP (NC), eee, 
REAL KCCAPN(NCX) , YCCAPN(NCX) , ZCCAPN(NCX 
REAL X(NCECP), Y(NCECP), Z(NCECP) 
INTEGER ICBDY (NCX 
COMMON /CGDATA/ ICG,COMPCG,XCCAPP,YCCAPP,ZCCAPP, 
x XCCAPN, YCCAPN, ZCCAPN 
COMMON /OBDATA/ XOBS, YOBS, ZOBS, ROBS 
COMMON /PIDATA/ PI, RTD 
COMMON /VUDATA/ PHI, THETA, RADIUS 
WRITE (2, 
WRITE (2, ve SUBROUTINE CYEBDY* 
‘5 WRITE De LCYL =a Lcye 
ttt 


C+++ TRANSFER THE LOCUS OF THE POSITIVE END CAP TO 
C+++ OBSERVATION CO-ORDINATE 


C+t++ 
DO 10 I=1,NCECP 
oo eas 
1) =YCYLP( ICY ,I 
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pie eGR CUCL I) 
10 CONTINUE 

CALL ROTRAN(NCECP ,X,Y,Z,XOBS,YOBS ,ZOBS , ROBS) 
C+++ 


C+++ TRANSFER THE CENTERS OF THE TWO END CAPS TO 
C+++ OBSERVATION CO-ORDINATE 


C+t++ 

XCP=XCCAPP(ICYL 

YCP=YCCAPP(ICYL 

ZCP=ZCCAPP( ICYL 

XCN=XCCAPN(ICYL 

YCN=YCCAPN (ICYL 

ZCN=ZCCAPN( ICYL 

CALL OUR cr PTGHT Sani oe cir GEC eye 
5 CALL ROTRAN(1,XCN,YCN,2ZCN,XOBS,YOBS ,ZOBS,ROBS 
att 


C+++ THE PLANE FORM BY X(OBS)-AXIS AND CEs IS THE PLANE 
C+++ THAT WILL BE DISPLAYED ON THE SCREEN 

C+++ HENCE, WE CAN REMOVE THE Z(OBS) COMPONENTS AND WORK ONLY WITH 
ei X(OBS$ AND Y(OBS) COMPONENTS. 

+++ 

C+++ THE POINT THAT THE LINES (BODY OF THE CYLINDER) MEET THE ELLIPSES 
C+++ OF THE END CAPS WILL HAVE MAXIMUM PERPENDICULAR DISTANCE H 
C+++ FROM THE LINE JOINING THE CENTERS OF THE TWO END CAPS. 

Cio 

C+++ ONLY ONE POINT (ICBDY) NEED TO BE COMPUTED BECAUSE THE OTHER 
C+++ POINTS ARE SYMMETRICAL TO IT AND CAN BE COMPUTED EASILY 


C+++ 
H1=0 
D=S Ree CN ats RCN cee YCN)*(YCP-YCN) ) 
IF 0.00001) GOTO 
DO et NCECE 
Dl= =S0RT (4 eee Be Ree ey my eee ey agen IG > 
Pen Be Rte I )-XCN)*(X(I)-XCN)+(Y(I)-YCN)*(Y(I) -YCN 
(D1 .LT. 0.000001) GOTO 20 
DUMMY=(BLeb1+D*D- D2*D2)/2.0/D/D1 
er aecCnne 2GT.. is 0001) GOTO 70 
IF (ABS( DUMMY GT. 1.0 .AND. ABS(DUMMY) .LE. 1.0001) THEN 
DUMMY=SIGN(1.0, DUMMY) 

END IF 
ALPHA=ACOS (DUMMY ) 
H3=D1*SIN(ALPHA) 
IF 2 -LE. 2) GOTO 30 
IF (H2 .GE. Hl .AND. H2 .GE. H3) GOTO 60 
GOTO 30 

20 CONTINUE 
H3=0.000001 

30 CONTINUE 
H1=H2 
H2=H3 

40 CONTINUE 
GOTO 60 


50 CONTINUE 
ICBDY(ICYL)=1 

60 CONTINUE 
en GG 
WRITE (2,*) 'SUBROUTINE CYLBDY END' 
RETURN 

70 PRINT *, 'ARGUMENT FOR ACR COS EXCEEDS LIMIT OF +/- Q! 


tae ae | peg a ae COS EXCEEDS +/- it 
WRAL (2, “a! Dis =" 7pi 

WRITE (2,*) ' BZ. =" Dz 

WRITE. (255) * ARG =' ,DUMMY 

CALL DONEPL 

STOP 

END 


CARA 
CAR 
CAR 


Ct 
erat 
erat 
C+++ 
ertt 
C+++ 
C+++ 
C+++ 
rat 


Crt + 
erst 
Saas 


COMMON 
* 


SUBROUTINE CYLCAP 


THIS SUBROUTINE COMPUTES THE POINTS NEEDED TO DRAW THE ELLIPSES 


ANGLES NEEDED TO DEFINE THE CYLINDER CO-ORDINATE 
THETA TO ZI AXIS : ICLZ 


PHI TO ZT AXIS. : PCLZ 
THETA TO XT AXIS : TCLX 
PHI TO XT AXIS: PCLX 


SUBROUTINE Saree meine ZCYLP, XCYLN, YCYLN, ZCYLN) 
PARAMETER eae =100) 
PARAMETER (NCECP=201) 
Boat, COMPEG(3,2),ROT(3,3),R(3,3) 
REAL RCELE (NCK ,NCECE), ,YCYLP NCX, eee a EYEN NCR. Reece 
REAL ACYLN(NCX ,NCECP) , YCYLN(NCX ,NCECP) , ZCYLN(NCX ,NCECP 
REAL eee CEI SCCAEP(NCK), ee NCS) 
REAL XCCAPN(NCX) , YCCAPN(NCA) , ZCCAPN (NCX 
REAL XP NCECE) «20 (NCECE) -ZEINCECP 
REAL XN(NCECP) , YN(NCECP) ,ZN(NCECP 
CGDATA/ ICG,COMPCG,XCCAPP ,YCCAPP,ZCCAPP, 
XCCAPN, YCCAPN, ZCCAPN 
COMMON /PIDATA/ PI, RTD 
COMMON /RTDATA/ XCR, YCR, 2CR, ROL 
COMMON 7 aid FREQ, IUN, IUS, sCALEF 


WRITE 

WRITE (2,*) ‘SUBROUTINE CYLCAP'! 

WRITE (2 a ICG =' ICG 

PEAD (1.9%) XCCYL, YCCYL,ZCCYL 

READ (1,.*) TCLZ,PCLZ,TCLX, PCLY 

READ (1,% 

READ (1,% ZCN, THTN,ZCP,THTP 

WRITE (2, XCCYL =', XCCYL 

WRITE (2,*) | YCCYL =! YCCYL 

WRITE %2,*) | ZCCYL =! .ZCCYL 

WRITE (2,*) ! Teuze—' TenzZ.§ PELZ =' PCLZ 
WRITE (2,%*) ! TCLX =',TCLX,' PCLX =". PCLX 
WRITE (2,%*) ! AC =' ac ' BC =! BC 

WRITE (2,*) ! ZCN =!,ZCN,' ZCP =',ZCP 
WRITE (2,%*) ! THTN =',THIN,' THTP =',THTP 


COMPUTE THE LOCI OF THE ELLIPSES IN CYLINDER CO-ORDINATE 


NHALF=(NCECP+1) /2 
RTHP=THTP/RTD 
RTHN=THTN/RTD 
XSTEP=2.0*AC/(NHALF-1) 
XSTART=-1.0*AC-XSTEP 
DO 10 I=1,NHALF 
ayes =XSTART+REAL (1) *KSTEP 
ae BP 00, 90-0 ee” 
IF (THIP 

UEC) =2 


E 
ZP(I)=ZCP+XP(1I)/TAN(RTHP) 

END IF 

oe a 

YN(I)=YP(I 

If REN a 90.0) THEN 
ZN(I)=ZC 


5 2N(1)=2CN+XN (1) /TAN(RTHN) 
IF 

XP (NCECP-I+1)=XP(I 

YP (NCECP-I+1)=-1.0*YP(I) 


AmMiGeCP=I+1 )=ZP( 1 
XN (NCECP-I+1)=XN(I 


Sl 


eae Roe ot) 
Z2N(NCECP-I+1 )=ZN(I) 


10 CONTINUE 


Cis 
Cres 
Cicer 


Crt: 


C+++ 
Cs 


C+++ 


YN (NCECP )=YN(1 
YP(NHALF )=-1.0 
YN(NHALF )=-1.O*YN 


TP (NCECP )=VEP 1} 
Te 


eae 
NHALF 


COMPUTE THE AXES ROTATION MATRIX FOR CYLINDER CO-ORDINATE 


CALL ROTATE(TCLZ,PCLZ,TCLX,PCLX,R) 


TRANSFER THE POINTS TO GLOBAL CO-ORDINATE 


CALL ROTRAN MRE eT oH ot REE Hee Mee 
CALL ROTRAN(NCECP,XN,YN,2ZN,XCCYL, YCCYL, ZCCYL,R 
CALL ROTRAN Nees aN HAY ac eee 
CALL ROTRAN(NCECP, XN, YN,2ZN,XCR, YCR, ZCR, ROT 
XCCAPP (ICG 
YCCAPP (ICG )=0 
ZCCAPP (ICG )=ZCP 
XCCAPN (ICG )=0 
YCCAPN (ICG )=0 
ZCCAPN( ICG CN 
CALL ROTRAN(1,XCCAPP(ICG) , YCCAPP(ICG) ,ZCCAPP(ICG), 
XCCYL , eC iL, ZCCuE R 


) 
4 CALL ROTRAN(1 ,XCCAPN (TCG) ,YCCAPN(ICG), ZCCAPN(ICG), 


Wi cecrl, COCR 
CALL Sea ee Keener ye) cae Bearer t cay, ee soe 
CALL ROTRAN(1,XCCAPN( ICG) , YCCAPN( ICG) ,ZCCAPN( ICG) ,XCR, YCR, ZCR, ROT 


Ctt+ SCALE THE Forts 


C+++ 


ZCYLN 
20 CONTINUE 


C+++ 


CALL SCALES(NGECE {G7 Yee 
CALL SCALES (NCECP, 'G' XN pL NG ZN 
CALL SCALES (13. G.; ROCABP | ICG) sae ae 


CALL SCALES co XCCAPN (ICG) , YCCAPN( ICG) ,ZCCAPN( ICG 
DO 20 I=1,NCECP 
ACYLP (LCG) —Ab ‘ 


YCYEP( (CG. = 7a 
Z2CYLP(ICG,I)=2P (1 
ACYLN(ICG,I)=AN(1 
YCYLN(ICG,I)=YN(I 
ICG,I)=2N(1 


Sie COMPUTE THE MIN. & MAX. X, Y & 2 VALUES FOR CYLINDER GEOMETRY 
eet 


CALL eee esc CONREG, ARy we ae 
CALL COMPAR(NCECP ,COMPCG,A&N, LYN, aN 
WRITE (2,*) ‘SUBROUTINE CYLCAP END! 
RETURN 

END 


SUBROUTINE DISSET 


THIS SUBROUTINE SETS UP THE GRAPHIC DIMENSIONS, AAESSAND 011 


SUBROUTINE DISSET(COMPAL ,LENGTH, TITLE) 
REAL COMPAL(3,2), LENGTH (3 

CHARACTER*7 PHIS, THETAS ,RADIS 
CHARACTER*60 Tl ney LABEL1, LABELZ , LABEL3 
COMMON /VUDATA/ PHI, THETA ,RADIUS 

WRITE (2, 
VRE CZ; , PIUBReULINES Dosa 
OPEN (3, STATUS='NEW' 

WRITE (3, 100) PHI,THETA, RADIUS 
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100 FORMAT (1X, F7.2,1X,F7.2,1X,F7.2) 
CLOSE (3) 


OPEN (3, STATUS='OLD' ) 
READ (3,200) PHI$,THETAS,RADIS 

200 FORMAT(1X,A7,1X,A7,1X,A7 
LABELI='VIEWING POINT : ?F) = '//PHI$//' DEGS' 
LABEL2=' 20) = HETAS//' DEGS' 
LABEL3=' ()R = '//RADI$//' M$! 


CLOSE (3) 

CALL AREA2D(8.0,8.0) 

CALL VOLM3D(LENGTH(1),LENGTH(2) ,LENGTH(3) ) 

CALL HEADIN(TITLE,100,1.0,4) 

CALL HEADIN(LABEL1,100,1.0,4 

CALL HEADIN(LABEL2,100. i. 0/4 

CALL HEADIN(LABEL3,100,1.0, "4 

CALL X3NAME('(X) IN mgt 100 

CALL Y3NAME('(Y) IN MS'‘100 

CALL Z3NAME('(Z) IN MS'.100 

THETAV=90.0-THETA 

CALL VUANGL ( PHI -THETAV , RADIUS) 

CALL GRAF3D(COMPAL(1,1), 'SCALE' EOL, 2), 
et. A. 'SCALE' , COMPAL (2,2), COMPAL(3,1), 'SCALE' ,COMPAL(3,2)) 


WRITE (2,%) 'SUBROUTINE DISSET END! 
RETURN 
END 
ES ne ea 
Cxxx SUBROUTINE DISSIN 
CxKX 
C+++ 
CH THIS SUBROUTINE INITIALIZE DISSPLA 
+4+ 
SUBROUTINE DISSIN 
CALL COMPRS 
CALL PAGE(8.5,11.0) 2 
CALL FUTURA 
CALL BASALF('L/CSTD') 
CALL MIXALF('STANDARD! ) 
CALL MX3ALF('L/CGREEK', '?') 
CALL NOBRDR 
CALL YAXANG(0.0) 
RETURN 
END 
preeroses 9 fll ees 
Cxx* SUBROUTINE INVERS 
CxKK 
C+t+ 


C+++ THIS SUBROUTINE INVERSES THE 3 X 3 MATRIX RIN 


SUBROUTINE Benen ROUT ) 
REAL RIN(3,3),R 


1 ope) 

DETERM=RIN(1, 1 SRIN Zee ori 5) Paine, 2) *RIN(Z,3)*RIN( 351 )+ 
os RIN 2, '3)*RIN Zouk ee ZR ENGL, '1)*RIN 2 '3)*RIN(3 12 = 
* RIN(1,2)*RIN(2,1)*RIN(3,3)-RIN(1,3)*RIN(2,2)*RIN 3; i 

IF (DETERM .EQ. 0.0) GOTO 10 

CU, 1 )=(RINCZ, 2 ein 3,3)-RIN(2,3)*RIN(3,2))/DETERM 

ROUT(1,2)=(RIN 1,3 *RIN(3,2)°RIN(1,2)*RIN(3,3) )/DETERM 

ROUT(1,3)=(RIN(1,2)*RIN(2,3)-RIN(1,3)*RIN(2,2))/DETERM 

ROUT(2,1)=(RIN(2Z,3)*RIN(3,1)-RIN(2,1)*RIN(3,3))/DETERM 

ROUT(2,2)=(RIN(1,1)*RIN(3,3)-RIN(1,3)*RIN(3,1))/DETERM 

ROUT(2,3)=(RIN(1,3)*RIN(2,1)-RIN(1,1)*RIN(2,3))/DETERM 

ROUT(3,1)=(RIN(2,1)*RIN(3,2)-RIN(2,2)*RIN(3,1))/DETERM 

ROUT(3,2)=(RIN(1,2)*RIN(3,1)-RIN(1,1)*RIN(3,2) ) /DETERM 

See =(RIN(1,1)*RIN(2,2)-RIN(1,2)*RIN(2,1))/DETERM 


10 CONTINUE 
PRINT *, 'THE DETERMINANT OF RIN IS EQUAL TO ZERO! 
WRITE en *) "SUBROUTINE INVERS' 
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WRITE (2,%*) ' THE DETERMINANT OF RIN IS EQUAL TO ZERO' 
WRITE “C2553 THE MATRIX RIN HAS THE FOLLOWING VALUES! 
WRITE (2,%) ' ', RING] ,1),° “©, RONC D2 eee 
WRITE (2,*) | ', RIN(2Z,1),° *, RENR Zee) eee 
WRETEGZa) ' RIN(S,1),' ',REING@32Z) eer 
CALL DONEPL 

SLOP 

END 


SUBROUTINE NOCYL 


THIS SUBROUTINE RESETS ALL DATA BELONG TO CYLINDER GEOMETRY 


SUBROUTINE NOCYL(XCYLP,YCYLP,ZCYLP, XCYLN, YCYLN, ZCYLN) 
PARAMETER Neate 

PARAMETER (NCECP=201) 

REAL COMPCG(3,2),R(3,3) 

REAL RCYLP (NC NCECE |, 22? (NCE CECE es 
REAL XCYLN(NCX,NCECP) , YCYLN(NCX,NCECP) , ZCYLN(NCX,NCECP 
REAL SEE en eer ee eee ted 
REAL XCCAPN(NCX) , YCCAPN(NCX) , ZCCAPN(NCX 


4g, COMMON /CGDATA/ ICG,COMPCG,XCCAPP,YCCAPP,ZCCAPP, 


ACCAPN , YCCAPN , ZCCAPN 
GALL REESE Gt TNR Neon Geetn eeu comiee nines 
CALL RESET1(1,NCX,;NCECP, XCYLN , YCYLN (2ZGVEh Gereec 
CALL RESEEZ MERE CARH NCH UTCCR DHT Oa aaemy 
CALL RESET2(NCX, ACCAPN ,NCX, YCCAPN ,NCX, ZCCAPN 
ae 
E 


SUBROUTINE NAT 


THIS SUBROUTINE RESET ALL GEOMETRY DATA 


SUBROUTINE NXT(XPC,YPC,ZPC,XCYLP, YCYLP,ZEYLP , XCYENP ye’ Lieae vie 


“XSC,YSC ,2SGj45a6, YSAC, ZoAC, LGP | Tair) 


PARAMETER (NPX=50 
PARAMETER (NEX=12 
PARAMETER (NCK=100) 
PARAMETER (NCECP=201) 
PARAMETER (NSX=1000) 
PARAMETER (NSAX=5 
PARAMETER (NESAX=500) 

REAL COMPPG(3,2),COMPCG(3,2) , COMPSG(3,2) , COMPSA(3, 2) 

REAL ROT(3,3),ROBS(3,3) 

REAL RCELP INCH NCECP} , YCYLP (NCR NGECE) ,ZOuLE (Ne Ncoge) 

REAL XCYLN(NCX,NCECP) . YCYLN(NCX.NCECP) , ZCYLN(NCX, NCECP 

REAL KCCAPE (NCK) , ZCCAPP (NCK) ,ZCCARE NCX 

REAL XCCAPN(NCX) , YCCAPN(NCX) , ZCCAPN(NCX 

REAL XPC(NPX,NEX) , YPC(NPX,NEX) , ZPC(NPX, NEX) 

REAL XSC(NSX,4) , YSC(NSX,4) ,ZSC(NSX, 4 

REAL XSAC(NSAX, NESAX, 4) , YSAC(NSAX, NESAX, 4) , ZSAC(NSAX,NESAX, 4) 
INTEGER MEP(NPX) ,MSAX(NSAX) 

CHARACTER*60 TITLE 

LOGICAL LGP 


4, COMMON /CGDATA/ ICG,COMPCG,XCCAPP ,YCCAPP,ZCCAPP, 
ane 


CAPN , YCCAPN, ZCCAPN 
COMMON /OBDATA/ XOBS,YOBS,ZOBS,ROBS 
COMMON /PGDATA/ IPG,MEP,COMPPG 
COMMON /RTDATA/ XCR,YCR,ZCR,ROT 
COMMON /SADATA/ ISA,MSAX,COMPSA ,MSAXT 
COMMON /SCALER/ FREQ, IUN,IUS,SCALEF 
COMMON /SGDATA/ ISG,COMPSG 
COMMON /VUDATA/ PHI,THETA, RADIUS 
CALL RESET1(1,NPX,NEX,XPC,YPC,ZPC, COMPPG) 
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CALL RESET1 eee 2 2CeL? .YCILP.ZC¥LF  COMPCG) 
CALL RESET1(1,NCX,NCECP , XCYLN, YCYLN, ZCYLN , COMPCG 
CALL RESET1(1,NSX,4,XSC,YSC,ZSC, COMPSG) 

CALL RESET1(NSAX ,NESAX,4,XSAC,YSAC,ZSAC, COMPSA) 
CALL RESET2(NPX ,MEP,NSAX,MSAX,1, FREQ) 

CALL RESET2(1,XCR,1,YCR,1,ZCR) 
CALL RESET2(1,X0OBS,1,YOBS,1,ZOBS) 

CALL RESET2 NCK, RCCAPP ,NCK, ZCCAPE ,NCX, ZCCAPP} 
CALL RESET2(NCX , XCCAPN , NCX, YCCAPN ,NCX, ZCCAPN 
CALL ROTATE{0.0,0.0,90.0,0.0,ROT) 

CALL ROTATE(0.0,0.0,90.0,0.0,ROBS) 


LGP=.FALSE. 
TITLE='$' 
SCALEF=1.0 
PHI=60 .0 
THETA=30.0 
RADIUS=1.0 
ICG=0 
IPG=0 
ISA=0 
ISG=0 
LUN=1 
IUS=1 
RETURN 
END 
CxARK 
ons SUBROUTINE OBSERV 
Ec 
or 


C+++ THIS SUBROUTINE COMPUTES THE AXES ROTATION MATRIX AND 
C+++ ORIGIN TRANSLATION CO-ORDINATES FOR THE CONVERSION FROM 
C+++ GLOBAL CO-ORDINATE TO OBSERVATION CO-ORDINATE 


Sao Ll ASSUMES THAT THE CENTER OF OBSERVATION IS AT THE CENTER 
Ctt+ OF THE WORK BOX DEFINED BY DISSPLA 


SUBROUTINE OBSERV(COMPAL) 

* REAL COMPAL(3,2),ROBS(3,3),R(3,3) 
CHARACTER ANSS 
CHARACTER*2 THETAS,PHIS 
CHARACTER*20 MSSG1,MSSG2 
COMMON /OBDATA/ XOBS, YOBS ,ZOBS,ROBS 
COMMON /OSDATA/ POSMAX, TOSMAX, PHIOS , THETOS 
COMMON /PIDATA/ PI,RTD 
COMMON /VUDATA/ PHI,THETA,RADIUS 
WRITE (2.3) 

WRITE (2,*) ‘SUBROUTINE OBSERV'! 
C+++ 


iret lhe OFFSET ANGLES ARE NEEDED TO COMPENSATE FOR THE DIFFERENCES 
C+++ BETWEEN THE VIEWING ANGLES INTERPRETED BY THE PROGRAM AND THAT 
St SEEN BY DISSPLA 





C+++ 
Be ESE AUG IOS COstnL CIN OREGRI EDGR 
‘ PHIOS=REAL (INT(-1.0*POSMAX*SIN(THETA/RTD) ) ) 
+++ 
C+++ DETERMINE THE CENTRE OF THE WORK BOX 
C+++ 
XCWB=(COMPAL(1,2)+COMPAL(1,1))/2.0 
YCWB=(COMPAL(2,2)+COMPAL(2,1))/2.0 
: ZCWB=(COMPAL(3,2)+COMPAL(3,1))/2.0 
+++ 


Soot COMPUTE THE OBSERVATION POINT WITH RESPECT TO THE CENTER 
Grt+ OF THE WORK BOX 


C+t+ 
RPHI=(PHI+PHIOS)/RTD 
RTHETA=(THETA+THETOS )/RTD 
Oe eae on cesar een acne) 
YCOBS=RADIUS*SIN(RTHETA) *SIN(RPHI 
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C++t 
C+++ 
C+++ 


Ct+++ 
Cra 
Crt 
C+t++ 
C+++ 


C+++ 
C+++ 
C+++ 


C+++ 
C+++ 
C+++ 


C+++ 
C+++ 
C+++ 


ZCOBS=RADIUS*COS (RTHETA) 
DETERMINE THE OBSERVATION POINT IN GLOBAL CO-ORDINATE 


XOBS=XCWB+XCOBS 
YOBS=YCWBt+YCOBS 
ZOBS=ZCWB+ZCOBS 


DETERMINE THE NEW AXES ANGLES CENTERED AT THE OBSERVATION 
POINT WITH RESPECT TO THE GLOBAL CO-ORDINATE WITH THE 
Z(OBS) AXIS POINTING TOWARDS THE CENTER OF THE WORK BOX 


THETAZ=180 .O- (THETA+THETOS ) 

PHIZ=180. 0+ (PRIFEHTOS) 

IF (PHIZ .GT. 360.0) PHIZ=PHIZ-360.0 
THETAX=ABS (90. O-THETAZ ) 

IF (THETAZ ~0G. Gr i PHIX=180.0+PHIZ 
TE (THETAZ ~CER. 7070) PHIA=PHIZ 

IF (PHIX .GT. 360. 0) PHIX=PHIX-360.0 


FORM AXES ROTATION MATRIX 
PHIZ, THETAX, PHIX,R) 


CALL igueRS(n Boes 

CALL INVERS BS} 

WRITE (2,% THE MATRIX ROBS HAS THE FOLLOWING VALUES' 

WRITE (2.%*) | | ROBS(1,1),' ",ROBS( ie aye =" RoBs Tum 

WRITE (2.*) | '” ROBS(2.1 | ROBS(2,2),' ‘' ROBS(2,3 

WRITE (2.*) ! '” ROBS(3,1 '’ROBS(3,2),' ' ROBS(3,3 
COMPUTE THE ORIGIN TRANSLATION 

XOBS=-1.0*XOBS 

YOBS=-1.0*YOBS 

ZOBS=-1.0*ZOBS 

CALL ROTRAN(1, XOBS , YOBS ,ZOBS ,0.0,0.0,0.0,ROBS) 

WRITE (2,% THE ORIGIN TRANSLATION HAS THE FOLLOWING VALUES! 

WRITE (2.%*) | XOBS =', XOBS 

WRITE (2.%*) | YOBS =',YOBS 

WRITE (2,*) |! ZOBS =' ZOBS 

WRITE (2,*) 'SUBROUTINE OBSERV END' 

RETURN 

END 


a gqe@g@qgq@@ @@eqgee@newmea.e@e @ @&e ee eg se SF eet eeGgtaegG:@e@ge@e@aeqgese Geeeeqgeqg@egeqq@geqeaqqg@@@@geqgqq@ @ @& 


SUBROUTINE PLATEG 


THIS SUBROUTINE READS IN THE PLATE CORNERS AND TRANSFER THEM TO 
GLOBAL CO-ORDINATE 


SUBROUTINE PLATEG( KPC, YPC,ZPC) 
PARAMETER (NEH=50 

PARAMETER (NEX=12 

REAL XPC(NPX, NEX) , YPC(NPX,NEX), ,2PC(NPX, NEX) ,X(NEX) , ¥(NEX) ,Z(NEX) 
REAL COMPPG(3,2),R(3,3),ROT(3,3) 
INTEGER MEP NBS 

COMMON /PGDATA/ IPG,MEP ,COMPPG 
COMMON /RTDATA/ XCR,YCR,ZCR,ROT 
COMMON /SCALER/ FREQ, IUN, IUS, SCALEF 
WRITE o x 

WRITE (2,*) ‘SUBROUTINE PLATEG' 


READ IN DATA 
Ren oy Teun Bere 


WRITE (2,*) P(IPG) =' MEP(IPG) 
IF (HER ({P6) Ache NEX) GOTO 50 
IF (LSLAB .EQ. 0) G 20 


READ(1,*) NSLAB 
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ee 


DO 10 I=1,NSLAB 

READ(1,*) DUMMY,DUMMY, DUMMY , DUMMY , DUMMY 
10 CONTINUE 
20 CONTINUE 


DO 30 Seah 
REaw ls) x) Y(h) (1) 
30 CONTINUE 
Crt 
ae TRANSFER THE POINTS TO GLOBAL CO-ORDINATE 
+++ 
: CALL ROTRAN(MEP( IPG) ,X,Y,Z,XCR,YCR,ZCR, ROT) 
+++ 
Gitte oCALE THE POINTS 
C+++ 
CALL ES ieart tee AOC cONe Ges 
CALL COMPAR(MEP( IPG) ,COMPPG,X,Y,Z) 
WRITE (2,*) ' EeG = Ieee 
DO 40 I=1,MEP(IPG) 
REC( IPG, 1I)=X(1 
YPC(IPG,I)=Y(1 
ZPC(IPG,I)=Z(1 





40 CONTINUE 
WRITE (2,*) 'SUBROUTINE PLATEG END' 
RETURN 


50 PRINT *, 'NO. OF CORNERS EXCEEDS NEX' 
CALL DONEPL 
“S20P 


END 
ChRK ee e@e@eee@eeee@ @ ee @ ee@@oege=aqrFee@eeeee@ere@g@eeeaee@e@gee@eerE@aee@@ ee eB Bets FFB SF FSF FSF FB se ae a2 = = 
a SUBROUTINE PLOT 
C 
C+++ 


C+++ THIS SUBROUTINE PLOTS THE GEOMETRY OF THE INPUT DATA 


SUBROUTINE PLOT(XPC,YPC,ZPC,XCYLP, YCYLP , ZCYLP , XCYLN, YCYLN, ZCYLN, 
*XSC,YSC,ZSC, XSAC, YSAC, ZSAC, LGP, TITLE) 

PARAMETER (NPX=50 

PARAMETER (NEX=12 

PARAMETER (NCX=100) 

PARAMETER (NCECP=201) 

PARAMETER (NSX=1000) 

PARAMETER (NSAX=5) 

PARAMETER (NESAX=500) 

REAL XPC(NPX,NEX) , YPC(NPX,NEX) ,ZPC(NPX,NEX) 

REAL re ne oe) DLE (HCA TEE) CELE (NCH, NCECP } 

REAL XCYLN(NCX,NCECP) , YCYLN(NCX,NCECP) , ZCYLN(NCX,NCECP 

REAL este nes NC 

REAL XCCAPN(NCX) , YCCAPN(NCX) , ZCCAPN(NCX 

REAL XSC(NSX,4) , YSC(NSX,4) ,ZSC(NSX, 4 

REAL XSAC(NSAX,NESAX,4) , YSAC(NSAX, NESAX Pee a eoee 4) 

REAL COMPAL(3,2),COMPCG(2,3),COMPPG(3,2),COMPSA(3,2) , COMPSG(3, 2) 

REAL ROT(3,3),ROBS(3,3 

REAL LENGTH(3 

INTEGER MEP(NPX) , MSAX(NSAX) 

CHARACTER*60 TITLE 

CHARACTER ANSS 

LOGICAL LGP 

COMMON /CGDATA/ ICG,COMPCG,XCCAPP,YCCAPP ,ZCCAPP, 
* XCCAPN , YCCAPN, ZCCAPN 

COMMON /OBDATA/ XOBS,YOBS,ZOBS,ROBS 

COMMON /OSDATA/ POSMAX, TOSMAX, PHIOS, THETOS 

COMMON /PGDATA/ IPG,MEP,COMPPG 

COMMON /PIDATA/ PI,RTD 

COMMON /RTDATA/ XCR,YCR,2ZCR,ROT 

COMMON /SCALER/ FREOQ,IUN,IUS,SCALEF 

COMMON /SGDATA/ ISG, COMPSG 

COMMON /SADATA/ ISA,MSAX,COMPSA ,MSAXT 

COMMON /VUDATA/ PHI ,THETA, RADIUS 
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DATA COMPAL/3%*10. 0E10,3*-10.0E10/ 
WRITE ae 
WRITE (2,*) ‘SUBROUTINE PLOT! 

DETERMINE OVERALL MIN. & MAX. X, Y & Z VALUES 


CALL COMP( COMPAL ,COMPCG,COMPPG,COMPSA,COMPSG, LENGTH) 


10 CONTINUE 


C+++ 
C+++ 
C+++ 


C+++ 
C+++ 
C+++ 


C+++ 
Cah 
C+++ 
Crete 


C+++ 
C+++ 
C+++ 


DETERMINE THE VIEWING ANGLES 
CALL VANGLE (TITLE, LENGTH) 

SET UP DISSPEASRCR S40 oF LOLTiic 
CALL DISSET(COMPAL ,LENGTH, TITLE) 


COMPUTE ROTATION MATRIX AND TRANSLATION POSITION FOR 
OBSERVATION CO-ORDINATE 


CALL OBSERV(COMPAL) 
PLOT INFINITE GROUND PLANE ? 


IF (.NOT. LGP) GOTO 20 
CALL PLOTGP (COMPAL) 


20 CONTINUE 


C+++ 
C+++ 
C+++ 


PLOT PLATE(S) ? 


(IPG EQ, 0) GOTO 30 
aie PLOTPG(XPC, YPC, ZPC) 


30 CONTINUE 


Ct+t++ 
C+t++ 
C+++ 


PLOT CYLINDER(S) ? 


IF (ICG .EQ. 0) GOTO 40 
CALL PLOTCG(XCYLP, YCYLP, ZCYLP, XCYLN, YCYLN, ZCYLN) 


40 CONTINUE 


Cr 
Cat 
Cres 


PLOT SOURCE(S) ? 


IF (ISG .EQ. 0) GOTO 50 
CALL PLOTSG(XSC, YSCsZsc) 


50 CONTINUE 


C+++ 
C+++ 
Girt 


PLOT SOURCE ARRAY(S) ? 


IF (ISA .EQ. 0) GOTO 60 
CALL PLOTSA(XSAC, YSAC,ZSAC) 


60 CONTINUE 


CALL ENDPL(0) 


70 CONTINUE 


PRINT *, 'DO ae WANT ANOTHER VIEW OF THE PLOT (Y OR N)' 
NS 


READ(5,'(Al)') ANSS 

IF (ANSS .NE. 'Y' .AND. ANS$ .NE. 'N') GOTO 70 
IF (ANSS .EQ. 'N') GOTO 80 

GOTO 10 


80 CONTINUE 


wee ale) ‘SUBROUTINE PLOT END' 
U 


SUBROUTINE PLOTCG 


THIS SUBROUTINE PLOTS THE CYLINDER(S) 
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SUBROUTINE PLO EARCYLE,, VevEr, ZCYLP Xey ih vey LN, ZCYLN) 
PARAMETER eee ) 
PARAMETER (NCECP=201) 
REAL ROT(3,3),ROBS(3,3) ,COMPCG(3,2) 
REAL Cece NCECP), ET CICEEENE SooucL NEES 
REAL XCYLN(NCK ,NCECP) , YCYLN(NCX,NCECP) ,ZCYLN(NCX ,NCECP 
Repiben2(2),%2\2),22 
REAL RCCAEE (NCK /ECCABP (NCK) , Been 
REAL XCCAPN(NCX) , YCCAPN(NCX) ,ZCCAPN(NCX 
INTEGER ICBDY (NCX) 
LOGICAL LBLNK1,LBLNK2 
COMMON /CGDATA/ TeG COHPCG, ACCAP EM GCAPP,ZCCAFPP, 

* XCCAPN, YCCAPN , ZCCAPN 

COMMON /OBDATA/ XOBS, YOBS, ZOBS, ROBS 
COMMON /PIDATA/ PI, RTD 
COMMON /RTDATA/ XCR, Wek, ZCR KOT 
COMMON /SCALER/ FREQ, IUN, IUS, SCALEF 
COMMON /VUDATA/ sg ol THETA, RADIUS 
WRLEIE (2, x 
WRELE (2, /*) ‘SUBROUTINE PLOTCG'! 
DO 10 I= i, ICG 

C+++ 


Soe EL LOL THE LINES JOINING THE ELLIPSES 


C+++ 
CALL CYLBDY(I,XCYLP, YCYLP , ZCYLP, ICBDY) 
1)=XCYLP(I; ICBDY(I 
1)=YCYLP(I, ICBDY(T 
1)=ZCYLP(I, ICBDY(I 
X2(2)=XCYLN(I, ICBDY(I 
2\=YCYLN(I, ICBDY(I 
2)=ZCYLN(1. ICBDY(I 
CALL CURV3D(X2,Y2,22,2,0) 
NHALF=(NCECP-15/2. 
1)=XCYLP(I, ICBDY(I)+NHALF 
1)=YCYLP(I, ICBDY(1)+NHALF 
1)=ZCYLP (I, ICBDY(1)+NHALF 
X%2(2)=XCYLN(1, ICBDY(1I)+NHALF 
2\=YCYLN(1I, ICBDY(I)+NHALF 
2)=ZCYLN(1, ICBDY(I )+NHALE 
CALL CURV3D(X2,Y2,Z2,2,0) 
C+++ 


Sumte FLOL THE POSITIVE END CAP ELLIPSE 
C+++ BLANK OFF THE REGION THAT IS BLOCKED 


C+++ 
NHALF=(NCECP-1) /2 
CALL ae ieee cola) CUNEO na eee 
2 CALL PLOTEP(I,XCYLP,YCYLP,ZCYLP, ICBDY(I) ,LBLNK1,LBLNK2) 
+++ 


Sot © LOT THE NEGATIVE END CAP ELLIPSE 
C+++ BLANK OFF THE REGION THAT IS BLOCKED 


ae 
CALL ee GE ee CU CEU nERDe a LeLeL MRL ee 
CALL PLOTEP(I, XCYLN, YCYLN, ZCYLN, ICBDY(I) , LBLNK1, LBLNK2) 
10 CONTINUE 
eee! ‘SUBROUTINE PLOTCG END' 


ET 
END 
UT UCTEONSS DN tthe 
C*x*x SUBROUTINE PLOTEP 
CrKK 
ae 


Sars) THIS SUBROUTINE PLOTS THE END CAP ELLIPSE OF CYLINDER 


SUBROUTINE PLOTEP(TCYL, XCYL, YCYL, ZCYL, ICBDY , LBLNK1, LBLNK2) 
PARAMETER a 100 

PARAMETER (NCECP= aO1) 

REAL XCYL(NCX,NCECP) , YCYL(NCX,NCECP) ,ZCYL(NCX,NCECP) 

REAL X(NCECP) ,Y(NCECP) ,Z(NCECP) 
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LOGICAL LBLNK1,LBLNK2 


WIGELE <( 27.5 

WRITE (2,*) ‘SUBROUTINE PLOTEP' 
WRITE (Zo) LBLNK1=' , LBLNK1 
WRITE (2,*) ' LBLNK2=' , LBLNK2 


NHALF=(NCECP-1)/2 
IF (LBLNK1) GOTO 20 
C+++ 


C+++ PLOT THE ONLY FIRST BAER OF THE ELEIESE 


C+++ 
DO 10 J=1,NHALF+1 
WRITE (2,*) | DO LOOP 1! 
K=ICBDY+NHALF+J-1 
IF (K .GT. NCECP) K=K-NCECP+1 
X(J)=XCYL(ICYL,K 
¥(J)=YCYL 
Z(J)=ZCYL 
WRITE (2,*) J,° x0) 7° *¢Y (de eee cree 
10 CONTINUE 
CALL CURV3D(X,Y,Z,NHALF+1 ,0) 
GOTO 60 
20 CONTINUE 
IF (LBLNK2) GOTO 40 


C+t+ 
C+++ PLOT ONLY THE SECOND HALF OF THE ELLIPSE 
+e 

DO 30 J=1,NHALF+1 

WRITE 12%) 1 DO LOOP 2! 

K=ICBDY+J-1 


IF (K .GT. NCECP) K=K-NCECP+1 





X(J)=XCYL(ICYL,K 
v(J)=YCYL(ICYL,K 
Z(J)=ZCYL(ICYL,K 
WRITE. (2,*) J,' ',X(J),'. ', YCQ}R aes 
30 CONTINUE 
CALL CURV3D(X,Y,Z,NHALF+1,0) 
GOTO 60 
40 CONTINUE 
eae 
67-7 Ee PLODSIHED BULERE EE IPSE 
C+++ 
DO 50 J=1,NCECP 
WRITE (2,*) | DO LOOP 3! 
X(J)=XCYL(ICYL,J 
yY(JJ=YCYL(ICYL, J 
Z(J)=ZCYL(ICYL,J 
WRITE (27) Jp 8) XCS)5 ey ee econ 
50 CONTINUE 
CALL CURV3D(X,Y,Z,NCECP,0) 
60 CONTINUE 
WRITE (2,*) ‘SUBROUTINE PLOTEP END! 
RETURN 
END 
Cxx* 
Cx** SUBROUTINE PLOTGP 
Cx 
C+++ 


C+++ THIS SUBROUTINE PLOTS THE INFINITE GROUND PLANE 


SUBROUTINE PLOTGP ( COMPAL) 

REAL COMPAL (3,2) ,X(5),¥(5),2(5) 
WRITE (2, 

WRITE (2, eo ‘SUBROUTINE PLOTGP' 
1)=COMPAL(1,1)*1.2 
2\=COMPAL(1,1)*1.2 
3)=COMPAL(1,2)*1.2 
4)=COMPAL(1,2)*1.2 

5 )=X(1) 


6 PG OS OGG 
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alee 
mle 


lee 


IN NIN WN G&G FG EG EG 
a ale 
OD WDM OLB Whe 
il 
K; 
7“ 
i 
eee” 
NI AO BAO RO 


OOOO 


A>, 
WRITE (373) ee OF INFINITE GROUND PLANE (GLOBAL) ' 

I 20) De 2({I)' 
DO 10 I=1,5 


WRITE G5) Tench v1). sin 


10 CONTIN 


Shc CURV3D(X,Y.Z, 5,0) 

WRITE (2,*) ‘SUBROUTINE PLOTGP END! 
RETURN 

END 


SUBROUTINE PLOTPG 


THIS SUBROUTINE PLOTS THE PLATE(S) 


SUBROUTINE PLOTPG(XPC, YPC, ZPC) 
PARAMETER ee 

PARAMETER (NEX=12 

REAL XPC(NPX,NEX) , YPC(NPX,NEX) 
REAL X(NEX+1), ¥(NEX+1) , Z(NEX+1 5 
INTEGER MEP(NPX) 

COMMON /PGDATA/ IPG,MEP , COMPPG 
WRITE (2,% 

WRITE (2.% SUBROUTINE PLOTPG | 


PLATE NO. =',I 
CORNERS OF PLATE oo 
r'., X(T)! ¥(I)',' 2{I)' 


ZPC(NPX,NEX) 


= 

i 

tI 

Hj 

m4 
INK OX 
hdd Qe oO 


ANAIUNDNDN il 


Y(J)=Y 
M=MEP 

X(M)= 
Y(M)= 
a(M 


DO 20 
WRITE G x) K,X(K),¥(K),Z(R) 





20 CONTIN 


Cun CURV3D(X, Y,Z,M,0) 


30 CONTINUE 


WRITE (2,*) ‘SUBROUTINE PLOTPG END! 
RETURN 
END 


SUBROUTINE PLOTSA 


THIS SUBROUTINE PLOTS THE SOURCE ARRAY(S) 


SUBROUTINE EOS RSAC YSAC,ZSAC) 

PARAMETER NeeATe 

PARAMETER (NESAX=500) 

REAL ET ae NESAX,4) , YSAC(NSAX, NESAX,4) ,ZSAC(NSAX , NESAX, 4) 
REG) (5) a (5 

INTEGER MSAX (NESAX 


on 


10 


10 


COMMON 
WRITE 


ey ISA,MSAX,COMPSA,MSAXT 
can PLOTSA' 
a , SOURCE ARRAY NO. =' 


, SOURCE ELEMENT NO. ao 
CORNERS OF SOURCE ELEMENT (GLOBAL) ' 


ve 


= 
a 
Hi 
J 


1 NP II 
~~ js 


Z 
CONTINUE 
X(5)=K(1 
5)=¥(1 
2(5)=Z(1 
WRITE (2, Toe 
DO 20 
WRITE (2, ai K,X(K),¥(K),Z(R) 
CONTINUE 
CALL CURV3D(X,Y,Z,5,0) 
CONTINUE 
CONTINUE 
WRITE (2,%) 
URN 





eee. 








"2 Act) Gy, | EAC) eee Z(I)' 


YSUBROUTINE PLOTSA SEND 


SUBROUTINE PLOTSG 


THIS SUBROUTINE PLOTS THE SINGLE SOURCE(S) 


SUBROUTINE PLOTSG(XSC,YSC,ZSC) 
PARAMETER (NSX=1000) 

REAL XSC(NSX,4), YSC(NSK, 4) ,ZSC(NSX,4) 
REAL X(5),¥(5),Z 

COMMON /SGDATA/ 186, COMPSG 


WRITE 
"SUBROUTINE PLOTSG'! 


WRITE 
DO. SCrI 

SOURGE NO >= .7 21 
WRITE CORNERS OF SOURCE eee 
GREE raein a 
DO 10 J 
ACJ 


i 
1 
YJ =¥5¢\ 


+ OFF 


~ ses s Pw~wH 
W— 


Cy Cy Cy 
ee ase 


WRITE 
I I 


Y(l) Z(1)' 


eee 
MQ IW PODD IM 


i 

~*~ 

WN 
Hi 


Z(J)= 
CONTINUE 

X(5)=K(1 

v(S5)=V(1 

Z(5)=Z 1) 

DO 20 K 

WRITE (2. ty K,X(K),¥(K),Z(R) 


20 CONTINUE 


CALL CURV3D(X,Y,Z,5,0) 


30 CONTINUE 


WRITE (2,*) ‘SUBROUTINE PLOTSG END! 
URN 


SUBROUTINE RESETI 


THIS SUBROUTINE INITIALIZES THE GEOMETRICAL CO-ORDINATES 
AS WELL AS THE COMPARISON VARIABLE 


SUBROUTINE RESETI(II-12,13,%, 47 e) 
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Ree aeiizelene V( Tl, 12,13), 2@i1,12,13),C(3,2) 
DO 30 J1=1, il 

DO 20 J2=1,12 

efey 6, 33=1 13 

Awe eS iO © 
Y(J1 32 ,93)=0-0 
Z(J1,J3J2,J33)=0.0 
10 CONTINUE 

20 CONTINUE 

30 CONTINUE 


A=10.0E10 
DO 50 I=1,2 
DO 40 J=1,3 
(J,1)= 
40 CONTINUE 
A=-10.0E10 
50 CONTINUE 
RETURN 
END 
co I IE 
C**X SUBROUTINE RESET2 
CxxK 
C+++ 


Siete tis SUBROUTINE RESETS 3 SINGLE DIMENSIONED ARRAYS TO ZERO 


SUBROUTINE RESET2(11,X1,12,X2,13,X3) 
REAL X1(I1),X2(I2), X3(13) 
DO 10 cee el 


X1(J)=0.0 
10 CONTINUE 
DO 20 J=1,12 
X2(J)=0.0 
20 CONTINUE 
DO 30 J=1,1I3 
X3(J)=0.0 
30 CONTINUE 
RETURN 
END 
Cxxx 
C*x*x SUBROUTINE ROTATE 
Cx 
Ct++t 


C+++ THIS SUBROUTINE COMPUTES THE AXES ROTATION MATRIX 
+ 


SUBROUTINE ROTATE (THETAZ, PHIZ,THETAX, PHIX,R) 
REeC R( 3,3 

COMMON /PIDATA/ Pak nD 

THZ=THETAZ/RTD 

PHZ=PHIZ/RTD 

THX=THETAX/RTD 

PHX=PHIX/RTID 

R(1,1)=COS ( PHX aes 


R(2.1)=SIN(PHX) *SIN( THX 
R(3,1)=CcOS (THX 
R(1,3)=COS(PHZ ASIN(THZ) 
R(2.3)=SIN PHZ kSIN(THZ 
R(3.3)=Cos 
TERA a FRCL, “ADR (2 3)AR(2, L)+R(3, 3)*R(3,1) 
IF (ABS(D1) . ee 
C R(1,1)=R 21}: Si Seni 
C R(2.1)=R 21} “A a 3)*D1 
C R(3.1)=R(3,1 “A 3,3)*D1 
C D2=SORT(R(1,1)*R ie eaeL 1)*R(2,1)+R(3,1)*R(3,1)) 
C R(1,1)=R(1,1)/D2 
C R(2.1)=R(2.1)/D2 
C R(3,1)=R a! ae 
R(1,2)=R(2,3)*R(3,1)-R(3,3)*R(2Z,1) 


a5 


RCS SNOREL ARG 2} -BOB RG 
RETURN 











10 PRINT *, 'AXES SPECIFIED ARE NOT ORTHOGONAL! 
WRITE (2,%) "SUBROUTINE ROTATE ' 
WRITE (2,* AXES SPECIFIED ARE NOT ORTHOGONAL ' 
WRITE (2,*) ! THE INPUT PARAMETERS ARE AS FOLLOWS' 
WRITE (2,*) ! THETA-Z =',THETAZ,' DEG! 
WRITE (2,%*) ! PHI-Z =',PHIZ,' DEG! 
WRITE (2,*) | THETA-X =',THETAX,' DEG! 
WRITE (2,*) ! PHI-X  =',PHIX,!' DEG! 
CALL DONEPL 
STOP 
END 
CxAKK 
C**kx SUBROUTINE ROTRAN 
CxARK 
C+++ 
C+++ THIS SUBROUTINE ROTATES AND TRANSLATE CO-ORDINATES 
C++ 
SUBROUTINE ROTRAN(N,X,Y,Z,XC,YC,ZC,R) 
REAL K(N), Y(N) ,Z(N) ,R(3,35 
DO 10 I=1,N 
X1=X 1 
YI=V(I 
Z1=Z(1I 
X(I)=R(1,1)*X1+R(1,2)*Y1+R(1,3)*Z1+XC 
Y(I)=R(2,1)*K1+R(2,2)*¥Y1+R(2,3)*Z1+YC 
Z(I)=R(3,1)*X1+R(3,2)¥*V14+R(3,3)*Z1+ZC 
10 CONTINUE 
RETURN 
END 
C on an as ap an an on an an lan anon an lanl anton anlanlenlenlanl ae lanl anianlaslantanl asian lanl an ania! an! an lam) aul auam| amt am! am aut aman laa? cnt inl cash can] cnn) cont ceria fcr eee cer tere eee ee cee eet cee ee 
CxARK 
Cx*x* SUBROUTINE SCALES 
CARR 
C++t 


Ct+++ THIS SUBROUTINE SCALES THE GEOMETRIES 
Ctt++ ALPHABET “Sv ORs sOURGHS 
Cotta AND 'G' FOR ALL OTHER GEOMETRIES 


C+++ 
SUBROUTINE SCALES(N,AS,X,Y,2Z) 
REAL X(N), eG Z(N) 
CHARACTER A 
COMMON / SCALER/ FREQ, 1 LUN; 1US; SCALEE 
IF oe ae 
TUN GSE Oe SCALE= SCALEF 
iF IUN .EQ. 2) SCALE=0.3048*SCALEF 
IF (IUN .EQ. 3) SCALE=0.0254*SCALEF 
END IF 
IF (AS .EQ. 'S') THEN 
IF (IUS .EO. 0) SCALE=0. 2991925 / FREQ SoDue 
IF (1US .EQO. 1) SCALE=SCA 
TP OCLUS BE Oim 22 SGALE—O, SO48*SCALEF 
LP (1US 2 bOn SCALE=0 .0254*SCALEF 
END IF 
DO 10 I=1,N 
L =SCALE*X ZI 
Y¥(I)=SCALE*Y(I 
Z(1I)=SCALE*Z(I 
10 CONTINUE 
RETURN 
END 
CxAK* 
Cx*xx SUBROUTINE SORCEA 
CxxX 
C+++ 
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Cos eee 
at 
Cao 
C+++ 
ret 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 


ent 


C+++ 


THIS SUBROUTINE READS IN THE DATA FOR SOURCE ARRAY AND COMPUTES 
THE CORNERS OF EACH ELEMENTS. 
IT ALSO COMPARE THE MIN. & MAX. X, Y & Z VALUES 


CENTER Ob eaeem eLeMENT = XCSA,YCSA,ZCSA 
ANGLES NEEDED TO DEFINE THE SOURCE ARRAY CO-ORDINATE 
THETA TO™2ZTl AXIS : THSAZ 
PHI TO ZT AXIS : PHSAZ 
THETA TO AT AXIS : THSAX 
Polo m aan los) PHOAK 
AXES ROTATION MATRIX FOR SOURCE ARRAY CO-ORDINATE : RSA 


SUBROUTINE SORCEA(XSAC, YSAC,ZSAC) 
PARAMETER (NSX=1000) 

PARAMETER (NSAX=5) 

PARAMETER (NESAX=500) 

REAL BE ESA oC Oe ua ce neoan 
REAL XCSA NESAX) , YCSA(NESAX) , ZCSA(NESAX) 
REAL Oe s yn 

REAL COMPSA(3,2),R(3,3),RSA(3,3) ,COMPSG(3,2) 
INTEGER MSAX(NSAX) 

COMMON /RTIDATA/ XCR,YCR,ZCR,R 

COMMON /SCALER/ FREQ,IUN,IUS,SCALEF 

COMMON /SADATA/ ISA,MSAX,COMPSA,MSAXT 
COMMON /SGDATA/ ISG,COMPSG 





4) ,ZSAC(NSAX , NESAX, 4) 


WRITE ca) : 
WRITE (2,.*) ‘SUBROUTINE SORCEA' 
C+++ READ IN DATA 
WRITE (2,%*) ! ISA =',ISA 
READ(1,*) MSAX(ISA) 
WRITE (2,%*) | MSAX(ISA) =! ,MSAX(ISA) 


IF (MSAX(ISA) .GT. NESAX) GOTO 40 
MSAXT=MSAXT+MSAX(ISA) 
IST=ISA+tISG+tMSAXT 


WRITE (2,*) ! ISG =',ISG 
WRETE (2,*) | MSAXT =' ,MSART 
Wrime (2,~) ' IMS =" IST 


Helo ol. NSX) GOTO 50 
DO -10 Pe cee ay a 
READ(1,*) XCSA(I),YCSA(I),ZCSA(T) 


10 CONTINUE 


C+++ 
C+++ 
C+++ 


te 
C+++ 
C+++ 
C+++ 


CALL SCALES (MSAX(ISA),'G',XCSA,YCSA,ZCSA) 
READ(1,*) THSAZ,PHSAZ,THSAX, PHSAX 

WRITE (2,%*) THSAZ=', THSAZ,' PHSAZ=' , PHSAZ 
WRITE (2.%*) | THSAX=' | THSAX,' PHSAX=' | PHSAX 
READ(1,*) DUMMY,HS,HAWS 

WRITE (2.*) | HS =' HS,' HAWS =' HAWS 
READ(1,*} DUMMY, DUMMY 


COMPUTE THE AXES ROTATION MATRIX FOR SOURCE ARRAY 


CALL ROTATE (THSAZ,PHSAZ,THSAX, PHSAX, RSA) 


DETERMINE THE CORNERS OF EACH SOURCE ELEMENTS IN SOURCE ARRAY 
CO-ORDINATE 


DO 30 I=1,MSAX(ISA) 
=HAWS/2.0 
=HAWS /2.0 
=-1.0*HAWS/2.0 
=-1.0*HAWS/2.0 


ra 


“OOOO 


NF BWR PONE 
On 


Hu ud ue a u 
1roooo 


NN MRR HX OX 
mie « 


0 
fis 2.0 


o> 


Z ase Q*HS/2. 0 
4)=HS/2.0 
CALL SCALES(4,'S',X,Y,Z) 


C+++ 
EE: TRANSFER THE CORNER LOCATIONS TO GLOBAL CO-ORDINATE 
+++ 
CALL Sonnet ae op eons eae ZCSA(I),RSA) 
CALL ROTRAN(4,X, ae Z,XCR,YCR,ZCR, RD 
DO 20 J=1,4 
ASAC(ISA, I ,J)HK(I 
YSAC(ISA,I, we 
ZOAC(ISA,1I,J)=2(J 
20 CONTINUE 
C+++ 


C+++ COMPARE THE MIN. & MAX. X, Y & Z VALUES 


erat 
CALL COMPAR(4,COMPSA,X,Y,Z) 
30 CONTINUE 
Were 2 ‘SUBROUTINE SORCEA END' 


ETU 
40 CONTINUE 
PRINT *, 'THE NO. OF SOURCE ELEMENTS ="! MSA isa) 
PRINT *, 'NO. OF SOURCE ELEMENTS PER ARRAY EXCEEDS NESAX' 
Snes DONEPL 
SG 


50 CONTINUE : 
PRINT *, 'THE TOTAL NO. OF SOURGES)=—3 ten 
PRINT *, 'TOTAL NO. OF SOURCES EXCEEDS NSX'! 
eaee DONEPL 


TOP 
END 
Cree 
Cx*xx SUBROUTINE SORCEG 
CxrK 
Greet 


C+++ THIS SUBROUTINE READS IN THE SOURCE DATA AND COMPUTES THE CORNERS'! 
C+++ LOCATIONS 
C+++ IT ALSO COMPARES THE MIN. & MAX. XK, Y & Z VALUES 


C+++ 
SUBROUTINE SORCEG(XSC, YSC,ZSC) 
PARAMETER (NSX=1000) 
REAL XSC(NSX,4) ¥SC(NSX, 4), ZSC(NSK, 4),X(4),¥(4) ,2(4) 
REAL COMPSG(3,2),R(3 3) 
COMMON /RTDATA/ XCR, von men R 
COMMON /SCALER/ FREQ, IUN,IUS,SCALEF 
COMMON /SGDATA/ISG, COMPSG 
WRITE (2, 
WRITE Ss ‘SUBROUTINE SORCEG'! 
C+++ 
C+++ READ IN DATA 
C+++ 
WRITE (2,*) ! ToCu— ehce 
READ (1, *} RCS , ¥CS ,ZCS 
WRITE (2,%) XCS ='_xCS,' YCS =',YCS,' ZCS =',ZCS 
CALL scheaecl 'G' xcs, oe Es) 
READ(1,*) THSZ, PHSZ, THSK, PHSX 
WRITE (2,*) ' THSZ=' | THSZ,' PHSZ=' , PHSZ 
WRITE (2 oe THet=! THs | PHSX=' | PHSX 
READ(1,*} DUMMY,HS ,HAWS 
WRITE (2,*) | HS =',HS,' HAWS =!,HAWS 
; READ(1,*) DUMMY, DUMMY 
+++ 
C++ COMPUTE THE AXES ROTATION MATRIX FOR SOURCE CO-ORDINATE 
bat 
: CALL ROTATE (THSZ,PHSZ,THSX, PHSX,RS) 
ot 


on COMPUTE THE CORNERS Of THE SOURGE 
b+ 
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X(1)=HAWS/2.0 
X(2)=HAWS/2.0 
X(3)=-1. O*XHAWS/2.0 
X(4)=-1.0*HAWS/2.0 
Y(1)=0.0 
Y¥(2)=0.0 
Y¥(3)=0.0 
Y(4)=0.0 
2(1)=HS/2.0 
Z2(2)=-1. O*HS/2.0 

g Se ee pees 2 -0 


4)=HS/2.0 
CALL SCALES (4, US Ge a) 


C+++ 
C+++ TRANSFER THE CORNER LOCATIONS TO GLOBAL CO-ORDINATE 
p++ 
CALL aan): Kean apmcsuMecmacs, RS) 
- CALL ROTRAN(4,X,Y,Z,XCR, YCR,ZCR, R) 
+++ 
Ct++ COMPARE THE MIN. & MAX. X, ¥ & Z VALUES 
p++ 
DO 10 I=1,4 
XSC(ISG,I)=x(I 
YSC(ISG,I)=¥(T 
ZSC(ISG,1I)=2(I 


10 CONTINUE 
CALL COMPAR(4,COMPSG,X,Y,Z) 
WRITE (2,%*) ‘SUBROUTINE SORCEG END! 


RETURN 
END 
ae 
Cx** SUBROUTINE VANGLE 
CrxK* 
C+++ 
At THIS SUBROUTINE DETERMINES THE VIEWING ANGLES 
+++ 
SUBROUTINE VANGLE (TITLE, LENGTH) 
REAL LENGTH(3) 
CHARACTER ANSS 
CHARACTER*60 TITLE 
COMMON /VUDATA/ PHI ,THETA, RADIUS 
WRITE ie x 
WRITE (2,*) ‘SUBROUTINE VANGLE! 
C+++ 


C+++ COMPUTE THE MINIMUM DISTANCE RADMIN AWAY FROM THE CENTER OF 
C+++ THE WORK BOX 
C+++ ALSO ENSURE THAT THE RADIUS CHOSEN IS GREATER THAN RADMIN 


C+++ 

RADMIN=1. 0. 
DO 10 I=1, 

R He (ueeTH(2)4s, .0 .GT. RADMIN) RADMIN=REAL(INT(LENGTH(I)*5.0)) 
IF (RADIUS .LT. RADMIN) RADIUS=RADMIN 

20 CONTINUE 
CALL FRTCHS ( * CLRSCRN ') 
PRINT *, 
PRINT ; 
PRINT *, 'DO YOU WANT TO CHANGE THE VIEWING ANGLE?! 
PRINT *, 'THE PRESENT SETTINGS (FROM CENTER OF WORK BOX) ARE:! 
PRINT *, ' PHI IN DEG’ =',PHI 
PRINT *, ' THETA IN DEG =',THETA 
PRINT *, ' RADIUS IN M =!,RADIUS 
PRINT *, ANSWER Y OR N PLEASE! 
READ(5, (AL)") (2 NS 
IF os . 'Y! (AND. ANSS$ .NE. 'N') GOTO 20 


IF (ANSS £0. 'N') GOTO 40 
PRINT *, TER ao IN DEG! 
READ(S5,*) PHI 
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PRINT *. 'ENTER THETA IN DEG! 
READ(5,*) THETA 
30 PRINT *, 'ENTER RADIUS IN M! 
READ(5,*) RADIUS 
IF (RADIUS .GE. RADMIN) GOTO 40 
PRINT *, 'RADIUS MUST BE GREATER OR EQUAL TO ',RADMIN,! M' 
GOTO 30 
40 CONTINUE 
CALL FRTCMS('CLRSCRN_') 


WRITE (2,% PHL = rae 

WRITE (2,0 )) | THETA =!,THETA 
WRETE (2,7) | RADIUS =' ,RADIUS 
WRITE (2,*) 'SUBROUTINE VANGLE END' 
RETURN 

END 
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4, SA CARD GENERATING CODE (GTDSA) 
a. Flow Chart 


start 


reeos tn 





selects 
tapering type 


- YY 
C-fleld eap 
teper ing retio 


Jt 0 

oer es 
E-¢feld 

pecnit for 
aC 


2 | ome 


writes output 
Cats file 
according to 





C+++ 
Cr+ 
Crt 
C+++ 
C+++ 
C+++ 
C+++ 
Cr+ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 
C+++ 


b. Program Listing 


KAKKAKKKAAKKARKKRRAAKKKA AK 
* GTDSA * 


* BY KOH WEE JIN * 
‘ ON OCT 10 1986 i 


AT NPS 
KAKKKKKKKKKKKAKKKKKKKKK 


THIS PROGRAM 15 TO CREATE THE SA GARD NEEDED BY NEG=-BS@ Cope 
FOR PLANAR ARRAY 


Ligis Soba of producing the following Taperings for both the 
Xt and Zt directions : 


Uniform 

Cosine 
Cesithessquare 
Triangular and 
Binomial 


with either TE 10 mode or Uniform distribution for all elements 
The code uses the following assumptions: 


Simple rectangular array grid 

All elements are confined to Xt-Zt plane 
All elements have the same length, width and orientation 
E field is orientated along Xt direction 


The variabel names used and their meanings are listed below: 
Maximum Dimensions 


NESAX 


: Max. No. of elements per array 
NELSAX 


: Max. No. of elements jpommeuce 00 
Note that NELSAx as tieeded@towtimit the 
dimensions of the array variables 


Units and Wavelength 


IUN : Integer used to define geometry unit 
IUS : Integer used to define source unit 
UNITS : Dimension of unit specified 

UNITS +: Unit cConversieneratio sco m 

FREQ): acl ilar in GHz 

LAMBDA : Wavelength in m 


Rotation Angles and Matrix 


THSZ : Theta angle of Zt axis w. r. t. global eo=ord, 
PHSZ : Phi angle of Zt axis w. Ff. €. glebatecorance 
THSX : Theta angle of sxtjaxisow. ©. tigiobal co-encr 
PHSX : Phi angle of @iteacicm et. t- Gloso co loua. 
ROT : Axes rotation matrix 

LORTHO : Logic to check 1f Kt and Zt are orthogonal 


Array Area 


XALENG : Width of array area including edge spacing 
ZALENG : Length of array area including edge spacing 
XCAREA : X co-ord. of the center of array area 
YCAREA : Y co-ord. of the center of array area 
ZCAREA : Z co-ord. of the center of array area 
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e+ 
eit + 
e+ 
crrt 
err+ 
Cttt 
Gttt 
Ca++ 
et 
Chu 
Catt 
eEtt+ 
ett 
Cana 
er++ 
C+++ 
C+++ 
oa 
C+++ 
C+++ 
ea 
er++ 
ett 
ert 
or++ 
Crr+ 
ert 
ots 
G++ 
ent 
ort 
ee. + 
t+ 
++ 
octt+ 
et 
ot 
C+++ 
Osea 
Girt 
ot 
Csparae 
ot 
Cr+ 
ett 
cart 
c++ 
Cr ++ 
cert 


C++4+ 


Element parameters 


NX : No. of elements along Xt direction 

NZ : No. of elements along Zt direction 

A : Length of element in Zt direction 

B : Width of element in Xt direction 

XCE : X co-ord. of the center of each element 

VG : Y co-ord. of the center of each element 

ZCE : Z co-ord. of the center of each element 

XELENM : Element width in Xt direction in m 

ZELENM : Element length in Zt direction inm 

XELENG : Element width in Xt direction in the same unit 


specat tee for geometry (LUN) . 
ZELENG : Element length in Zt direction in the same unit 
specified for geometry (IUN) 


XESPAC : Element spacing in Xt direction 
ZESPAC : Element spacing in Zt direction 
Tapering 
IXTAP : Integer used to define the type of element 
pape rug Beene Re cebbec len 
ieee Uthcegem Usedmco celine tne type of element 
tapering along Zt direction 
AMPPED : Amplitude ratio of pedestal |. . 
XAMP : Amplitude tapering along Xt direction 
ZAMP : Amplitude tapering along Zt direction 
Power and E Field 
PTOTAL : Total input power in watts 
E : E field for each element in V/m 
ERATIO : E field amplitude ratio for each element 
BeOoUMe-moun Of themsquanesor np tield ratio 
EMAX ; Max. E field in V/m 
Constants 


: Speed of light divided by 1E9 |. 
DTR : Be aeee to radian conversion ratio 
ZO INtrinsic wave impedance in ohm 


PARAMETER Peaeageso). 
PARAMETER (NELSAX=270) 

REAL XCE(NELSAX,NELSAX) , YCE (NELSAX , NELSAX) , ZCE (NELSAX ,NELSAX) 
REAL XAMP(NELSAX) , ZAMP(NELSAX) 

REAL E(NELSAX,NELSAX) , LAMBDA, ROT(3,3) 

CHARACTER*1 ANSS 

LOGICAL LORTHO 

COMMON /EANGLE/ THSZ,PHSZ,THSX,PHSX 

COMMON /GENDAT/ C,PI,DTR,ZO 

COMMON /NECDAT/ RAM GROEN GGnenE oe 

COMMON /RTDATA/ ROT, XCAREA, YCAREA, ZCAREA 

COMMON /XZDATA/ NX,NZ,XALENG,ZALENG,A,B,IDIST 

C-07299/925 

PI=3.1415926 

DTR=P1I/180.0 

Z0=376.7 


C+++ Get NEC-BSC commands on units and frequency 


asa ae 


CALL NECCMD 
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Compute the center locations of each elements 
CALL ELMCEN(XCE,YCE,ZCE) 

Select the tapering types and pedestal height 
CALL TAPSEL(IXTAP,IZTAP,AMPPED) 

Compute the amplitude taperings 


CALL TEES cs 
CALL TAPER(IZTAP,'Z' ,ZAMP,ZCE 


Computer the E-field magnitude 
CALL EFIELD(XAMP,ZAMP,E,AMPPED) 
Write out the SA card data 


CALL WRITED(XCE,YCE,ZCE,E) 
STOP 


KAKRKKKAKKKKKKKKKKKAKAKKKK 


* SUBROUTINE EFIELD * 
KAKAKKAAKKAAKKAKKKAKKKAKA 


This subroutine computes the — field ine meron Cdeneclcmene 


SUBROUTINE EFIELD(XAMP, ZAMP,E,AMPPED) 

PARAMETER (NESAXE 45000) 

PARAMETER (NELSAX=270 

REAL XAMP(NELSAX) , ZAMP(NELSAX) , LAMBDA 

REAL ERATIO(NELSAX, NELSAX) , E(NELSAX, NELSAX) 

REAL UNITS(0:3 

CHARACTER*1 ANSS 

COMMON /GENDAT/ C,PI,DTR,ZO 

COMMON /NECDAT/ IUN,IUS,FREO, LAMBDA 

COMMON /XZDATA/ NX,NZ, XALEN /ZRLENG ABe] IST 

DATA UNITS(1),UNITS(2),UNITS(3)/1.0 3048, 0.0254/ 

UNITS(0)= =LAMBDA 

WRITE (2,% 

WRITE (2.% 

WRITE (2.% 

DO 10 I=1,NX 
WRITE (2,*) I , XAMP (1) 


‘Subroutine EFIELD' 
; . Xamp ! 


10 CONTINUE 


WRELE 27 iT 
DO’ 20° I= 7 NZ 
WRITE Zs I ,ZAMP(TI) 


Zamp | 


20 CONTINUE 


Crt+ 
C+++ 
C+t+ 


Cia 
Oar rte 
C+t++ 


30 


Convert element lengths to m 


SEE eee iarel easy 
ZELENM=A*UNITS(IUS 


Compute the relative amplitude 


ESQSUM=0.0 
DO 40 Taiz 
DO 30 J=1,N 
ERATIO(!, = “KAMP (J) AZAMP (1) (1. .0-AMPPED )+AMPPED 
ESQSUM= =ESQSUM+ERATIO(I,J) 
CONTINU 


40 CONTINUE 
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WRITE (2,*) '  ESQSUM = ',ESQSUM 


C+++ 
Senne Enter total input power 
eo 
PRINT * | 
PRINT *, 'Enter total input power in Watt' 
PRINT * 
READ (5,*) PTOTAL 
PRINT * 
WRITES(2,™) ' PTOTAL = ',PTOTAL,!' W'! 
C+++ 
C+++ Enter element distribution type 
C+++ TE 10 of Uniform 
C+++ 
50 CONTINUE 
PRINT * — 
PRINT *, ‘Enter element distribution type’ 
PRINT *, 2 Uniform distribution! 
PRINT *, '1) TE 10 mode! 
PRINT * 
READ (5,*) IDIST 
IF ee .NE. O .AND. IDIST .NE. 1) GOTO 50 
: IF (IDIST .EQ. 0) GOTO 55 
tee 


eett Compute maximum E field 
+++ 
Be eke azo! CEReNty MEVENMESOSUM) 
GOTO 5 
55 CONTINUE 
ee oe eee, 20/ ZELENM/AELENM/ESQSUN) 
58 CONTINU 
WRITE (2,*) ! EMAX = ' EMAX,' V/m' 
C+t++ 
C+++ Check if Emax exceeds 2.9 MV/m 


C+++ 
IF (EMAX .GT. 1.0631E6) THEN 
60 CONTINUE 


PRINT * 
PRINT *, 'Emax = ',EMAX,!' V/m' 
PRINT *, 'Emax is greater than the max. power flux constrain! 


PRINT *, ‘of 300 KW/sq-cm or equivalent to 10.631 KV/cm' 
PRINT *, 'Do you want to continue? ' 


PRON 
READ (5,'(Al1)') ANSS 
IF Nee .NE. 'N' .AND. ANSS .NE. 'Y!') GOTO 60 
IF (ANSS .EQ. 'N') GOTO 930 
END IF 
C+++ 
C+++ Compute individual excitations 


C+++ 
DO 80 I=1,NZ 
DO 70 J=1,NX 
E(1I,J)=EMAX*ERATIO(I,J) 

70 CONTINUE 

80 CONTINUE | 
WRITE (2,*) ‘Subroutine EFIELD end! 
RETURN 


90 CONTINUE 
SOE 


TO 

END 
C 
Cc KRAEKKKKKKKKKKKKKKKKKAKKKK 
& * SUBROUTINE ELMCEN * 
Cc KAKKKAKKAKKKKKKKKKKKKKKK 


C+++ This subroutine reads in the geometry data of the array elements 
Ctt+t+ and computes the center locations of each element 
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C+++ 


SUBROUTINE ELMCEN(XCE,YCE,ZCE) 

PARAMETER ES eA eeocoMy 

PARAMETER (NELSAX=270 ) 

REAL Tete Shee ,NELSAX), Le ees NELSAX) , ZCE (NELSAX , NELSAX) 
REAL ROT(3,3),LAMBDA ,UNITS (0: =) 

CHARACTER*10 UNITS(0: :3) 

CHARACTER*1 ANSS 

LOGICAL LORTHO 

COMMON /EANGLE/ THSZ,PHSZ,THSX,PHSX 

COMMON /GENDAT/ C,PI,DTR 20 

COMMON /NECDAT/ IUN, IUS, ,EREQ, , LAMBDA 

COMMON /RTDATA/ ROT, /XCAREA CAREA , ZCAREA 

COMMON /XZDATA/NX, NZ, XALENG , ZALENG, A, OP oor 

DATA cheno ,UNITS(2), UNITS (3)/1. 0, 0. 3088, ,O% 0254/ 


DATA UNITS 'wavelength', 7 en 
UNITS (0)=LAMBDA 

WRITE (2, } , 

WRITE (2,%) 'Subroutine ELMCEN' 


C+++ Enter array element orientation and compute axes rotation matrix 


Ct+++ 


CALL FRTCMS('CLRSCRN ') 


10 CONTINUE 


PRINT *, ‘Enter Array Element orientation (in deg) that' 

PRINT i 'defines Zt and Xt directions as used in NEC-BSC code'! 
PRINT x 'i.e. THSZ, PHSZ, THSx, PHSx" 

PRINT x! 

READ (5,%*) THSZ,PHSZ,THSX,PHSX 

CALL ROTATE (LORTHO ) 

IF (LORTHO ) GOTO 20 


PRINT 
PRINT *, ' THSZ = ',THSZ,' deg! 
PRINT *, ! PHSZ =? FHsa,' deg: 
PRINT *, ' THSX = ',THSX,! deg! 
PRINT *, | PHS. = ! /PHSX, ' deg' 
PRINT *, 'Axes specified are not orthogonal ' 
PRINT *, 'Re-enter the angles' 
PRINT * 
GOTO 10 
20 CONTINUE 
WRITE (2,*) ' THSZ = ',THSZ,' deg! 
WRITE (2,*) '! PHS4 = ',PHSZ,' deg' 
WRITE (2, 2 THSA = ',THSX,' deg' 
Whe (2. ee Pao = ' PHox,. dec 
C+++ 
ae Enter center location of the Array Area (in Global co-ord) 
bee 
CALL ERTCMS ( ' CLRSCRN ‘) 
PRINT 
PRINT 'Enter center co~ -ordinate of Array Area (X,Y,Z)' 
ERaNe Bs 'The unit is in “/UNEDSGnuN) 
READ (5, a Se mae Bee ZCAREA 
WRITE AREA = ',XCAREA, UNITS(IUN 
WRITE 5" YCAREA = ' YCAREA , UNITS IUN 
ze WRITE (2,*) '! ZCAREA = ',ZCAREA,UNITS(IUN 
+++ 
C+++ Enter number of elements 
C+t++ 
CALL FRTCMS('CLRSCRN ') 
30 CONTINUE 


PRINT a 'Enter No. of elements along Xt axis' 
PRINT * 
READ (5,*) NX 
IF (NX GI. NELSAX) THEN 
PRINT 
PRINT - 'No. of elements along Xt axis = ',NX 
PRINT *, 'No. of elements along Xt axis exceeds NELSAX' 
PRINT *, ‘re-enter the No. of élements' 
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PRINT * 
GOTO 30 
END IF 
85) CONTINUE 
PRINT * 
PRINT *, 
READ (5, i 


‘Enter No. 


) NZ 
IF (NZ -GT. NELSAX) THEN 


PRINT 
PRINT . 
PRINT x 
PRINT x 
PRINT * 
GOTO 35 
END IF 
NTOTAL=NX*NZ 


'No. 


of elements along Zt axis' 


'No. of elements along Zt axis = ',2X 
of elements along Zt axis exceeds NELSAX'! 
're-enter the No. 


of elements' 


IF (NTOTAL, .GT. NESAX) THEN 


PRINT * 
PRINT *, 'Total No. of elements = ',NTOTAL 
PRINT *” 'Total No. of elements exceeds NESAX'! 
PRINT * 're-enter the No. of elements’! 
PRINT * 
GOTO 30 

END IF 

WRITE ie x} ' NX = ' NX 

WRITE (2,*) ! NZ = ' NZ 


C+++ 
C+++ Enter the size of element 
ett 
40 CONTINUE 
CALL FRTCMS('CLRSCRN ') 


PRINT * 

PRINT *, 'Enter the size 

PRINT * 

PRINT *, ‘the same size' 

PRINT * 

PRINT *, ‘Enter the element len 
PRINT *, 'The unit is in 

PRINT * 

READ a x) B 

PRIN 

PRINT *, ‘Enter the element len 
PRINT * 'The unit is in 

PRINT a 


READ (5,*) A 


of the element! 


, 'Note that all elements, are assumed to have' 


th along Xt axis direction' 
“UNITS (I (IUS) 


gtr along Zt axis direction! 
' UNITS(IUS) 


c IF (B GT. A) THEN 
c PRINT * 
C PRINT *, 'THE ELEMENT LENGTH = ',A,UNITS(IUS 
c PRINT *, 'THE ELEMENT WIDTH = ',B,UNITS(IUS 
C PRINT *, 'THE WIDTH MUST BE SHORTER THAN THE LENGTH! 
C PRINT *, 'RE-ENTER ELEMENT SIZE' 
c PRINT * 
C GOTO 40 
S END IF 
WRITE $3 3 A = Panes 
WRITE (2,*) ' B = ',B,UNITS(IUS 
et 


C+++ Check if TE 10 mode exist 
C+++ 
ZELENG= AXUNITS(IUS) 


IF (ZELENG .LT. LAMBDA/2.0) THEN 
PRINT * 
PRINT *, ‘The operating freq = ',FREQ,' GHz' 
PRINT *, 'The corresponding weeMicaceh = | LAMBDA , | m' 
PRINT *, 'The length of array element = ',ZELENG,' m' 
PRINT *, ‘For TE 10 mode to exist, the ieee of the element' 
PRINT *, ‘must be greater than or equal to half wave length' 
eeu ni 'Re-enter element size'! 
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GOTO 40 


END IF 
C+++ 
ou Convert the element lengths to the geometry units 
+++ 
XELENG= -segerts( ius) forts {2a 
ZELENG=A*UNITS(IUS )/UNITS(IUN 
WRITE (2,% XELENG = ',XELENG, Rieeraeiny 
: WRITE (2,7%0"' ZELENG = ',ZELENG, UNITS(IUN 
+++ 


C+++ Enter gap size between elements 


Cit 
CALL FRICMS('CLRSCRN ') | 
PRINE @,  Hiitewerne aga. size between elements injxt dtrece7cny 
PRINT *, 'The unit is in ',UNITS(IUS) 
PRINT * 
SOT (3, *) XGAP 
ERINE *, (ENGemer de gap size between elements in Zt direction! 
PRINT *, 'The unit is in ',UNITS(IUS) 
PRINT * 
READ (5,*) ZGAP 
OnE b= =xGh Ores US| te 
ZGAP1=ZGAP*UNITS(1IUS)/UNITS(IUN 
XESPAC=XGAP1+XELENG 
ZES PAC=ZGAP1+ZELENG 
XALENG=XESPAC*REAL(NX+1 ) 
IF (NX .EQ. 1) XALENG=XELENG 
ZALENG=ZESPAC REAL(NZ+1 ) 
IF (NZ .EQ. 1) ZALENG= ZELENG 
Z XGAP 


WRITE oe = ',&8GAP, NTT ee 
WRITE ,e) | ZGAP = 13 AP ,UNITS(IUS 
WRITE (2,4) XESPAC = ', XESPAC,UNITS(IUN 
WER EOE 02,3) am ZESPAC = ! ,ZESPAC, UNITS ( IUN 
WRITE (2,%) | XALENG’ = | /XALENG, cere IUN 
Weitere (2 ae ZALENG = | ZALENG , UNIT IUN 


Grr 
C+++ Determine the center locations of each element in Xt-Zt plane 
Cra 
CXXX WRITE (2,*) ' I J XCE!, 
CXXXK * TGs ZCE! 
DO 70 I=1,NZ 
Z=-1. oe .O+FI*ZESPAC 
IF (NZ .EQ Aa 
DO 60 J=1,NX 
X=-1. ag 4% .0+J*XESPAC 
IF (NAGE 1) X=0.0 
ACEC Ig 
1CE a - =0.0 
Z2CE =Z 
CXXXK WRITE Gai ) 1,d,XCE(i, J), YGE( 1, 0) Gennes 
60 CONTINUE 


70 CONTINUE | 
WRITE (2,*) ‘Subroutine ELMCEN end! 


RETURN 
END 
C aaa nmr ann awa ew ae eee ae eee ee eee eee ee ee eee eee eee ee ee Oe eee eee ennaaesaeannaanaaas & 
KAEAKKAAKKKKKKKKKKAKAKAAKK 
c * SUBROUTINE NECCMD * 
Cc KAARKARAKKKKKKAKKKAKKKKKK 
SUBROUTINE NECCMD 
REAL LAMBDA 
COMMON /GENDAT/ C,PI,DTR,ZO 
COMMON /NECDAT/ IUN, LUS FREQ, LAMBDA 
WRITE (2, im 
WRITE ae 'Subroutine NECCMD'! 
C+++ 


C+++ The following NEC-BSC commands are needed to determine 
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C+++ the units used by the whole data set 
C+++ They should be consistance with the rest of the data 


C+++ 
CALL FRTCMS ( ' CLRSCRN 4 
PRINT ‘Enter Bae fo lowing NEC-BSC commands' 
PRINT ae rune Us. 
PRINT *, ‘where UN mane US are eects. 
BRING =)! =1 unit in 
oGeit Ee UN = 2 eas, abyg ft! 
PRINT *, ! UN = 3 Une eed in: 
PRINT * Aad 
PRINT *, ! US = 0 unit in wavelength' 
PRINT *, ! US = 1 Uiieee lia, m 
BRONe *, | US = 2 unit in ft! 
Pein: *<, ' US = 3 Ware in in’ 
PRINT * | 
PRINT *, 'FR is frequency in GHz' 
PRINT * 
READ (5,%*) IUN,IUS,FREQ 
LAMBDA=C/FREO 
WRITE (2,*) ' IUN = ', TUN, 1US = * -1Us 
WRITE (2,*) ' FREQ = NEREO,’ GHz 
WRITE (2,*) ' LAMBDA = ',LAMBDA,' m' 
WRITE (2,*) 'Subroutine ween Bnd! 
RETURN 
END 
c 
C KAKKKKKKKRKKERKKKKKKRKKRK 
C * SUBROUTINE ROTATE * 
C KRARKKKKKAKKRKKKKRKKKKKA 
C+++ | 
Set This subroutine computes the axes rotation matrix 
S++ 
SUBROUTINE ROTATE (LORTHO) 
REAL ROT(3,3) 
LOGICAL LORTHO 
COMMON /EANGLE/ THSZ,PHSZ,THSX,PHSX 
COMMON /GENDAT/ C,PI,DTR, Z0 
COMMON /RTDATA/ ROT, XCAREA, YCAREA ,ZCAREA 
LORTHO=.TRUE. 
WRITE (2, H | 
WRITE (2,*) ‘Subroutine ROTATE! 
THZ=THSZ*DTR 
THX=THSX*DTR 
PHZ=PHSZ*DTR 
PHX=PHSK*DTR 
ROT (1,1 )=COS( PHX pera 
ROT(2,1)=SIN( PHX )*SIN( THX 
ROT (3 ,1)=COS (THX 
ROT (1,3 )=COS (PHZ eee 
ROT fe =SIN(PHZ)*SIN(THZ 
ROT (3 ,3)=COS ( THZ 
D1i= Soil, Me ce Gat eee or 1)+ROT(3,3)*ROT(3,1) 
IF (ABS(D1 Gi 0. GOT 
eee ROV(1,1)=ROT(1,1 -ROT 3 roe 
Gee? §6ROT 7a =ROT(2,1)-ROT(2,3 Boe 
Maes) ROT(3,1)=ROT(3,1)-ROT(3,3 
C??? D2=SORT(ROT(1,1) ROT (1,1 1 )+ ROr(2, 1)*ROT(2,1)+ROT(3,1)*ROT(3,1)) 
C??? ROT(1,1)=ROT(1,1)/D2 
C??? ROT(2,1)=ROT(2,1)/D2 
C??? ROT(3,1)=ROT(3,1)/D2 
ROT(1,2)=ROT(2,3)*ROT(3,1)-ROT(3,3)*ROT(2,1 
ROT(2,2)=ROT(1,1)*ROT(3,3)-ROT(1,3)*ROT(3,1 
ROT(3,2)=ROT(1,3)*ROT(2,1)-ROT(2,3)*ROT(1,1 
WRITE (2,*) ‘AXIS ROTATION MATRIX' 
Wkeebee(2,*) ROT(1,1),ROT(1,2),ROT(1,3 
Winner 2) ROT(2,1),ROT(2,2),ROT(2,3 
Wieeeec ~ ROL(S,1),RO1(3,2),ROT(3,3 





WRITE (2,%*) ‘Subroutine ROTATE end! 
RETURN 


10 CONTINUE 
LORTHO=.FALSE. 
WRITE (2,%) ‘Subroutine ROTATE end! 
Sa 
EN 


RRARKRRRRRRKRRRRRKRRRRRRRRE 


* SUBROUTINE ROTRAN * 
KARAKKKARKKARAKRRAKKKAKK 


C+++ This subroutine rotates and translate co-ordinates 


SUBROUTINE ROTRAN(X,Y,Z) 
PARAMETER (NELSAX=270) 
REAL X(NELSAX,NELSAX) , Y(NELSAX,NELSAX), Z(NELSAX, NELSAX) ,ROT(3,3) 
COMMON /RTDATA/ ROT, XCAREA, YCAREA, ZCAREA 
COMMON, / XZDATA/ NX,NZ, XALENG, ZALENG,A,B,IDIST 
N 


*Y1+ROT *XZ1+YCAREA 


HL tROt 
*Y1+ROT AZ1+ZCAREA 


1 
tee (22 
S72 


siLRon 2/3} 2LHYCAREA 
a, 
SAL tROT 23 





Zid 
10 CONTINUE 
20 CONTINUE 
RETURN 
END 


RREREAKARKRARRRRRRRKRARRRERE 


* SUBROUTINE TAPER ~ * 
KRAKRKAKAARKKAAKRKARKAAK 


C+++ This subroutine computes the distribution in Xt or Zt axis 


SUBROUTINE TAPER(ITAP,XZ$,AMP, CENTER) 
PARAMETER ce 
PARAMETER (NELSAX=270) 
REAL AMP(NELSAX) , CE (NELSAX) , CENTER(NELSAX, NELSAX) 
CHARACTER*1 XZ$ 
COMMON /GENDAT/ C,PI,DTR,ZO 
COMMON /XZDATA/ NX,NZ,XALENG, ZALENG,A,B,IDIST 
AMPMAX=0 . 0 
IF (XZ$ .EQ. 'X') THEN 

N=NX 


ALENG=XALENG 
DO 10 I=1,NX 
CE (1)= ‘CENTER (1,1) 
10 CONTINUE 
ELSE 
N=NZ 
ALENG=ZALENG 
DO 20 I=1,NZ 
CE(I)=CENTER(I,1) 
20 CONTINUE 
END IF 
IF (ITAP .EQ. 1 .OR. ITAP .EQ. 2) GOTO 40 
ITAP .EO. 3) GOTO 50 
ITAP .EO. 4) GOTO 60 


C+++ 


C+++ Uniform tapering 
C+++ 


DO 30 I=1,N 
AMP = ia 
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30 CONTINUE 
RETURN 


C+++ 
C+++ Cosine and cosine square tapering 
C+++ 
40 CONTINUE 
J=ITAP 
DO 45 I=1,N 
Beda dope ane ay) AEM 
IF (AMPMAX .LT. AMP(I)) AMPMAX=AMP(T) 
45 CONTINUE 
GOTO 70 
Ct+t++ | 
C+++ Triangular tapering 
C+++ 
50 CONTINUE 
DO 55 I=1 


,N 
AMP Der tent ma eT MY one 
IF (AMPMAX .LT. AMP(I AMPMAX=AMP (T ) 
55 CONTINUE 


GOTO 70 
CHt+ | 
C+++ Binomial tapering 
C+t+ 


60 CONTINUE 
DO 66 I=0,N-1 
COEFF=1.0 
pee Ca geen OR. I .EQ., N-1) GOTO 64 
DO 62 J=1,I 
COEFF=COEFF*(N-1l-J+1)/J 
62 CONTINUE 
64 CONTINUE 
12 lanes .LT. COEFF) AMPMAX=COEFF 
AMP I+1)=COEFF 
66 CONTINUE 
C+++ 
C+++ Normalization 
C+++ 
70 CONTINUE 
DO 80 I=1,N 
AMP ( )=AMP (I) /AMPMAX 
80 CONTINUE 
URN 


RET 

END 
c 
Cc RAKKKKKKKAKKKARKKKKKKRER 
c * SUBROUTINE TAPSEL * 
Cc KAKKKKKKKKKKKKKKKKKKKKA 


Sette Lhis subroutine selects the Xt and Zt axis tapering type 
C+++ and the tapering pedestal height 


SUBROUTINE TAPSEL(IXTAP,IZTAP,AMPPED) 
WRITE (2,* | 
WRITE (2,*) 'Subroutine TAPSEL' 
10 CONTINUE 
CALL FRTCMS('CLRSCRN ') | } | 
PRINT *, 'Select one of the tapering type in Xt direction' 


PRINT * | 

PRINT *, '0) Uniform! 

PRINT *, '1) Cosine! 

PRINT *, '2) Cosine square' 

PRINT *, '3) Triangular’ 

PRINT *, '4) Binomial' 

PRINT * | 

PRINT *, ‘Enter the number before the tapering type’ 


PRINT *- 
READ (5,'(I1)'!) IXTAP 
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TF (CIXTAP .LT. O° OR] IXDAP =. Giwee) GOnOmsO 
20 CONTINUE 


PRINT * 

PRINT * } | | 
ERENT er 'Select one of the tapering type in Zt direction' 
PRI 

PRINT *, '0) Uniform! 

PRINT *, '1) Cosine' 

PRINT *, '2) Cosine square' 

PRINT *, '3) Trrangquiam 

PRINT *, '4) Binomial! 

PRINT * . 

PRINT *, 'Enter the number before the tapering type' 
PRINT x! 


READ (5, SE ) IZTAP 

IF (IZtae -LigsGngente IZTAP.Gr. 4) Gonome) 

CALL FRTCMS ( 'CLRSCRN ) 

PRINT *, ‘Enter pedestal nedehe relative to max amplitude’ 
PRINT a 'Enter 0.0 if no pedestal is needed! 


PRINT 

READ (5,*) AMPPED 

WRITE. (2. : IXTAP = | TXTAP 
WRITE (2,*) ! IZTAP = ',IZTAP 
WRITE (2,%) ! AMPPED = '!  AMPPED 
WRITE (2,*) 'Subroutine TAPSEL end! 
RETURN 

END 


KREKKKREKKKKKRKRKRKRRKAKKKKE 


* SUBROUTINE WRITED * 
RRAEKKKKKKKKKKKKKAKKKKAKK 


C+++ This subroutine writes the output data in the format 
Ctt+ sé@t by NEC-bSelcede 


SUBROUTINE WRITED(XCE,YCE,ZCE E) 
PARAMETER (NESANS 43000) 
PARAMETER (NELSAX=270) 
REAL XCE(NELSAX,NELSAX), YCE (NELSAX, NELSAX) ,ZCE(NELSAX, NELSAX) 
REAL E(NELSAX,NELSAX) , ROT 
COMMON /EANGLE/ THSZ,PHSZ, mek, PHSX 
COMMON /RTDATA/ ROT, XCAREA, YCAREA, ZCAREA 
COMMON /XZDATA/ NX, NZ, XALENG, ZALENG,A,B, IDIST 
NTOTAL=NX4*NZ 
WRITE (1,1000) NTOTAL 
CALL ROTRAN(XCE, YCE,ZCE) 
DO 20 I=1,NZ 
DO 10 J=1, 
WRITE ep 1010) XCE(1,J),YCE(1 go) eZenar aD) 
10 CONTINUE 
20 CONTINUE 
WRITE (1,1020) THSZ,PHSZ,THSX,PHSX 
IF fiDi8t EQ. 0} IMS=1 


LES IDIST, 2-20. A) eiise=s 
WRITE (1,1030) IMS,A,B 
PHASE=0.0 
DO 40 I=1,N2 

DO 30 J=1 


WRITE ce 1040) E(1,J),PHASE 

30 CONTINUE 

40 CONTINUE 
1000 FORMAT (I6) 
1010 FORMAT(F12.7,1H,,F12-7,1H,,F9-4) 
1020 FORMAT(F6.1,3(1H, ,F6.1)) 
1030 FORMAT(12,1H,,F9.5,1H,,F9.5) 
1040 FORMAT(E12.6,1H, ,F6.1) 

RETURN 


END 
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APPENDIX B 
RADIATION PATTERN PLOTS 
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(c) Element Length = 1.0A 





Figure B.l Element Width = 0.25, Element Length Variable, 
Z, Spacing = 1.0A and X, Spacing = 0.5. 
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Figure B.2. Element Width = 0.25A, Element Length Variabie, 
X, Spacing = 0.5A, Z, Spacing = 0.75A for (a) & (b) 
and Z, Spacing = 0.5A for (c). 
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Figure B.3. Element Width = 0.75), Element Width Variable, 
X, Spacing = 0.75A, Z, Spacing = 0.75A. 
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Figure B.6 Vertical Cut at @ = 0° with Fixed Array Area of 
7.5. X 3.75. Uniform Tapering in both Directions. 


{ 


‘, 


| 
/ 





a0 46 
; : 
g | \ § | 
S 1 g - ‘| 
7) | | | i 
8 4 
3 = 
.a Ni, $ -” i 
Po a Mt it. ugh - Wa 
f 
_ é, Mi 
| | 
| ~ 
“110 een 110 
a a ma me re eee $0 10 160 200 300 330 400 6 $90 100 0 200 230 360 350 400 
Angle 6 tieg) Angle 8 Geg 
(a) Element Size = 0.75A x 0.375 tb) Element Size = 0.75 x 0.25A 
No. of Elements = 10 x 10 No. of Elements = 10 x 15 
-10 
-20 
-30 
04 | / 
2 i| om (| 2 
3 -40 3 
we we 
-80 
3 3 
3 -- Vig. sace ial _ z 
e | "\A | } tl . if im 
00 | , 2 
S \ 
=100- 
@10 
50 co 196 209 250 390 356 400 o 50 160 139 269 330 390 330 400 
Angie 6 Keg) Angle 6 Weg 
(¢) Element Size = 0.54 x 0.3754 (d) Element Size = 0.5A «x 0.254 


No. of Elements = 


15 x 10 


No. of Elements = 


15 x 15 


Figure B.7 Vertical Cut at @ = 90° with Fixed Array Area of 


7.54 x 3.75k. Uniform Tapering in both Directions. 
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Figure B.8 Vertical Cut at @ = 0° with Fixed Array Area of 
7.54 x 3.75k. Cosine Tapering in both Directions. 
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Figure B.9 Vertical Cut at @ = 90° with Fixed Array Area of 
7.5. x°3.75. Cosine Tapering in both Directions. 
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Figure B.10 Far Field and Near Field Plots with Uniform Tapering. 
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Far Field and Near Field Plots with Cosine Tapering. 
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Figure B.12 Far Field and Near Field Plots with Cosine Square Tapering. 
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Figure B.13 Far Field and Near Field Plots with Triangular Tapering. 
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Figure B.14 Far Field and Near Field Plots with Binomial Tapering. 
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Figure B.15 Uniformly Tapered Linear TE,, Array with 
20 Elements of Sizes 0.5A x 0.25 each. 
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Figure B.16 Single Uniformly Distributed Slot Antenna with Size 
10.00 x 0.25). 
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Figure B.17 Cosine Tapered Linear TE,, Array with 
20 Elements of Sizes 0.54 x 0.25A each. 
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Figure B.18 Single TE,,) Mode Slot Antenna with Size 10.54 x 0.25A each. 
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Figure B.19 Cosine Tapered Linear Array with 20 Uniformly Distributed 
Elements of Sizes 0.5 x 0.25X each. 
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Figure B.20 Far Field Patterns of Simplified Medel with 


Uniformly Distributed Elements and Triangular Tapering. 
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Figure B.21 Near Field Patterns of Simplified Model with 


Uniformly Distributed Elements and Triangular Tapering. r 
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(c) Horizontal Cut at @ = 90° 


Figure B.22) Near Fieid Radiation Patterns of Uniform - Uniform Tapered Array. 
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Figure B.23. Near Field Radiation Patterns of Cosine - Uniform Tapered Array. 
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Figure B.24. Near Field Radiation Patterns of Uniform - Cosine Tapered Array. 
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Figure B.25_ Far Field Radiation Patterns of Array using 
Uniform - Cosine Tapering with Reflector Plate Tilted at 45°. 
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Figure B.26 Far Field Radiation Patterns of Array using 


Uniform - Cosine Tapering with Reflector Plate Tilted at 22.5°. 
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Figure B.27 Far Field Radiation Patterns of the Full Svstem Model 
with Reflector Plate Tilted at 45°. 
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Figure B,28 Far Field Radiation Patterns of the pul ovcteni™ ocd) 


with Reflector Plate [Tilted at 22.5°. 


Preax (dBWAq re) 
& 


Magnitude 


6 a 100 


180 
Angle 


(a) Vertical Cut 





Poeax (dBWAgq mw) 
3 3 
ae 








3 
%, oO 
2 
-20 
30 : | t | 
g “40 t 
200 230 300 330 408 0 30 wo 150 200 250 300 230 400 
6 Weg) Angle 6 deg 
at @ = 0° (b) Vertical Cut at @ = 90° 
as : 


Perea (dBWAgq m) 


Magnitude 





s 10 30 200 236 300 330 406 


Angle g (deg) 
(c) Horizontal Cut at 6 = 990 


Figure B.29 Near Field Radiation Patterns of Full Modei with Refiector Piate 


Facing @ = 0° and Tilted at 45°. 
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Figure B.30 Near Field Radiation Patterns of Full Model with Reflector Piate 


Facing @ = 0° and Tilted at 22.5”. 
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Figure B.32 Near Field Radiation Patterns of Full Model with Reflector Plate 
Facing @ = 90° and Tilted at 45°. 
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APPENDIX C 


GAR AND MAXIMUM APERTURE 


DERIVATION OF 
S BETWEEN ACTUAL AND 


RELATIONSHIP 
The derivation of the gain and maximum aperture E-field relationships between 

the actual linear array using either Uniform or Cosine (TE,,) distribution elements and 
those simplified single element models of equivalent length are shown in this appendix. 
Based on the different configurations of the array elements, we subdivide them into 3 


sections, namely, 


¢ Uniformly tapering TE,, linear array vs single uniformly distributed slot 
antenna. 

° Cosine tapering TE), linear array vs single TE,, slot antenna. 

e Cosine tapering uniformly distributed linear array vs single TE,, slot 
antenna. 


For the case when the array is formed by uniformly tapered uniformly distributed 
element, the simplified model of single uniformly distributed aperture will yield exact 
solution and hence no derivation is needed. 

I. UNIFORMLY TAPERED TE10 LINEAR ARRAY VS SINGLE UNIFORMLY 

DISTRIBUTED SLOT ANTENNA 

a. Comparing Maximum Aperture E-Field 

For single Uniform Aperture, 


es 
E a,b 


Where E, and a,” are the maximum field strength and the length of the uniform 
aperture respectively. b is the width of each element and it is the same for both the 
actual array and the simplified model. Zp) is the intrinsic impedence. 


As for the 20 elements uniformly tapered linear area, we have 


Tr 2 
| ie ab 
n 
! 4Z, 


ley 


for E, and a, are the maximum field strength and the length of each element of the 
Cosine distributed (TE,,) aperture respectively. With the same input power and 


ay = Na, we have, 


E_ = 2E 


Cc u 


Hence, the simplified model gives only half the field strength of the actual 


array. 


b. Comparing Gain 
The boresight far field field strength of an uniformly distributed aperture can 


easily be obtained from the available literature [Refs. 5,6], 


Ea 


Be ee 
FE upl 7 rh 


where r is the far field distance and i is the wavelength. The boresight far field field 
strength of a TE,, mode rectangular aperture on the other hand, can be derived from 
the two-dimensional Fourier transform of the aperture field given by [Ref. 5] 
b/2 .a/2 A 

P=0, P= ©  E. exp(jBr.r’) dx,’ dz,’ 

y Je Le u ExPUIBE-r) dx,’ dz, 
where r is the unit vector in the direction of propagation and r’ is the vector pointing 
towards the aperture. The boresight far field field strength from one element is then 


found to be equal to 
2E.ab 
trA 


Assuming that the distance r is so far that all fields arriving at that point from 


all the N elements are in phase, we then have 


EI 2NE.a.b 
oe trh 
Hence, the gain ratio between the two models 1s 


al oe Nees 


LE ap! TE ae 
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Since E, = 1/2E, anda, = Na,, assuming both apertures are fed with the same input 
power 


Eu! 4 


ea T 

This means that the gain of the actual array is about 2.098 dB higher than that of the 
simplified model. 
2. COSINE TAPERED TEIO LINEAR ARRAY VS SINGLE TE10 SLOT 

ANTENNA 

a. Comparing Maximum Aperture E-Field 

The maximum aperture E-field of a single TE,, slot antenna with length a’, is 
given by, 
ALP 


E’ 2 ame —— SS ae 
Cc , : 
a (> 


while that of a Cosine tapered TE,, linear array is 


E — aL Pia 


2. 
ab YE 20) 
n= 


where E(n) is the relative E-field of the n“ element with maximum relative E-field 


9 


equal to unity. Since Cosine tapering is used, the relative E-field may be expressed in 


terms of n as follow: 





E(n) = sin 


3 


iNet 


and the total sum over N elements is 


¥ E20) = —— 
n=1 * 2 


Therefore, the maximum aperture E-field is expressed as, 





Dale 8 LoP in 


: a.b(N+1) _ 


The ratio of the maximum aperture E-field for the same input power 1s 


is? 


Ee co 
EF. vila (N+) 





By equating a. = (N+1)a,, we have 
NOE 


The actual array maximum aperture E-field is equal to that of the simplified model 
multiplied by ./2. 


b. Comparing Gain 
As derived before, the boresight far field field strength of the single TE,, 
aperture antenna with length a.’ is given by 
2 one 


Ea 
IE ob mrr 


For a Cosine tapered TE,, linear array with far field distance r long enough 
to guarantee in-phase E-field addition, we have 


aes 
Eel = ag lo 


ass sin{t/(N + 1)] 
Trh l - cos[t/(N + 1)] 
2a bE : 
= —&—& cot{n/[20N+1)]} . 
trh 


Hence, the gain ratio is, 


E aE 
ou = ——< cot{m/[2(0N+1)]} . 
ie cb! a oe 


From above, we have E, = RI DEX. Esa fOf deen cea n) ae 





E Z 
JE pl = v cot{m/[2(N+1)]} . 
Eo! pial 


Table 4 shows the gain ratio as N varies. 


From the table, it shows that for N > 8, the boresight gain difference 1s 
reduced to | dB or less. 
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TABLE 4 


GAIN RATIO BETWEEN COSINE TAPERED TE10 LINEAR ARRAY AND SINGLE TE10 
SLOT ANTENNA 


po | ati |i ference in aB 
0.8165 “1.7610 


a 
[To SidSSSCSC~—tsg SiC 
[oi Stes |S 
[so S* SSS 000 Sig 
0 [0.8123 





COSINE TAPERED UNIFORMLY DISTRIBUTED LINEAR ARRAY VS 
SINGLE TE10 SLOT ANTENNA 
a 


. Comparing Maximum Aperture E-Field 
The maximum aperture E-field of a Cosine tapered uniformly distributed linear 


array with element length a,’ can be expressed as follows: 


ZoP.. 


a,b > E7(n) 


2Z, P.. 
ancGNis 1) 


and the maximum aperture E-field of a single TE,, slot antenna is 


E = 


U 


E2= 4ZoP 
. a’.b 


Hence, the maximum aperture E-field ratio is, 


Ee Ja’. 
E J {2a (N+ 1) 


Cc 





14] 


By letting a’, = a, (N+ 1), 





Eon 
Bs af 2 


Cc 


b. Comparing Gain 
The boresight far field field strength of a single element uniformly distributed 


aperture is, 


u_uUu 


ri 


as given before. With N Cosine tapered elements, the far field field strength is 


E-aeb 


E ab 
E =|- UU 
| ue rr nat 
E ab 
= a cot (R/[2IN+ 1}, 
r 


E(n) 


assuming r is large enough to guarantee in-phase E-field addition. 
The boresight far field field strength of the single TE), aperture is, 


2E’ ab 
IE el ~ ae 
Ttrh 


Hence, the gain ratio is, 


E REA 
uel = ete cot (aN +1] 
|E cb 2E ct | 


With a’, = a (N+1) and E, = (1//2)E’,, we have 


IEW 50 pee 
Ee 2f/2(N+1) 


Table 5 shows the boresight gain ratio as N varies. It shows that for large N, 


cot{m/[20N+1)]}} . 





the difference in dB approaches -3.01 dB. 







TABLE 5 


GAIN RATIO BETWEEN COSINE TAPERED UNIFORMLY DISTRIBUTED LINEAR 
ARRAY AND SINGLE TE10 SLOT ANTENNA 











og | mb | citmatne 
es 
ee 
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APPENDIX D 


INPUT DATA FILE FOR GTD2P51 


CT (FWRP1W1F)S 
CE 


UN 
US 
FR 
Ri 


PG 


EG 


PG 


PG 


PG 


PG 


Fe 


CG 


CG 


CG 


GEQUETRY UNIT IN uM 


SOURCE UNIT IN WAVELENGTH 


5 FREQUENCY IS 10.0 GHZ 


ROTATE THE ANTENNA REFLECTOR 


0,08 

45 ,0, 30270 

nae REFLECTOR PLATE 
1,35;,0 (05m=2500 
=1.35, 0 USneaeoc 

=1 30 Umno ae. 


1.357 01052200 
REDLECTOR ELUATE 


4,0 

LoS oa O07 sees 

le 3e7 50 Oo eos 

=i Oe 

-1.35,- ° ° 
REFLECTOR]: Ua Er 


4,0 
Wel So. OOO 2 OO 


“On 90, 270 


FRONT PLATE 


BACK PLATE 


TOE VE EATE 


BOLTON PLATE 


2) ee 2 FV 


RIGHT FUATE 


Dae 
ROTATE’ THE ARRAY ANTENNA 


PIVOT FOR REFLECTOR 


-1545 0,0 
SACL, 90, 90 
O. 14, 0.14 
-0.11,90,0.11,90 


PIVOT FOR REFLECTOR 


1.46,0,0 

90,0 90, 30 

Oma 0.14 

“0. il, 30,0. LL, 30 


SIDE SUPPORTING CYLINDER 1 
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CG 


CG 


CG 


CG 


CG 


PG 


PG 


PG 


PG 


PG 


SA 


moa 0 ,-1.37 
Poeo.e 0, 108. 314,0 
OOS, 0). 05 
amie 1 35.660,1020, 73.686 
Baie SUPPORTING CYLINDER 2 
ioe ,0,<1.37 
16.314, 180° FOGeore, 150 
Oe 05,0. Ue 
geez?) /3.686,1.28, 73.686 
SIDE SUPPORTING CYLINDER 3 
meee ,0.595,-1.35 
e270 ,112,-90 
gees, 0.05 
eee, 65,1.35,100 
SIDE SUPPORTING CYLINDER 4 
gomeeo, -0.595,-1.35 
ae), 112,90 
05 


ao, 06,1.35,100 
SIDE SUPPORTING CYLINDER 5 
ero, 0.595,°1.35 
gem 20,112,-90 
woo, 0.05 
goet,086,1.35,100 
SIDE SUPPORTING CYLINDER 6 
mies, -0.595,-1.35 
227 0),112,90 
oFO5,0.05 
aot, 06,1.35,100 
ARRAY ANTENNA GROUND PLATE 


Sle 

ARRAY ANTENNA SIDE PLATE 1 
Bi 

ARRAY ANTENNA SiDE PLATEeZ 


5. -3.1 
fa ARRAY ANTENNA SIDE PLATE 3 
52;25,1.5,52,6 


Pee) .5,-2. 
ee. 5, 73.1 
mec eD ,=3.1 
ARRAY ANTENNA SIDE PLATE 4 


4,0 
eareorn 1. 5.72.6 
ae, 1.5,-2.6 
=2.25, 1. car = Sic 
22.25, 5, ace 
ARRAY ELEMENTS CONFIGURATION 


S53 

mOw00GS, -0.9893153, <-2.5999 
-0000003, -0.9743257, -2.5999 
-0000003, -0.9593360, -2.5999 
.0000003, -0.9443464, -2.5999 
meewenvos, =-0.9293568, <-2.5999 
moeweoos, —-0.9143671, -2.5999 
muUNV003, -0.8993775, 2.5999 
-0000003, -0.8843879, <-2.5999 


ee rs rs rs rs re ee eG ee es | 
DOOD DCD DOOD DODD OCOD OC OCOCOOCCODODDOOODOOCODOOCOCOCOOOOOOOODOOOOOOOOOOOCOCOOOOOO000 


.0000003, 
.0000003, 
.0000003, 
.0000003, 
.0000003, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000002, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000000, 
.0000000, 
.0000000, 
.0000000, 
.0000000, 
.0000000, 
.0000000, 
.0000000, 
.0000000, 
.0000000, 
.0000000, 
.0000000, 
.0000000, 
.0000000, 
.0000000, 
.0000000, 
-0000000; 
.0000000, 
.0000000, 
.0000000, 
.0000000, 


=O S05 O57 
-0.8544087, 
-0.8394190, 
-0.8244294, 
=-0 8034353, 
-0.7944502, 
=O, (794005, 
-0.7644709, 
-0.7494813, 
-0.7344916, 
—O@ tl 2020 
-0.7045124, 
“Oe bo g52acr 
=). 674555 .— 
=('.60 72655" 
-0 .6445539, 
-0.6295643, 
-0.6145747, 
=O oes ur 
-0.5845954, 
=O 0960s s) 
=O, 240lo2, 


moO) Sele V4eis) 


-0.5246369, 
=O25096473, 
-0.4946576, 
=-0.4/796650) 
-0.4646784, 
-0.4496888, 
-0.4346991, 
=0.4197095, 
-0.4047199, 
“0756975057 
-0.3747407, 
=e oor ou, 
-0.3447614, 
=O6 35 2770 lc. 
=O. 2047eZ2t, 
-Og2e5 i925" 
-0.2848029, 


=Qm2o eld ; 


=O42 545230, 
=0.2576540), 
-0.2248444 , 
~G 72096545, 
-0.1948651, 
=O- ioe o°:, 
=O. letee59, 


-0.1498963, 
-0.1349066, 
=02 22291707, 
-0.1049274, 


=O20So9577, 
-0.0749481, 
“07 05995857, 
-0.0449688, 
=O2U2997o72, 
=0.01429696- 
Q0.0000000, 


0 
0 
Q 
0 
0 
0 
0 
0 
0 
0 


0149889, 
O2Z3997e8s, 
0449686, 
U5I9505);, 
0749474, 
Users 7 Z- 
1049271, 
Pio hoo:, 
1349058, 
1498957, 


OVDWDIWDOOOO0O00 


.0000001, 
.0000001, 
.0000001, 
.0000001, 
-0000001, 
-0000001, 
.0000001, 
.0000001, 
.0000001, 
seueoodl, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000001, 
-0000001, 
-0000001, 
-0000001, 
.0000001, 
.0000001, 
.0000001, 
.0000002, 
-0000002, 
.0000002, 
-9000002, 
.0000002, 
.0000002, 
-0000002, 
.0000002, 
.0000002, 
-0000002, 
.0000002, 
-0000002, 
.0000002, 
.0000002, 
-0000002, 
.0000002, 
.0000002, 
.0000002, 
-0000002, 
-0000002, 
-0000002, 
-0000002, 
.0000003, 
-0000003, 
.0000003, 
-0000003, 
-0000003, 
.0000003, 
-0000003, 
-0000003, 
-0000003, 
.0000003, 
-0000003, 
"0900003, 
-0000003, 
0 one 


DODD WDOOOD DOOD DOBAWWDAWBDBDWDWDAWDIAIIADOMDWDWDDWVWWDWDOVOWOWWDVWWMVDWVWVWVVVOVCOCCCOCOO 


‘Oo 
oO 


ics. s910l, 

2224507E+05, 
-248939E+05, 
m7o2s9et+O5, 


w49/329E+05, 
moZLlolE+O5, 
-744631E+05, 
moet oOO7ETO5, 
RoJOZOIIE+OS, 
2111234E+06, 
mass) Oft06., 
»135453E+06, 


.1648855, 
me O1 Sa, 
.1948643, 
.2098541, 
- 2248440, 
eo 96336, 
- 2548227, 
wee J01Z26, 
- 2848024, 
mag g1 923, 
~ 3147821, 
woco77 10, 
- 3447609, 
roe oOd, 
-3747406, 
WoO97 295. 
-4047194, 
-4197092, 
-4346991, 
-4496880, 
-4646778, 
-4796677, 
-4946575, 
.59096464, 
79240363, 
elo ceL 
-5546160, 
.5696049, 
- 5845947, 
.5995846, 
-6145744, 
-6295643, 
-6445532, 


6745329, 
.6895227, 
.7045116, 
.7195015, 
.7344913, 
.7494812, 
.7644701, 
.7794600, 
.7944498, 
8094397, 
8244286, 
8394184, 
8544083, 
8693981, 
8843870, 
.8993769, 
.9143667, 
.9293566, 
.9443464, 
.9593353, 
.9743252, 
.9893150, 
0:0, 90.0 
0.50000 


ODODDODODDDCDOD DODD DVDVDOVDOCOCOCOCOODOD COO CD ODO COOCOOOCOOCOCOOCOCOOCOCOCCOCOOCOOCOOCOOCOCOCO 


OOWDVWIWIOO0000°0O 
OWDWOWO0O00000 


-6595430, | 
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147455E+06, 
15937 584.06, 
171Z208E+06, 
182946E+06, 
194584E+06, 
206116E+06, 
217 S35Enn06, 
s2Z0GS2E706, 
-240004E+06, 
~-251045E+06, 
-261948E+06, 
«2121 06GE06:. 
2635158406, 
-293768E+06, 
.304059E+06, 
.314184E+06, 
~324135E+06, 
-333909E+06, 
~343499E+06, 
~-352900E+06, 
-362108E+06, 
-S/ LEIGH .OG, 
-3/99Z0ET06, 
~388516E+06, 
-396898E+06, 
.405062E+06, 
-413003E+06, 
-420717E+06, 
~-428200E+06, 
.435448E+06, 
-442456E+06, 
-449221E+06, 
-455739E+06, 
-462007E+06, 
-4680Z21E+06, 
-473777E706, 
-479273E+06, 
-484506E+06, 
-489472E+06, 
.494170E+06, 
-498595E+06, 
~502747E+06, 
-506622E+06, 
-510219E+06, 
51355652406. 
~516570E+06, 
.519320E+06, 
, S27 eon r0o., 
-SZ3965E+06. 
-525854E+06, 
~527455E+06, 
-528766E+06, 
-529786E+06, 
-530516E+06, 
»-530954E+06, 
.531099E+06, 
~530954E+06, 
.S30516E206. 
-529786E+06, 
-528766E+06, 
~-527455E+06, 
~525854E+06, 
»523963E+06, 
-521785E+06, 
.519321E+06, 
~516570E+06, 
ol 35Se6E4067 
-510219E+06, 
-506623E+06, 


DOD DODD DOD ODDO DODODODDODO DOOD DODO DOODO DODO DODODOOODOOOODOOOOOOOOOOOOOOOOOOOOOOO0O0000 


DOOD DODD DDODDODODODODODODDDOODDODOODDOOCDOODODOOCOCOOOCOODDOVOOOD OOOO OOO OOOO 9COO9O9O000000 
eleleleleleleleleleleleleleloleloleleleleleleleleleleleleleleleleleleleleleolelelelolelelejeleleloelelelelolelolelelolelololeloyoloyeloleile 
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0.502747E+06, 
QO. 498595E+06, 
O. 494170E+06, 
0.489473E+06, 
0.484506E+06, 
0. 479273E+06, 
0.473777E+06, 
0.468021E+06, 
0. 462007E+06, 
0.455739E+06, 
0.449221E+06, 
0.442456E+06, 
0.435448E+06, 
0. 428200E+06, 
0.420717E+06, 
0.413003E+06, 
0.405062E+06, 
0.396898E+06, 
0.388516E+06, 
0.379920E+06, 
Sell 7E+06, 
0.362108E+06, 
ee >2901E+06, 
0.343499E+06, 
Oes527LOE+06, 
0. 324136E+06, 
0.314184E+06, 
0.304059E+06, 
O22 /OSkt+06, 
0.283315E+06, 
0.272706E+06, 
0.261948E+06, 


.228832E+06, 
~217534E+06, 
.206116E+06, 
~194585E+06, 
.182947E+06, 
~171208E+06, 
mo I575E+06, 
~147455E+06, 
~135454E+06, 
eizsovortO6, 
»111234E+06, 
se 0zZ97e+05, 
Roo?T7UZE+05, 
~-744631E+05, 
moetloLle+Os, 
-497329E+05, 
~373244E+05, 
-248944E+05, 
.124507E+05, 

RESET T 

0,0,0 
Bo), 907,0 
VEM@enr BODY SIDE PLATE 1 


4,0 
Met. 5,-3.2 
eee, 2.0, -3.2 
aaa j= 6, —4.29 
oy ceo, 4625 
PG VEHICLE BODY SIDE PLATE 2 


ee 


j- 


RIGINAL GLOBAL COORDINATE 


dO 
Q 


= al 
PG VEHICLE BODY SIDE PLATE 3 
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